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Prefazione 


Gia da alcuni anni sono stati introdotti, in Italia, presso quegli Atenei con gruppi di 
ricerca attivi nel settore, i corsi di Ingegneria del vento, la cui denominazione è mu- 
tuata in modo diretto da quella della relativa disciplina e cioè lo studio delle azioni 
e degli effetti del vento atmosferico sulle costruzioni e sull’ambiente. Questa nuova 
materia, seppur antichissima nelle origini, si è andata delineando, come avanguardia 
internazionale, negli anni ’60 (soprattutto in Inghilterra, Danimarca, Canada e USA) 
ed è poi entrata in Italia all’inizio degli anni ’80, con la fondazione dell’ ANIV (As- 
sociazione Nazionale di Ingegneria del Vento). 

In realtà tale nuova area di studio nasce decisamente multidisciplinare e spinta 
dal fabbisogno delle società moderne di edificare e pianificare infrastrutture sempre 
più imponenti ed ardite: i materiali sempre più leggeri e le dimensioni sempre più 
imponenti hanno reso molte delle opere dell’uomo decisamente sensibili all’azione 
del vento. Inoltre, a seguito di numerosi e clamorosi insuccessi (crolli o gravi danni, 
spesso poco tempo dopo la costruzione e messa in servizio di un’opera), l’approfon- 
dimento delle tematiche di fluido-dinamica applicata e dell’ aerodinamica dei corpi 
tozzi (prima) e dell’aeroelastica (poi) è divenuto necessario anche nel settore classico 
dell’ingegneria civile e non più solo dell’ingegneria industriale. Infine, i nuovi corsi 
di ingegneria dell’ambiente e gestione del territorio hanno introdotto quegli aspetti di 
fluido-dinamica ambientale tipici dell’interazione del vento atmosferico con l’orogra- 
fia terrestre, ivi compresi i fenomeni di dispersione degli inquinanti e visualizzazioni 
di flusso ai fini di un miglior comfort delle nuove aree edificate. 

Le presenti Lezioni di Ingegneria del Vento, raccolte nell’ambito dell’omonimo 
corso creato ed attivato da pochi anni presso il Corso di laurea specialistica di Inge- 
gneria civile dell’ Università di Firenze, costituiscono uno dei primi tentativi (forse 
il primo in forma così organica e completa) di raccolta delle nozioni di cui sopra, a 
livello di corso universitario; esse pertanto hanno, come obiettivo primario, quello di 
fornire supporto di base agli studenti di corsi di livello avanzato, inclusi i corsi di Ma- 
ster e di Dottorato. In secondo luogo, esse possono altresì offrire un aiuto importante 
a quei professionisti la cui attività prevalente si svolga nella progettazione, gestione e 
sviluppo di opere di grandi dimensioni (sensibili all’azione del vento). Molto spesso 
tecnici (anche ad alto livello) non hanno avuto, nel corso della propria formazione, 
alcun elemento informativo sulla materia e si trovano pertanto ad operare con gravi 
difficoltà, essendo costretti a ricorrere in modo continuo ad esperti e consulenti. 


xviii PREFAZIONE 


I primi capitoli (1, 2 3 e 4) sono dedicati alla descrizione dettagliata del vento 
atmosferico e della sua struttura turbolenta, inquadrando così l’origine delle correnti 
eoliche e le loro caratteristiche pricipali. In generale, le strutture modificano il campo 
di moto delle correnti eoliche incidenti, cosicché le azioni su di esse indotte dipen- 
dono sensibilmente dalle loro prestazioni aerodinamiche. Nel capitolo 5 vengono 
introdotti i concetti basilari dell’aerodinamica dei corpi di forma qualsiasi (tra cui le 
strutture), più propriamente detta aerodinamica dei corpi tozzi. La natura turbolenta 
del vento impone delle analisi di tipo statistico, pertanto nel capitolo 6 vengono ri- 
chiamate le teorie delle variabili aleatorie e dei processi stocastici. La simulazione 
numerica delle storie temporali della velocità eolica è discussa nel capitolo 7. I capi- 
toli 8 e 9 richiamano i concetti fondamentali per la valutazione della risposta dinamica 
delle strutture rispettivamente in campo deterministico e aleatorio. Tali concetti so- 
no propedeutici alla valutazione della risposta strutturale along-wind e across-wind 
trattate rispettivamente nei capitoli 10 e 12. Le prescrizioni dell’ Eurocodice 1, per 
le suddette tipologie di risposte strutturali, sono discusse nei capitoli 11 e 13. La 
risposta aeroelastica delle strutture e le relative prescrizioni dell’ Eurocodice 1 sono 
riportate nei capitoli 14 e 15. Infine, il capitolo 16 prende in esame le tecniche consi- 
gliate per una corretta simulazione in galleria del vento dei fenomeni che interessano 
le strutture reali soggette al vento atmosferico. 

Gli Autori desiderano ringraziare 1 numerosi colleghi che hanno contribuito, nel 
modo più vario, allo sviluppo del materiale grezzo alla base di queste lezioni ed hanno 
fornito utilissime indicazioni sui singoli aspetti; poiché il loro numero è decisamente 
elevato, il ringraziamento viene esteso pertanto a tutta la comunità scientifica dell’ 
ingegneria del vento in Italia, nella quale gli Autori si riconoscono ed alla quale sono 
orgogliosi di appartenere. 


Claudio Borri e Stefano Pastò 


Introduzione 


L’azione del vento ha una notevole importanza nella progettazione delle strutture. 
Un’adeguata progettazione di strutture resistenti alle azioni eoliche si articola secon- 
do gli anelli della cosiddetta catena di Davenport che ne introdusse il concetto. Il 
collasso di un’intera struttura o di parte di essa è legato al mancato adempimento 
delle varie fasi progettuali (anelli) della catena. Pertanto una struttura è tanto più 
resistente all’azione del vento quanto più affidabili sono i vari anelli, ciascuno dei 
quali coinvolge parametri aleatori da inquadrare secondo metodi statistici. Vediamo 
in dettaglio i vari anelli della catena: 


e Vento globale 
e Vento locale 
e Risposta aerodinamica della struttura 


e Risposta meccanica della struttura 


Criteri di progetto 


1.1 Vento globale 


Il vento varia in base alle aree geografiche, tuttavia le velocità eoliche di riferimento 
si ottengono mediando le registrazioni su un periodo non minore di dieci minuti e che 
non superi un’ora. 


2 LEZIONI DI INGEGNERIA DEL VENTO 


1.2 Vento locale 


La rugosita del terreno esercita una grande influenza sul vento. La velocita eolica 
media diminuisce per effetto del terreno, ma allo stesso tempo il vento diviene tur- 
bolento e puo essere descritto soltanto in termini statistici. In particolare, il vento 
atmosferico è un processo stocastico multicorrelato nello spazio e nel tempo. La 
velocità media del vento, per effetto della frizione esercitata dal terreno, aumenta 
con l’altezza secondo un profilo esponenziale o logaritmico. La tipologia del terreno 
influisce su tale profilo in maniera così complessa che risulta difficile quantificarne 
l’effetto con esattezza. Pertanto, si introducono le cosiddette lunghezze di rugosità. 


1.3 Risposta aerodinamica della struttura 


Il carico eolico dipende dalla forma della struttura. La valutazione del comportamen- 
to aerodinamico viene generalmente effettuata in galleria del vento, dove è possibile 
misurare come la velocità eolica si trasformi in un campo di pressioni, multidimen- 
sionale e multicorrelato, agente sulla strutura. 


1.4 Risposta meccanica della struttura 


Noto il campo di pressione intorno alla struttura è possibile valutarne la risposta mec- 
canica con analisi dinamiche, quasi-statiche, o statiche equivalenti. 

Le strutture possono comportarsi in diversi modi quando sono investite dal vento. 
In particolare, le strutture moderatamente rigide possono vibrare nella direzione del 
vento (risposta along-wind) per effetto della turbolenza. 

Le strutture snelle sono principalmente vulnerabili alle vibrazioni nella direzio- 
ne ortogonale al flusso (risposta across-wind) indotte dal fenomeno del distacco dei 
vortici. Inoltre, entro un certo intervallo di velocità (intervallo di sincronizzazione) 
le strutture possono andare in risonanza. 

I ponti sospesi possono presentare oscillazioni che portano al collasso quando il 
modo di vibrare verticale e quello torsionale si accoppiano. Il fenomeno è noto come 
flutter. 


1.5 Criteridi progetto 


Nota la risposta meccanica della struttura è quindi possibile procedere al dimensio- 
namento dei vari elementi strutturali. 


Circolazione atmosferica 


La circolazione dell’aria rispetto alla superficie della terra è indotta dall’irraggia- 
mento della atmosfera terrestre. In particolare, il vento si genera per una differenza 
di pressione tra due punti aventi la stessa quota. Tale differenza si genera per effetto 
di complessi fenomeni termodinamici e meccanici che avvengono nell’ atmosfera in 
maniera non uniforme nello spazio e nel tempo. Vediamo in dettaglio i vari aspetti 
che intervengono nella circolazione atmosferica. 


2.1 Termodinamica dell’atmosfera 


2.1.1 Temperatura dell’atmosfera 


Un semplice esempio (dovuto ad Humphreys [9]) per comprendere l’influenza della 
temperatura sulla circolazione terrestre è dato dal modello semplificativo in Fig. 2.1 
in cui si hanno due colonne di fluido, A e B, riempite fino al livello a. Inizialmente le 
due colonne sono chiuse ed la temperatura all’interno di esse è uguale. Se si aumen- 
ta la temperatura nella colonna A, mantenendo costante quella in B, il fluido in A si 
espande fino a raggiungere il livello b. L'espansione non comporta alcuna variazione 
del peso del fluido, e quindi la pressione in c rimane invariata. Se si apre la valvola 2 
non vi è deflusso tra le due colonne, l'apertura della valvola / comporta uno sposta- 
mento di fluido da A verso B, per effetto del dislivello b-a dovuto al riscaldamento di 
A.Inc la temperatura decresce mentre nella colonna B la temperatura sale. Pertanto, 
aprendo la valvola 2 si genera una circolazione di flusso da B verso A. La circolazione 
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Figura 2.1: Circolazione dovuta alla differenza di temperatura tra due colonne di 
fluido 


si mantiene fin tanto che vi é differenza di temperatura tra le due colonne. 

La terra è inclinata di 66°30’ rispetto al piano della sua orbita, pertanto la re- 
gione equatoriale è maggiormente irraggiata, e dunque riscaldata, rispetto a quelle 
polari. Analogamente al modello di Hemphreys, pur trascurando la complessità del 
fenomeno e i vari fattori che vi rientrano, si comprende l’origine del vento come il- 
lustrato in Fig. 2.2. 

La temperatura dell’ atmosfera è determinata da i seguenti processi: 


e Radiazione solare e terrestre 
e Radiazione nell’ atmosfera 
e Compressione ed espansione dell’aria 


e Condensazione ed evaporazione del vapor d’acqua 


2.1.2 Radiazione solare e terrestre 


L'energia necessaria per spostare le masse d’aria è fornita dal Sole, nella forma di 
calore irradiato. Tuttavia, il Sole non è l’unica fonte, semmai la fonte primaria. In- 
fatti, un ruolo fondamentale è svolto dalla superficie terrestre. In particolare, essendo 
l’atmosfera quasi trasparente all’irraggiamento solare, la quasi totalità dell’ energia 
raggiunge la Terra. Si può assumere che tale energia sia completamente assorbita dal- 
la Terra, che riscaldandosi la rilascia sotto forma di radiazione terrestre. La lunghezza 
d’onda della radiazione terrestre è però molto lunga rispetto a quella solare (dell’or- 
dine di 104), dunque soltanto una piccola parte riesce ad attraversare |’ atmosfera, 
mentre la maggior parte viene assorbita e di nuovo riflessa verso la Terra. 
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Aria Fredda 
Polo Nord 


Aria Calda 


Equatore 


Figura 2.2: Modello semplificato della circolazione atmosferica 


2.1.3 Radiazione nell’atmosfera 


Il calore emesso dalla superficie terrestre è assorbito dagli strati d’aria nelle sue im- 
mediate vicinanze (lo stesso accade sulla superficie degli oceani) che a loro volta 
lo radiano verso il basso e verso l’alto. Gli strati sovrastanti assorbono il calore e lo 
emettono ancora verso il basso e verso l’alto. Il processo si ripete in tutta l’ atmosfera. 


2.1.4 Compressione ed espansione dell’aria 


La pressione atmosferica è prodotta dal peso dell’aria sovrastante. Una piccola mas- 
sa d’aria secca che si muove verticalmente nell’atmosfera, subisce una variazione di 
pressione corrispondente ad una variazione di temperatura. Per valutare questa va- 
riazione di temperatura è necessario ricorrere all’equazione di stato dei gas perfetti 
e alla prima legge della termodinamica, che sono dati rispettivamente dalle seguenti 
relazioni: 


pv = RT (2.1) 
dq = cydT + pdv (2.2) 


dove p è la pressione, v il volume specifico, R è una costante che dipende dal tipo di 

gas (in questo caso aria secca), T è la temperatura assoluta, dq è la quantita di calore 

trasferita alla massa d’aria in moto, e c, è il calore specifico a volume costante. 
Differenziando la prima e sostituendo la quantità pdv si ottiene: 


dq = (cy + R)dT — vdp (2.3) 
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Confrontando la 2.3 con la prima legge della termodinamica dq = cpdT nel caso 
di processo a pressione costante si ottiene facilmente che c, + R = cp. Utilizzando 
ancora l’equazione di stato possiamo scrivere: 


d 
dq = cydT — RT (2.4) 
p 


Per processi adiabatici dq = 0, pertanto dalla precedente relazione si perviene 
alla seguente: 


=~ =0 (2.5) 


che integrata fornisce l’equazione: 


T R 

p \° 
a 2.6 
To (2) Co) 


nota come equazione di Poisson o equazione adiabatica. L’aria secca presenta R/cp = 
0.288. 

Se nell’atmosfera il moto verticale di una particella d’aria è sufficientemente ra- 
pido, lo scambio di calore della particella stessa e l’ambiente circostante può essere 
trascurato e l'ipotesi dg = 0 è accettabile. Pertanto, dall’equazione di Poisson si 
evince che una massa d’aria ascendente subisce una diminuzione di pressione cui 
corrisponde una diminuzione di temperatura. In particolare, nell’atmosfera terrestre 
la diminuzione di temperatura con l’altezza è di circa 1C°/100m e prende il nome di 
rapporto di diminuzione adiabatica, indicato con Va. 

Consideriamo una particella d’aria ascendente, inizialmente alla quota hı carat- 
terizzata da una temperatura T}. Alla quota ho la temperatura dell’aria è 79, mentre 
la temperatura della particella è T = Ti — (ha — h1 )Ya. Poiché la pressione dell’aria 
alla quota hə e quella della particella devono essere uguali, si ha dall’equazione di 
stato che alla differenza di temperatura T% < T> corrisponde una differenza di den- 
sità tra le due. Questo comporta la nascita di una forza che per Tə < TZ sospinge 
la particella verso l’alto, per Tz > 73 sospinge la particella verso il basso facendo- 
la tornare nella sua posizione iniziale. Nel primo caso si dice che la stratificazione 
dell’atmosfera è instabile, nel secondo è stabile. Infine, se Tz = T la stratificazione 
dell’atmosfera si dice neutra. 


2.1.5 Condensazione ed evaporazione del vapor d’acqua 


Secondo la legge di Dalton, la pressione dell’aria umida è data dalla differenza di 
pressione p e la pressione del vapor d’acqua, ossia p — e. Sperimentalmente si osser- 
va che se e supera un certo valore E, noto come pressione di saturazione del vapore, 
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si ha la condensazione del vapor d’acqua in eccesso. 
Una particella ascendente d’aria umida non satura subisce una diminuzione di tem- 
peratura cui corrisponde una diminuzione di pressione di saturazione. Nel caso in 
cui e = E, si ha la condensazione del vapor d’acqua. In questo processo si libera il 
calore latente di vaporizzazione che contribuisce all’ ascesa della particella fornendo 
ulteriore energia oltre a quella interna, che garantisce il movimento verticale dell’aria 
secca. 

Poiché viene liberato il calore latente di vaporizzazione, la diminuzione della 
temperatura della particella d’aria umida in ascesa è minore rispetto al caso d’aria 
secca, pertanto aiuta la convezione dell’aria verso livelli più alti dell’atmosfera. 


2.2 Idrodinamica dell’atmosfera 


Il moto di una massa elementare d’aria è determinato dalla seconda legge di Newton, 
X F = ma, dove m è la massa, a è l’accelerazione, ed X` F è la somma delle forze 
agenti sulla massa stessa. Le singole forze sono le seguenti: 


e la forza dovuta al gradiente di pressione orizzontale 
e la forza dovuta alla rotazione terrestre 
e la forza dovuta all’ attrito 


Descriviamo nel dettaglio tali forze. 


2.2.1 Forza dovuta al gradiente di pressione orizzontale 


Consideriamo un elemento infinitesimo di volume d’aria da - dy - dz, soggetto alle 
forze pe p+(0p/0z)dz agenti rispettivamente sulla faccia inferiore e superiore (Fig 
2.3). In assenza di altre forze, la forza per unità di volume agente sull’elemento sarà 
-—0p/0z. Analogamente, lungo le direzioni x e y, si avrà — 0p/0x e —Op/Oy ri- 
spettivamente. La risultante di queste forze prende il nome di gradiente di pressione 
orizzontale e si denota con —Op/On, con n normale ad una certa isobara. 

Il gradiente di pressione orizzontale è la forza che genera il moto dell’aria. La 
forza per unità di massa esercitata dal gradiente di pressione orizontale —(1/p)0p/On 
prende il nome di gradiente di forza di pressione e sospinge l’aria verso regioni aven- 
ti pressione minore. 
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dy 
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Figura 2.3: Pressione verticale su un elemento di massa d’aria 


2.2.2 Forza dovuta alla rotazione della terra 


La deviazione di una particella in moto da una traiettoria rettilinea è dovuta al mo- 
to della terra. In particolare, la traiettoria diviene curva per effetto di una forza 
apparente, nota come forza di Coriolis, data dalla seguente relazione: 


Fe = 2m(v A w) (2.7) 


dove m è la massa della particella, v è la velocità della particella relativa ad un sistema 
di coordinate solidale alla rotazione terrestre, w è la velocità angolare della terra. La 
forza di Coriolis risulta normale a v e w, pertanto il suo modulo sarà: 


|F.| = 2m|v||w| sin(d) (2.8) 


dove ¢ è l’angolo tra ve w. 

Sia N il polo Nord, e consideriamo un elemento d’aria che si muove lungo un 
meridiano NP. Mentre la massa d’aria si sposta muovendosi da Nord, lo stesso me- 
ridiano risulta spostato verso Est per effetto della rotazione terrestre (NP’). Pertanto, 
nell’emisfero Nord, la massa d’aria che si muove verso Nord è avvertita come un 
vento proveniente da Ovest, mentre l’aria diretta verso Sud vira verso Ovest, dunque 
è percepita come un vento proveniente da Est. Nell’emisfero Sud, si osservano feno- 
meni opposti (Fig. 2.4). 

Introducendo il parametro di Coriolis 


f = 2wsin(¢) (2.9) 
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Parallelo 


Figura 2.4: Moto apparente di una massa d’aria dovuto al moto della Terra 


dove ¢ è la latitudine del punto considerato, segue che la forza di Coriolis per unità di 
massa agente su un elemento d’aria che si muove su un piano parallelo alla superficie 
della terra (Fig. 2.5) è data da: 


F,=mfv (2.10) 


dove v è la velocità della massa d’aria relativa alla Terra, mentre f è un parametro 
tabulato in funzione della latitudine. 


2.2.3 Vento geostrofico 


Ad una quota sufficientemente alta, il vento non risente degli effetti della frizione 
con il terreno ed il moto dell’aria si organizza per filetti fluidi paralleli. Mantenendo 
costante l’accelerazione, il moto è garantito dall’equilibrio tra la forza di gradiente di 
pressione, la forza di Coriolis, e la forza centrifuga. 

L'effetto delle forze agenti su una massa elementare d’aria è mostrato in Fig. 
2.6. Se la particella si muove nella direzione della forza di gradiente (Fig. 2.6a) di 
pressione, P, viene deviata dalla forza di Coriolis, Fa. La sua traiettoria è dunque 
determinata dalla risultante delle suddette forze (direzione JI in Fig. 2.6b). Detta 
Fæ la forza che devia la particella, diretta normalmente alla sua traiettoria, la massa 
d’aria elementare si sposta lungo la direzione ZM data dalla risultante di Fe e P. 
Raggiunto l’equilibrio, la particella si sposterà lungo le isobare (Fig. 2.6c). 

Nel caso più generale in cui le isobare sono curve, dovrà essere tenuta in conto 
anche la forza centrifuga. 
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Figura 2.5: Componenti del vettore di rotazione w 
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Figura 2.6: Equilibrio del vento geostrofico nell’emisfero Nord 
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La velocità in corrispondenza della quale si ha l’equilibrio tra la forza di gradiente 
orizzontale e la forza di Coriolis, prende il nome di velocità geostrofica del vento e 
può essere espressa come segue: 


-yp _ NOD 
2wG sin(¢) = oon (2.11) 
da cui si ottiene: 
_ 1 Op 
G oF ôn (2.12) 


dove f è il parametro di Coriolis, e p è la densità dell’aria. 

Se le isobare sono curve, la particella è soggetta ad un ulteriore forza, quella cen- 
trifuga, G, e per l'equilibrio la sua traiettoria sarà comunque diretta come le isobare. 
Nella condizione di equilibrio, la velocità prende il nome di velocità di gradiente del 
vento che coincide con la velocità geostrofica se la curvatura r — oo: 


v2 1 Op 
V ci A 
orf E r pon C49) 


dove, ipotizzando che la massa d’aria sia nell’emisfero Nord, il segno positivo si ri- 
ferisce ad una circolazione ciclonica (ovvero intorno ad un centro a bassa pressione), 
ed il segno negativo ad una circolazione anticiclonica (ossia intorno ad un centro ad 
alta pressione). In generale, per raggi di curvatura finiti, si ha: 


LIA r Op | rf 3 

Var z T \ On H (4) (2.14) 
e EF rf\* rôp 

Vor = 4 5 \ (2) zan (2.15) 


rispettivamente per venti ciclonici e anticiclonici. Per gli stessi valori di f, re Op/On 
i venti anticiclonici sono più deboli di quelli ciclonici. 


2.2.4 Effetti dell’attrito 


La superficie terrestre esercita sull’aria in movimento una forza di resistenza che 
ne ritarda il moto. Tale effetto diminuisce con l’altezza rispetto al terreno fino a 
raggiungere una quota ò, nota come altezza dello strato limite dell’atmosfera, oltre la 
quale la velocità del vento segue le isobare (velocità di gradiente). L’ atmosfera al di 
sopra dello strato limite prende il mome di atmosfera libera. 

All’interno dello stato limite la velocità del vento non è diretta secondo le isobare. 
Dato un vento geostrofico, consideriamo due particelle che A e B che si muovono 


12 LEZIONI DI INGEGNERIA DEL VENTO 


orizzontalmente nello strato limite, ma a quote differenti. Se la particella A è posta 
ad una quota più alta rispetto a B allora la sua velocità, nonché la forza di Coriolis che 
agisce su di essa, saranno maggiori rispetto a quelle della particella B. L'angolo di 
deviazione tra la direzione del vento e quella delle isobare diminuisce all’aumentare 
della quota, mentre è massimo in prossimità del terreno. Pertanto, l’ andamento della 
velocità del vento con l’altezza diviene una spirale (Fig. 2.7). 


Atmospheric 
boundary layer 


Figura 2.7: Andamento a spirale della velocità del vento all’ interno dello strato limite 
atmosferico 


2.3 Moti atmosferici 
La maggior parte dei processi atmosferici coinvolgono: 
e Velocità del vento orizzontale e verticale 


Pressione 


Temperatura 
e Densità 
e Umidità 

Tali quantità sono governate dalle seguenti equazioni: 
e L'equazione di stato 


e La prima legge della termodinamica 
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e L’equazione di conservazione della massa 
e L'equazione di conservazione dell’ umidità 
e L’equilibrio del moto verticale e orizzontale 


Il moto dell’atmosfera può essere descritto come la sovrapposizione di flussi interdi- 
pendenti caratterizzati da diverse scale che variano da pochi millimetri a chilometri. 
Per analizzare questi moti è conveniente suddividerli in tre gruppi principali: 


e Microscala 
e Mesoscala (o Scala Convettiva) 


e Macroscala (o Scala Sinottica) 


4 
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Figura 2.8: Suddivisione dei moti atmosferici secondo dimensione e durata 


La scala sinottica include moti atmosferici con una dimensione caratteristica di al- 
meno 500km ed una durata di almeno 2 giorni. La microscala comprende que moti 
che hanno una dimensione che non supera i 20km ed una durata minore di un’ ora. 
In quest’ ultima categoria rientra anche la turbolenza atmosferica. Infine tutti i moti 
atmosferici con dimensioni e durata intermedie alle precedenti scale rientrano nella 
mesoscala (Fig. 2.8). 


2.3.1 Circolazione generale 


L'effetto combinato della rotazione terrestre e della frizione prodotta dalla superficie 
terrestre, rende la circolazione atmosferica più complessa rispetto al modello intro- 
dotto in Fig. 2.2. In particolare, si può pensare di schematizzare la circolazione 
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atmosferica in tre cellule (Fig. 2.9). Tale schema è coerente con l’esistenza di un 
fronte di bassa pressione intorno al polo (50°-60° di latitudine Nord) e di un fronte di 
alta pressione nella zona subtropicale (circa 30° di latitudine), che si formano in cor- 
rispondenza della superficie di incontro di venti che provengono da direzioni diverse. 

Questo modello risulta ulteriormente complicato dagli effetti stagionali, come la 
variazione della posizione e dell’intensità dei fronti di pressione e geografici, la diffe- 
renza delle proprietà fisiche del terreno, e la distribuzione non uniforme dei mari sul 
globo. I grandi continenti presenti nell’emisfero boreale provocano deviazioni spes- 
so importanti rispetto a tale modello. Complicazioni analoghe interessano l’emisfero 
australe, che ha però un comportamento meno variabile rispetto a quello boreale. 
Generalmente, la circolazione che interessa l’emisfero australe e quella che interessa 
l emisfero boreale possono essere considerate separatamente. 
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Figura 2.9: Modello di circolazione tricellulare 


2.3.2 Circolazione secondaria di natura termica 


Alla circolazione generale si aggiungono delle circolazioni secondarie, indotte dal- 
la differenza di temperatura. In particolare, queste circolazioni si generano quando 
1 centri di alta o bassa pressione vengono creati per effetto di un riscaldamento o 
raffreddamento dei livelli più bassi dell'atmosfera. Descriviamo in dettaglio i venti 
appartenenti a questa categoria. 


Monsoni 


Durante l’estate la superficie degli oceani si scalda meno di quella dei continenti. 
Pertanto, in corrispondenza dei continenti l’aria più calda tende a salire. Per la con- 
servazione della massa, si attiva dunque una circolazione dagli oceani verso i conti- 
nenti. Durante l’inverno il fenomeno procede al contrario. 

Vista la notevole estensione del continente asiatico, i Monsoni sono venti carat- 
teristici di queste zone. 
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Cicloni tropicali 


I cicloni tropicali, hanno origine sugli oceani tropicali dove la temperatura del’ acqua 
supera i 26°C. L'energia di questi venti deriva dal rilascio di calore latente. In parti- 
colare, per certe condizioni climatologiche, l’aria calda e umida forma una nuvola in 
cui il vapore condensa e rilascia calore latente. L'energia liberata è molto elevata, e i 
venti che ne conseguono possono essere disastrosi. Generalmente, si formano ad una 
latitudine tra i 5°-20°, ed il loro diametro raggiunge le centinaia di chilometri. 


500 200 250 25 200 500km 
e 


“The eye” 


Figura 2.10: Suddivisione dei moti atmosferici secondo dimensione e durata 


In Fig. 2.10 è schematizzata la struttura di un ciclone tropicale. Nella regio- 
ne I risiede /’occhio del ciclone ossia una regione caratterizzata da aria quasi secca 
relativamente calma. Lungo il contorno dell’occhio l’aria sale lentamente, e succes- 
sivamente cade in corrispondenza del centro. Nella regione II l’aria calda e umida 
sale, il vapore condensa rilasciando una grande quantità di calore latente (e dunque 
di energia) e causando piogge molto intense. Nelle regioni III e IV l’aria viene risuc- 
chiata dall’occhio lungo la superficie del terreno, ed incanalata nella regione II dove 
è sospinta ad un’altezza di circa 10km rispetto all’occhio del ciclone stesso. Inoltre, 
in corrispondenza della superfice oceanica l’aria calda e umida, in moto nella regione 
V defluisce nella regione II fornendo ulteriore energia all’uragano. In corrispondenza 
della terra ferma, la frizione con il terreno produce una perdita di energia, pertanto il 
ciclone tende gradualmente ad estinguersi. 

Nel Pacifico occidentale vengono chiamati tifoni, nell’ Oceano Indiano, nel Mar 
Arabico e nel Golfo di Bengala cicloni. A seconda dell’intensità del vento, i cicloni 
vengono classificati come depressione tropicale se il vento non raggiunge i 63km/h, 
tempesta tropicale se non si superano i 118km/h, uragano qualora il vento soffi oltre 
i 118km/h. I danni provocati dai cicloni sia in termini di vite umane sia di danni ma- 
teriali sono molto ingenti. Oltre al vento che con la sua forza distrugge tutto ciò che 
incontra sul suo cammino, i cicloni provocano spesso, nelle zone costiere, inondazio- 
ni generate dall’innalzamento del livello del mare, mentre nelle zone interne possono 
provocare piogge così violente da causare lo straripamento dei corsi d’acqua. Pur- 
troppo non è possibile prevedere la genesi di un ciclone ma, una volta avvistato, se 
ne può seguire l'andamento e il percorso. 
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Figura 2.11: Direzione degli uragani 


2.3.3 Cicloni extratropicali 


Questi venti sono generati o per l’azione meccanica prodotta da ostacoli (montagne 
e barriere) su correnti atmosferiche di larga scala, o per l’interazione di masse d’arie 
in corrispondenza di diversi fronti. 

Le zone di transizione tra diverse masse d’aria, prendono il nome di zone fron- 
tali. La variazione delle propietà fisiche dell’atmosfera in corrispondenza delle zone 
frontali è rapida ed avviene lungo le cosiddette superfici frontali che sono superfici 
di discontinuità. L’intersezione di una superficie frontale con una superfice avente la 
stessa quota prende il nome di fronte. 

Un fronte può essere freddo o caldo in base alla direzione del suo moto, verso 
arie più calde o più fredde. Generalmente, un fronte caldo si sposta lentamente e non 
è associato ad avverse condizioni meteorologiche, mentre un fronte freddo si muove 
più rapidamente e può causare maltempo. 

Le masse d’aria messe in moto dai cicloni extratropicali, hanno propietà presso- 
ché costanti lungo distanze orizzontali comparabili alle dimensioni degli oceani e dei 
continenti. Inoltre, possono essere distinte in tre diverse tipologie in base alla loro 
origine: 


e masse d’aria artiche 
e masse d’aria polari 
e masse d’aria tropicali 


Ciascuna di queste tre tipologie di masse d’aria può essere suddivisa ulteriormente in 
circolazioni continentali e marittime. Le prime sono generalmente fredde e secche, 
mentre le secondo sono calde e umide. 
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2.3.4 Venti locali 
Fohen e Bora 


L'influenza di correnti di piccola scala sulla circolazione generale è trascurabile, tut- 
tavia la loro intensità può essere talvolta considerevole cosicché non se ne può pre- 
scindere durante la progettazione di struttura resistente all’azione del vento. 

Venti forti e localizzati possono generarsi per effetto della variazione di altezza 
del terreno. Se una massa d’aria si sposta da una zona pianeggiante verso una regione 
montuosa o collinare, l’aria si raffredda durante |’ ascesa. Inizialmente la temperatura 
si abbassa di 1°C per ogni 100m di altitudine. Inoltre, quando la temperatura scen- 
de sotto un certo livello, il vapore d’acqua condensa portando piogge o neve, e la 
variazione di temperatura con l’altezza si assesta sul valore di 1°C per ogni 100m 
di altezza. Oltrepassato il picco della montagna, tale valore ritorna ad essere quello 
iniziale di 1°C/100m, e l’aria spostata dalla corrente atmosferica è secca. Se la massa 
d’aria si riscalda durante la discesa si ha il vento chiamato fohen, nel caso in cui la 
massa d’aria non riesce a riscaldarsi a sufficienza si ha la bora. 
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Figura 2.12: Genesi del vento fohen 


Uragani 


Un uragano si può formare in corrispondenza di un veloce moto ascensionale di aria 
calda e umida. La condensazione dell’aria che sale determina un rilascio di energia, 
sotto forma di calore latente, mentre la pioggia prodotta determina un moto del siste- 
ma centrale verso il basso. La scala dell’uragano è di circa 10 km, mentre la durata è 
dell’ordine 1 ora. 


Tornado 


Un tornado è un vortice che ruota intorno ad un asse inclinato, con un diametro di 
circa 300m, e si muove ad una velocità di 10-30m/s. La sua massima velocità tan- 
genziale raggiunge i 100m/s. Al di sotto del tornado si crea una forte depressione che 
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Figura 2.13: Sezione trasversale di un uragano (a sinistra) e vista di un tornado (a 
destra) 


puo agire anche in modo molto violento sugli edifici circostanti. 


2.4 Spettro della componente longitudinale della velocità 
del vento 


Un grafico che riassume i vari fenomeni coinvolti nella circolazione atmosferica è 
dato dallo spettro di Van der Hoven (Fig. 2.14). Lo spettro risale al 1957 e fu inizial- 
mente ottenuto dalla componente longitudinale della velocità del vento ad una quota 
di 100m. Il suo utilizzo è stato successivamente esteso alla descrizione della circola- 
zione atmosferica di qualsiasi sito, e supposto valido per tutto lo sviluppo in altezza 
dello strato limite. 
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Figura 2.14: Spettro della componente longitudinale della velocità del vento 
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Figura 2.15: Suddivisione della componente longitudinale della velocita del vento 
come somma del valor medio, Hu, e della componente turbolenta, u 


Nella zona caratterizzata dalle basse frequenze, si hanno i fenomeni macrome- 
tereologici. Tale regione dello spettro presenta due picchi, uno corrispondente alla 
periodicità del vento giornaliera (brezze caratterizzate da un periodo pari a 12 ore), 
l’altro relativo al normale periodo di sviluppo di una burrasca o tempesta, ossia circa 
4 giorni (100 ore). Lo spettro è stato esteso anche a sinistra oltre il campo di defini- 
zione dello spettro originale, e presenta un ulteriore massimo in corrispondenza della 
periodicità annuale. Alle alte frequenze è possibile osservare un ulteriore picco, in- 
torno a fenomeni della durata di 1-2 minuti, da attribuire alla turbolenza atmosferica. 
Quest'ultima non influenza la circolazione atmosferica, ma invece è importante nelle 
pratiche progettuali. 

Lo spettro è una misura della varianza statistica del vento turbolento. Nella zona 
centrale del grafico in Fig. 2.14, la varianza risulta minima e pressoché costante in 
un periodo di tempo compreso tra 10 minuti ed un’ora. Tale periodo di tempo prende 
il nome di gap spettrale e fornisce un’utile informazione per la valutazione della ve- 
locità di riferimento di un determinato sito. Essendo minima e costante la varianza, 
è minima l’escursione della componente fluttuante della velocità eolica rispetto alla 
sua media che pertanto può essere calcolata con buona approssimazione mediando le 
registrazioni su un periodo compreso tra 10 minuti ed un’ora. 

Poiché la velocità media è stazionaria all’interno del gap spettrale, è possibile 
considerare la componente fluttuante longitudinale del vento come somma del valor 
medio ottenuto su un periodo di 10-60 minuti e della componente fluttuante di origine 
turbolenta (Fig. 2.15). 


Strato limite atmosferico 


La superficie terrestre esercita sull’aria in movimento prossima ad essa una forza 
d’attrito che ne ritarda il flusso. Inoltre, un’ ulteriore frizione si genera tra strati 
d’aria adiacenti che si muovono a velocita differenti. L’effetto combinato di queste 
frizioni rende il flusso turbolento, cosicché la direzione della velocita del vento si 
discosta dalle isobare (turbolenza atmosferica). La regione dell’atmosfera in cui il 
vento risente di questi fenomeni prende il nome di strato limite atmosferico. Lo 
spessore di tale regione si estende per un’altezza che varia da poche centinaia di 
metri a qualche chilometro, in base all’ intensità del vento, alla rugosita del terreno, 
e all’angolo di latitudine. All’interno dello strato limite atmosferico la velocità varia 
con l’altezza, passando da valori più bassi in prossimità della superficie terrestre, 
dove si risente maggiormente dell’attrito da essa esercitata, a velocità gradualmente 
maggiori fino ad arrivare alla velocità di gradiente, che segna il confine tra lo strato 
limite atmosferico e la cosiddetta atmosfera libera, dove la direzione delle masse 
d’aria in movimento segue le isobare. 


3.1 Equazioni del moto dell’atmosfera 


Il moto dell’atmosfera è governato dall’equazioni fondamentali della meccanica dei 
continui, che includono l’equazione di continuità, come conseguenza del principio di 
conservazione della massa, e la seconda legge di Newton. 

Mediando l’equazioni di continuità e di equilibrio, il moto medio dell’atmosfera 
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Figura 3.1: Variazione del profilo della velocita eolica in base alla rugosita del terreno 


puo essere descritto dalle seguenti relazioni: 
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dove U, V e W sono rispettivamente le velocità eoliche nelle direzioni x, y e z di un 
sistema Cartesiano con z asse verticale. I parametri, p, p, f e g sono rispettivamente 
la pressione media, la densità dell’aria, il parametro di Corioli, e l’accelerazione di 
Coriolis. Infine 7, e 7, sono gli sforzi di taglio nelle direzioni x e y. L’ asse x è scelto 
per convenienza coincidente con la direzione dello sforzo di taglio sul piano, 7, (Fig. 
3.2). 

Differenziando la 3.3 rispetto ad x e y si può dimostrare che la variazione ver- 
ticale del gradiente di pressione orizzontale è funzione del gradiente di densità oriz- 
zontale. Trascurando tale variazione (flussi barotropici) si ha che all’interno dello 
strato limite atmosferico il gradiente di pressione orizzontale rimane invariato: 
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Figura 3.2: Assi coordinati 


dove V,, è la velocità di gradiente, r è il raggio di curvatura delle isobare, e n è la 
direzione del vento di gradiente. Se si assume un vento geostrofico si ottiene: 


10p _ 

LF fV; (3.6) 
1 Op 

222 = fu (3.7) 
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dove U, e V} sono le componenti della velocità geostrofica lungo gli assi coordinati 
xey. 

Le condizioni al contorno posso essere scelte ponendo che sulla superficie la 
velocità eolica è nulla, mentre ad un’altezza pari allo spessore dello strato limite 
atmosferico, gli sforzi di taglio svaniscono dal momento che la direzione del vento si 
allinea alle isobare (velocità di gradiente). 


3.2 Soluzione dell’equazioni di moto medio 


Per risolvere il problema del moto medio è necessario introdurre alcune relazioni che 
definiscono gli sforzi di taglio Tu € Ty. 


Tu = pK (x,y, z) 2 


ae (3.8) 
av 
Ty = pK (a, y, 2) 57 (3.9) 
aUN? avVN? 
K(x,y,z) = L?(a,y, 2) (7) + (=) (3.10) 


dove K e L sono rispettivamente la viscosita turbolenta e la lunghezza di mescola- 
mento. 


24 LEZIONI DI INGEGNERIA DEL VENTO 


3.3 Soluzione del moto turbolento 


Dall’equazioni di equilibrio del moto medio è possibile ottenere la seguente espres- 


sione: 
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dove il sovrasegno indica la media rispetto al tempo, e q = \/(u? + v? + w?) è la ri- 
sultante della velocita fluttuante le cui componenti lungo x, y e z sono rispettivamente 
u, ve w. L'equazione 3.11 è nota come equazione dell’energia cinetica della turbo- 
lenza ed esprime l’equilibrio dell’avvezione (primo termine), produzione (secondo 
termine), diffusione (terzo termine) e dissipazione (ultimo termine) dell’energia ci- 
netica. 

L'equazione 3.11 risulta essere vantaggiosa nello studio dei flussi nello strato li- 
mite atmosferico inteso come tridimensionale. In particolare, si hanno le seguenti 
relazioni: 


VTZ +72 = pag? (3.12) 
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dove a; = 0.16, ô è lo spessore dello strato limite atmosferico, e Qe è la velocità di 
gradiente nel caso di vento atmosferico. 


3.4 Profilo della velocità media 


Le tempeste di vento di larga-scala attraversano vaste superfici con differente rugosi- 
tà. Limitando il dominio di integrazione delle equazioni che governano il moto ad una 
regione relativamente piccola rispetto alle dimensioni della tempesta, e con rugosità 
approssimativamente uniforme, è possibile assumere che la velocità media nel piano 
x € y sia costante lungo tutto lo sviluppo della regione: 0U/0x = 0 e 0V/0y = 0. 
In questo caso il flusso si dice orizontalmente omogeneo. Se poi si tiene conto del 
fatto che il gradiente di pressione orizzontale ostacola la crescita dello strato limite 
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atmosferico, allora la componente verticale della velocità eolica può essere trascurata 
e l'ipotesi di flusso orizzontale omogeneo risulta ancora valida. Tenendo conto delle 
precedenti assunzioni, nell’ulteriore ipotesi di vento, geostrofico sono valide le 3.6 e 
3.7, e quindi l’equazioni 3.1 e 3.2 diventano: 


1 OT, 
V -V= At OE (3.16) 
197, 
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Figura 3.3: Spirale di Ekman 


3.4.1 La spirale di Ekman 


Se, in aggiunta a quanto appena detto, consideriamo le condizioni al contorno per gli 
sforzi di taglio (Eq.ni 3.8 e 3.9) e assumiamo che la viscosità turbolenta (Eq. 3.10) 
sia costante, il modello che ne deriva prende il nome di spirale di Ekman in cui le 
relazioni 3.16 e 3.17 sono equazioni differenziali a coefficienti costanti. Imponendo 
le condizioni al contorno U = V = 0 per z = 0 e U = Uy, V = V; per z + œ, la 
soluzione del sistema è la seguente: 


1 

U = —G|1-e “(cos(az) — sin(az 3.18 
30 [L €%(cos(az) — sin(a) 6.18) 
1 
V2 

cona = y f /2K. 
Tali equazioni descrivono la spirale di Ekman (Fig. 3.3) ma non forniscono sem- 
pre buoni risultati a causa dell’ipotesi di viscosità turbolenta costante con l'altezza. 


V = —=G |1 — e “ (cos(az) + sin(az)| (3.19) 
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Figura 3.4: Suddivisione dello strato limite in due regioni: strato superficiale e strato 
esterno 


3.4.2 Lo strato limite turbolento di Ekman 


Assumendo dei modelli in cui la viscosità turbolenta varia con l'altezza, si hanno 
risultati più plausibili. Appartengono a questa categoria di modelli quelli che ipotiz- 
zano lo strato limite atmosferico diviso in due regioni: lo strato superficiale e lo strato 
esterno. È plausibile che lo sforzo di taglio 7, dipenda dalla velocità del vento in pros- 
simita del suolo, dalla rugosità del terreno, zo, e dalla densità dell’aria p. Pertanto To 
dovrà essere espresso in funzione delle precedenti quantita: To = F(Ui+ Vj, z, Zo, p) 
dove i e j sono i versori diretti rispettivamente secondo gli assi coordinati x, y. La 
precedente relazione può essere convenientemente espressa in forma adimensionale : 


ian (=) (3.20) 


Ux o 
nota come legge del muro che descrive il flusso nel strato supericiale. La quantita 


tie = ec (3.21) 
p 


prende il nome di velocità di taglio, mentre fı è una funzione di 2/20. 
Anolagamente, nello strato esterno si può assumere che la variazione di veloci- 
tà [(Ugi + Vj) — (Ui + Vj)], alla generica quota z, sia espressa in funzione di To, 


CAPITOLO 3. STRATO LIMITE ATMOSFERICO 27 


dell’ altezza dello strato limite ô, e della densità dell’ aria p: 


Ui+ Vj Ugit Voi fo (£) 


(3.22) 


Ux Ux 


Esprimendo la 3.20 nella forma equivalente: 


nea iù (3) (*)| (3.23) 


uguagliando i primi membri della 3.22 e della 3.23 e scegliendo opportunamente le 
funzioni fı e fo: 


fil€) = (nev) i (3.24) 
Aa (me) i ži (3.25) 
(3.26) 


con k e B costanti, si perviene alle seguenti relazioni: 


Ug = La In o: (3.27) 
Ux ke Zo 
V; B 
si (3.28) 
Ux k 
da cui si ottiene la velocita geostrofica: 
9 Ô \ Us 
G= B2+Inî—|)- (3.29) 
Zope k 


Inoltre si può dimostrare che lo spessore dello strato limite atmosferico, 6 sia dato 
dalla seguente relazione: 


Ux 
6 =c— (3.30) 
f 
dove c è una costante. 
Infine, è dunque possibile definire il profilo della velocità media del vento gene- 
ralizzando la 3.27 per una velocità ad una quota generica z: 


U(2) = Lu, In (=) (3.31) 


oO 


dove k = 0.4 è la costante di von Karman, zo è la rugosita del suolo, e ux = \/To/p 
è la velocità d’ attrito. 
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Figura 3.5: Lunghezza di rugosità zo 


Tabella 3.1: Categorie del terreno e relativi parametri 


Categoria del terreno kp zo (mM) 2min(m) 
I Mare aperto, laghi con minimo 5 km di 0.17 | 0.01 2 
estensione sopravento e campagna senza 
ostacoli 
II Terreno agricolo con recinzioni, piccole fattorie | 0.19 | 0.05 4 
isolate, case o alberi 
III Aree suburbane o industriali o foreste 0.22 | 0.3 8 
permanenti 
IV Area urbana in cui almeno il 15% della 0.24 1 16 
superficie sia coperta da edifici la cui altezza 
media superi i 15m 


3.4.3 Lunghezza di rugosità 


La lunghezza di rugosità, zo può essere interpretata come la dimensione caratteristica 
dei vortici che si generano per effetto della frizione dell’aria in moto e il suolo (Fig. 
3.5). 

Nell’Eurocodice 1 sono state introdotte quattro categorie di terreno (Tabella 3.1), 
per le quali è definito un fattore di terreno ky proporzionale alla velocità d’attrito, 
e che aumenta con la rugosità del terreno. Inoltre, l’Eurocodice 1 introduce anche 
un'altezza minima entro la quale la velocità rimane costante. 

In particolare, l Eurocodice 1 fornisce un profilo di velocità logartimico: 


U(z) = Urefkrln (=) È S dra (3.32) 


Zo 


dove U,e; è la velocità di riferimento del sito in esame. 


Gli elementi del suolo che contribuiscono ad aumentarne la scabrezza sono detti 
elementi di rugosità. Tali elementi aumentano la resistenza al vento e pertanto con- 
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15 20 25 30 35 40 45 


Um (m/s) 


Figura 3.6: Profili di velocita media del vento secondo le quattro categorie previste 
dall’Eurocodice 1 


tribuiscono a produrre ulteriore turbolenza. 
Se gli elementi di rugosità sono uniformemente distribuiti sul suolo, è possibile 
stimare la lungezza di rugosità zo con la seguente formula empirica: 


Zo = 0.5h 2 (3.33) 
Ar 

dove h è l’altezza degli elementi di rugosità, A, è l’area degli elementi ortogonali 
alla direzione del vento, e A; è l’area d’influenza di ciascun elemento. Se A, e A; 
sono dello stesso ordine di grandezza (elementi di rugosità molto fitti), il flusso viene 
traslato verso l’alto di una quantità pari all’ altezza media degli elementi di rugosita. 
Tale fenomeno lo si registra ad esempio nelle foreste (Fig. 3.7). 

La traslazione della superficie di contatto viene tenuta in conto introducendo 
l'altezza media della rugosità fitta, d, nell’espressione del profilo logaritmico della 
velocità media del vento: 


(3.34) 


3.5 Terreni non omogenei 


Nel caso in cui il vento spiri in prossimità di un suolo, caratterizzato da una rugosità 
che nel complesso può essere intesa come uniforme, il profilo di velocità lo si può 
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U(z) 


Figura 3.7: Profilo di velocita media in corrispondenza di una foresta 


approssimare con una legge logaritmica. 

Consideriamo il caso in Fig. 3.8 in cui il vento lambisce dapprima un suolo 
con rugosità zo] e successivamente un suolo con 252. In corrispondenza del primo 
terreno, e ad un’altezza hz dopo la transizione, la velocità del vento è determinata 
soltanto dalla rugosità z01. La rugosità 202 influisce sulla velocità del vento solo 
entro l’altezza h1, mentre nella regione intermedia tra hı e ho si ha una graduale 
transizione in cui entrambe le rugosità influenzano la velocità del vento. La regione 
cha si estende dal suolo fino ad h2 prende il nome di strato limite interno, mentre al 
disotto di hı si ha lo strato di equilibrio. La lunghezza dello strato limite esterno è 
stato valutato da Elliot che formulò la seguente relazione: 


0.8 
hax) = 202 [ozs + 0.03 In (4)| (=) (3.35) 


202 202 


dove la distanza x è misurata, in direzione del vento, dal punto in cui cambia la 
rugosità. Si evince dalla 3.35 che hg aumenta più rapidamente dopo una transizione 
verso terreni più rugosi, piuttosto che dopo una transizione verso terreni meno rugosi 
(Fig. 3.9). 


3.5.1 Topografia 


La presenza di una discontinuità nel piano di campagna (Fig. 3.10) determina un’ac- 
celerazione del flusso, e conseguentemente un aumento della pressione eolica sulle 
strutture. Per questo motivo l’Eurocodice 1 introduce un coefficiente di topogra- 
fia, cr, che dipende dalla velocità del vento alla base della discontinuità, dalla sua 
pendenza, dalla sua lunghezza misurata nella direzione del vento, dalla sua altezza 
rispetto al suolo, e dalla distanza della struttura dal suo picco. 
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Figura 3.8: Profilo di velocità media in corrispondenza di terreni con diversa rugosità 
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Figura 3.9: hə dopo una transizione di rugosità 
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separation 


separation separation 


Figura 3.10: Discontinuità del terreno 
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Tabella 3.2: Tipici valori di ô e a per i profili esponenziali 


Area costiera Terreno aperto Terreno suburbano Centro città 
a dgr (m) | a dgr (m) a dgr (m) a dgr (m) 
- a .16 275 0.28 400 0.40 520 
1/10 215 1/7 275 1/4.5 400 1/3 460 


5(Z) 


U(Zz1,Z51) 


Figura 3.11: Profili di vento per differenti rugosita 


3.6 Profilo esponenziale della velocita del vento 


Il primo modello di profilo di velocita del vento ad essere introdoto fu quello espo- 
nenziale: 


a 
Zgl 
U(zg1) = U(2g2) (4) (3.36) 
Z92 
dove a è un coefficiente che dipende dal tipo di terreno, mentre 241 € 242 sono due 
quote rispetto al terreno. Ponendo 242 = Zref € zgi = Zz si ottiene quello che 
originariamente era il profilo di velocità media previsto dalle norme canadesi: 
a 
z 
U(z) =U (ees) ( ) (3.37) 
Zref 


E bene ricordare che il profilo esponenziale non ha un fondamento fisico ma 
fornisce valori accettabili. 


34 LEZIONI DI INGEGNERIA DEL VENTO 


Se si assume che la quota di riferimento sia l’altezza dello strato limite, zref = 
Ò = dgr, allora la legge esponenziale è funzione solo di a: 


zZ 


u(z) =G (+) (3.38) 
dar 


3.6.1 Discontinuita del terreno 


Consideriamo, come nel caso precedente, due terreni aventi differenti lunghezza di 
rugosità, siano queste Zo € zo] CON Zo1 < Zo. La differenza di rugosità induce due 
differenti profili di velocità (Fig. 3.11). 

Utilizzando il profilo esponenziale, per trovare la velocità U(z,, Zo) nota la velo- 
cità U (291, 201), si ottiene: 


a(Zo) 5 a(Zo1) 
U (zg, 20) = U(2g1; 201) | | = (3.39) 


Turbolenza atmosferica 


Il vento nello strato limite atmosferico è sempre turbolento. La turbolenza, è dovuta 
ad un treno di vortici che si genera per effetto dell’attrito dell’aria con il suolo e con 
tutti gli ostacoli presenti su di esso (elementi di rugosità). Quanto più grande è il 
valore della rugosità equivalente del sito in esame, tanto più turbolento è il vento in 
prossimità della superficie di rugosità. L'effetto si ripercuote su tutto lo strato limite 
atmosferico fino a raggiungere l'atmosfera libera dove la velocità del vento segue la 
direzione delle isobare. Vi è dunque un graduale attenuamento della turbolenza al- 
ľ aumentare della quota rispetto al suolo, dovuto alla minore influenza della rugosità 
superficiale sulla circolazione atmosferica. 

Il vento, in quanto turbolento, può essere inteso come un processo stocastico in 
cui ciascun vortice contribuisce a caratterizzarne l’energia (densità spettrale), in base 
alle proprie dimensioni e al proprio periodo (o alla propria frequenza). Generalmente, 
nei vortici più grandi prevalgono le le forze fluide inerziali, mentre quelli più piccoli 
sono caratterizzati dalle forze fluide viscose. Quest’ultimi tendono a dissipare mag- 
giore energia, che viene di volta in volta rifornita dai vortici più grandi. In definitiva, 
nel vento turbolento vi è un continuo trasferimento di energia dai vortici più grandi a 
quelli più piccoli. 

La valutazione della risposta strutturale diviene così molto complessa. Tuttavia 
alcune evidenze sperimentali permettono di semplificare in parte il problema. Nel ca- 
pitolo precedente, abbiamo visto che lo spettro della componente longitudinale della 
velocità del vento è caratterizzata dal gap spettrale che garantisce la stazionarietà 
del valor medio della velocità eolica su un periodo di tempo compreso tra dieci mi- 
nuti ed un’ora. È dunque possibile trattare separatamente la componente fluttuante 
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turbolenta del vento e quella media che di per sé induce una risposta strutturale di 
tipo statico. La risposta al vento turbolento viene valutata in campo dinamico ed è 
intesa come un processo aleatorio da studiare con opportuni modelli statistici. Un 
ulteriore constatazione facilita la soluzione del problema dinamico: la turbolenza at- 
mosferica può essere intesa come un processo il cui momento statistico del secondo 
ordine è stazionario. Inoltre, il vento turbolento è, con buona approssimazione, di tipo 
Gaussiano, ossia i primi due momenti statistici permettono di caratterizzarlo comple- 
tamente. Pertanto, essendo stazionari media e varianza del processo Gaussiano ad 
essa associato, la turbolenza può essere intesa anche come ergodica. L’ assunzione 
di ergodicità del processo sfoltisce notevolmente le pratiche sperimentali in sito o in 
galleria del vento: una singola registrazione del vento turbolento risulta essere rap- 
presentativa di tutto il processo stocastico ad esso associato. Se venisse meno questa 
proprietà, una corretta campagna sperimentale dovrebbe essere condotta su una serie 
numerosa di registrazioni. 

Riassumendo, il vento atmosferico può essere inteso come un processo Gaussiano 
stazionario ergodico. Tralasciando la componente media, una completa descrizione 
fisica della turbolenza è data dalla caratterizzazione dell’intensità di turbolenza, delle 
sue scale integrali, della sua densità spettrale di potenza e della sua coerenza che ne 
identifica la correlazione spaziale. 


4.1 Intensità di turbolenza 


Siano U, V e W le velocità eoliche medie nelle direzioni x, y e z di un sistema Car- 
tesiano con z asse verticale (Fig. 3.2), e u, v e w le rispettive componenti turbolen- 
te. Nell’ipotesi di terreno omogeneo e flusso orizontale, le deviazioni standard delle 
componenti fluttuanti, Cu, Cy € Cw, possono essere considerate costanti sui piani indi- 
viduati da differenti quote rispetto al suolo. Inoltre, è stato osservato che tali quantità 
decrescono molto lentamente entro le altezze tipiche delle strutture convenzionali, 
o addirittura fino ad un’altezza pari alla metà dello strato limite atmosferico, come 
postulò Armit. Fino ad un’altezza di 100-200 metri è possibile assumere che: 


o, = Au, (4.1) 
Oy © 0.750, (4.2) 
Cw X 0.50, (4.3) 


con A = 2.5 per zo = 0.05m e A ~ 1.8 per zo = 0.3m. 
A partire dalle deviazioni standard, di definiscono le intensità di turbolenza in 
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funzione della quota z: 


I,(z) = do (4.4) 
Le sa (4.5) 
Iu() = DO) (4.6) 


Fino ad altezze pari a 100-200 metri, si può assumere che le componenti fluttuan- 
ti della velocità del vento siano distribuite normalmente (distribuzione Gaussiana), e 
abbiano media nulla!. Tuttavia, qualora la deviazione standard di tali processi stoca- 
stici cade al di fuori di un certo intervallo, sia esso +3, l’ipotesi di Gaussianità può 
essere a sfavore di sicurezza. 


4.2 Scale temporali 


Nell’assunzione di processo stazionario, la funzione di autocorrelazione permette di 
stimare le informazioni che una registrazione della velocità del vento, ad un certo 
istante tı, fornisce sulla valore della stessa registrazione al temp to. Indicando con 
T = tı — t2 il periodo di tempo che intercorre tra due valori della stessa registrazione, 
è possibile definire la funzione di autocorrelazione nel seguente modo: 


Ru(z,7) = E [u(2,y,2,t1)u(x,y,2,t1+7)] (4.7) 


che è funzione della sola altezza z nell’ipotesi di flusso orizzontalmente omoge- 
neo che investe un sito uniformemente rugoso. La funzione di autocorrelazione può 
essere anche adimensionalizzata: 


E [u(a, y, 2 ty ules y, z, ti aS T)] = E lula, y, z, tıJulz, y, z, tı + T)] 
E[u? (x,y, 2, t)] o? (z) 


Pulz, T) = 
(4.8) 


dove o,,(z) è la varianza. La funzione di autocorrelazione, p,,(z,7), assume così 
valori compresi tra 0 ed 1. 
La scala temporale T(z) è definita come: 


+00 
T(z) = f Pulz, T)dT (4.9) 
0 


ed è un tempo di memoria: una misura di u al tempo tı fornisce diverse informa- 
zioni circa u(t; + T) se T << T(z), al contrario se 7 >> T(z). Un esempio di 


‘La media viene considerata a parte. 
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funzione di autocorrelazione é data dalla Fig. 4.1, che fa riferimento ad una buona 
approssimazione é data dalla seguente relazione: 


palT) = e TIT) (4.10) 


Figura 4.1: Funzione di autocorrelazione pu(T) 


4.3 Scale integrali 


Le scale integrali sono una misura della dimensione caratteristica dei vortici che co- 
stituiscono la corrente eolica, o in maniera equivalente, della dimensione media del- 
la raffica di vento. Considerando la componente u, le rispettive scale integrali si 
scrivono: 


+00 

Li =} Pulz, Tx) dre (4.11) 
0 
+00 

L? a Pulz, Ty)dry (4.12) 
0 
+00 

Li =} Pulz, rz)drz (4.13) 
0 


dove pu(z,rj) j = x,y,z è la funzione di cross-correlazione tra due punti della 
raffica in direzione j distanziati di r;. In generale è possibile definire 3 scale integrali 
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per ciascuna componente di velocità. 

Secondo l’ipotesi di Taylor, la turbolenza atmosferica può essere schematizzata 
come un treno di vortici congelati che si muovono alla velocità media del vento. In 
questa ipotesi, considerando j = x, si ha rx = U(z)r e quindi: 


L*(z) = U(z)T(z) (4.14) 


valida per tutte le 9 scale integrali di turbolenza. 
Le scali integrali dipendono dall’altezza, dalla rugosita del suolo, e anche dalla 
velocita del vento. Per altezze comprese trai 10-240 metri vale la seguente relazione: 


Dre) (4.15) 
dove C e m dipendono dalla rugosita del sito, zo (Fig. 4.2). 


Le rimanenti scale integrali sono spesso espresse in funzione della scala integrale 


1000 


0.001 0.01 0.1 1.0 10 


zo (meters) 


Figura 4.2: Valori di C e m in funzione di zo 


longitudinale, L2 (z): 


LY ~ 0.3L% (4.16) 
IZ x 0.2L" (4.17) 


4.4 Spettro della componente longitudinale della velocità 
del vento 


La turbolenza può essere intesa come la sovrapposizione di vortici caratterizzati da un 
periodo e dunque da una frequenza angolare w = 27n. L'energia totale del processo 
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è quindi la somma dei vari contribuiti energetici associati ad ogni singolo vortice. 
Dall’equazione di moto del flusso turbolento si evince che i termini inerziali sono 
associati ad un trasferimento di energia dai vortici più grandi a quelli più piccoli, 
mentre i termini viscosi si riferiscono ad una dissipazione di energia che interessa 
principalmente i vortici più piccoli. A causa dei questa dissipazione, la turbolenza 
tenderebbe ad estinguersi, tuttavia la corrente viene rifornita di energia dai vortici con 
dimensioni maggiori. Il flusso turbolento è dunque garantito da un equilibrio energe- 
tico tra l’energia introdotta nella corrente dai vortici più grandi, e l’energia dissipata 
da quelli più piccoli. 


Prima ipotesi di Kolmogorov 


Il moto dei vortici minori è governato dalla quantità di energia fornita, che per l’equi- 
librio è pari a quella dissipata, e dalla viscosità. Da tale ipotesi, deriva che i piccoli 
vortici sono indipendenti dalle condizioni al contorno, a tal punto che questi non 
hanno una direzione preferenziale, ossia costituiscono un flusso isotropo. 


Seconda ipotesi di Kolmogorov 


I vortici maggiormente responsabili della dissipazione energetica sono i vortici molto 
piccoli aventi cioè delle lunghezze d’ onda molto piccole. Conseguenza di quest’ipo- 
tesi è che i vortici più grandi, aventi lunghezze d’onda maggiori, si muovono indipen- 
dentemente dalla viscosità, e pertanto il loro moto è determinato soltanto dalla quan- 
tità di energia fornita, che per l'equilibrio sarà uguale a quella dissipata. La seconda 
ipotesi di Kolmogorov è valida nel cosiddetto sotto-intervallo inerziale e fornisce 
l’espressione generale della componente longitudinale della velocità del vento: 


Suul K) = a K 5/3 (4.18) 


dove a = 0.5, e è la quantita d’energia trasferita (o dissipata), e K = 27r/A è il 
numero d’onda con A lunghezza d’onda. 


4.4.1 Spettro nel sotto-intervallo inerziale 


Si può assumere che l’energia trasferita, e, sia data dalla seguente relazione: 


neo (4.19) 
p Oz 
con U (z) dato dal profilo di velocità media: 
1 
UG) =k (4.20) 


k zo 
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Si dimostra, utilizzando le precedenti relazioni, che lo spettro della componente 
longitudinale della velocità del vento può essere scritta nella forma: 


nSuu(z,n) nz yore 
= 0.26 4.21 
> T(z) (an 
dove la quantità f = nz/U (z) è nota come coordinata di Monin. 
Si riportano di seguito alcuni spettri presenti in letteratura: 
Spettro di Von Karman 
nSuu(z,n) nLî 1 
=4 4.22 
u2 U(z) penal’? we) 
h +70.8 (74) | 
Spettro di Davenport 
2 
12 1 
nin) a o 0) E (4.23) 
ta 10 14 (1200n 2 
Uio 
Spettro di Kaimal 
us (1 + 50f)5/8 
Spettro utilizzato dall’Eurocodice 1 
mSuu(2 0) _ gg Mr : (4.25) 
u2 U(z) nLz 13/3 
E + 10.274 5 


Per completezza si forniscono anche gli spettri di letteratura relativi alla compo- 
nenti trasversali della velocità del vento. 


Spettro di Kaimal 


NSyy(Z,N) _ 15 nz 1 (4.26) 


Ms CA E + 9.575] tà 
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= = = Von Karman 
Ti Kaimal 
EC1 

cl bh pa a Davenport 


=i 0 
10 L2/U(2) 10 10 10 


Figura 4.3: Confronto tra vari spettri della componente longitudinale della velocità 
del vento 


Spettro di Lumley e Panofsky 


NS ww(Z, n) nz 1 
= 3.36 - (4.27) 
ug U(z) Le 
1+10 (z5) 


4.5 Correlazione spaziale della turbolenza 


Il cross-spettro normalizzato, Sy(P;, Pa, n) della turbolenza longitudinale, fornisce 
la dipendenza statistica in funzione della frequenza n, tra due punti nello spazio P; e 
Ps. Lo si esprime come segue: 


— Suul Pi, Pan) 
V Suul Pi, n) Suul Po n) 


Sn(P1, Po, n) (4.28) 


dove Suu( P1, Pa, n) è il cross-spettro della componente turbolenta nei due punti con- 
siderati. Suu( Pj, n) j = 1,2 è lo spettro della j-esima componente turbolenta. 

La parte reale dello cross-spettro prende il nome di co-spettro, Yuu, ed è una fun- 
zione che decade all’aumentare della distanza. Questo decadimento dipende dalla 
dimensione dei vortici che generalmente può essere stimata come il rapporto tra la 
velocità e la frequenza, U/n. La parte immaginaria prende il nome di quad-spettro. 
Inoltre |Sy(P,, Pa, n)| = VCuu(Pi, Pa, n) dove Cuu(Pi, Pa, n) è la funzione di 
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coerenza. Il co-spettro e la radice quadra della coerenza coincidono quando la parte 
immaginaria del cross-spettro normalizzato è nulla o trascurabile. In generale nella 
valutazione della risposta di strutture soggette al vento turbolento, si assume che la 
parte immaginaria di Sy (P1, P., n) sia nulla, pertanto si ottiene: 


Suu(Pi, Pa, n) = VSuu(Pi, Nn) Suu(Pa, n) Cohuu(Pi, Pa, n) (4.29) 


dove Cohuu(Pi, Pa, n) = Cual Pi, Pa, n) = YPuul Pi, Pa, n). 

Date le distanze ry = |yi — y2| e rz = |z1 — 22] lungo y e z tra due punti nello 
strato limite atmosferico, una buona approssimazione della funzione di coerenza è 
data da: 


\/ Cara Py ra n) = Cohuu(Pi, Py, n) = PuulTy, Tz, n) = eu (Cury) (Cara)? 


(4.30) 


dove U = 3 [U(z21) + U(z2)], mentre Cy e C, sono due costanti che generalmen- 
te possono assumersi pari a 10. Tali costanti garantiscono il decadimento della 
correlazione spaziale con la distanza. 


Aerodinamica dei corpi tozzi 


Lo studio dei fenomeni aerodinamici nasce dalle pratiche aeronautiche, e successiva- 
mente fu esteso alle strutture immerse nello strato limite atmosferico. Di particolare 
interesse, nelle pratiche strutturali, è la conoscenza della turbolenza del vento in- 
cidente, oppure la turbolenza prodotta dalle stesse strutture (signature). Poiché, in 
generale, le strutture presentano forme qualsiasi, si fa riferimento all’aerodinamica 
dei corpi tozzi. 

In questo capitolo vengono riportati i concetti basilari dell’aerodinamica dei corpi 
di forma qualsiasi, e quindi delle strutture. 


5.1 Equazioni che governano il moto di una corrente fluida 


Una corrente fluida, come ad esempio il vento atmosferico, si muove rispettando le 
seguenti equazioni di equilibrio. 
5.1.1 Equazioni di moto 


Dato un volume elementare di fluido, dV, il vettore di velocità del fluido può essere 
scritto come segue: 


3 
u= 5 uiii (5.1) 
i=1 


dove u; (i = 1,2,3) sono le componenti della velocità lungo x; (i = 1,2,3) che 
costituiscono un sistema cartesiano di riferimento i cui versori sono rispettivamente 


46 LEZIONI DI INGEGNERIA DEL VENTO 


dati da i; (i = 1, 2,3). 
Le equazioni di moto del volume dV di fluido si scrivono come somma delle 
forze d’inerzia, delle forze di massa, e delle forze di superficie: 


pot = = pgi + 2 i= 1,2,3 (5.2) 


dove: 


Du; do; 
p n dV pgidV 3 md (5.3) 


sono rispettivamente (procedendo da sinistra verso destra) le forze d’inerzia, le forze 
di massa, e le forze di superficie. Nella 5.3, p è la densità del fluido, t è il tempo, g; è 
l’i-esima componente dell’accelerazione, e o;; una tensione agente su una particella 
di fluido (Fig. 5.1). Inoltre, 


3 
Pai Yu; 2 (5.4) 


a )dxgdxs 


de 
o,, 9x, 0x, nt è 


Figura 5.1: Forze agenti su un volume elementare di fluido 


5.1.2 Equazioni di continuità 


Per il principio di conservazione della massa, l aumento di massa fluida contenuta 
all’interno di una certa superficie chiusa deve essere uguale alla differenza tra la mas- 
sa che esce ed entra dal volume racchiuso dalla superficie. L'equazione di continuità 
si scrive come segue: 


3 
(pu) _ Op 
È CA (5.5) 
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Se il fluido è incompressibile, Op/Ot = 0, pertanto l’equazione di continuità 
diviene: 


vi (5.6) 


5.1.3 Equazioni di Navier-Stokes 


Al contrario di un solido, un fluido in condizioni statiche può resistere solo a pressioni 
normali. In condizioni dinamiche invece, il flusso può resistere anche a sforzi di 
taglio. Se tali sforzi sono dovuti alla viscosità, jz, allora si parla di fluidi Newtoniani. 
Ad esempio, lo sforzo di taglio in Fig. 5.2 è espresso come: 


lo) 
bid a (5.7) 


dove il fattore di proporzionalità è la viscosità fluida, u, la cui unità di misura è 


Kq/(m- sec). 
Conviene dividere il tensore delle tensioni, o, in una parte idrostastica, p, ed una 


x, 


Figura 5.2: Velocità che varia linearmente allontanandosi da una superficie 
deviatorica, d: 


p i=j 
i = 5.8 
Pij AES (5.8) 


3 
1 = 
dij = 2u | eij — gli Desk (i,j = 1,2,3) (5.9) 
k=1 
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dove: 


2 Ox; 0 if 


Con la precedente scomposizione si arriva alla seguente relazione: 


i ; he Besa 
€ij = 1 (= È 2u) di; = l : 1 (5.10) 


3 
1 
Cij = —pòij + 2u (a _ gli > aa) (5.11) 
k=1 


che sostituita all’ equazione di moto (Eq. 5.2) porta alla definizione dell’ equazioni di 
Navier-Stokes: 


Du; Op E; 1 3 
aor alla SS x; È (e = DN] i,j =1,2,3 (5.12) 
Jar k=1 


Se la viscosità, u, può essere considerata costante, la 5.12 diviene: 


3 3 
Du; Op Ou; 10 Our, 
data e Z Sa 5.13 
P Dt PIi Ox; = A x Ox? = 3 Ox; di, Ox, ( ) 
j=1 J kal’ 
che in forma vettoriale diviene: 
3 3 
Du Op . du 
Pay = PE- Dig tH as (5.14) 
i=l j=l CJ 


5.1.4 Equazioni di Bernoulli 


Nel caso in cui il fluido è incompressibile (Eq. 5.6) e inviscido (u = 0), l'equazione 
5.14 diviene: 


LPE (5.15) 


se xı e la coordinata lungo la quale avviene il moto del fluido, si ha per integra- 
zione della 5.15 che: 


1 
= |u|? + f= cost. (5.16) 
2 p 


in ogni punto di una linea di corrente fluida. 
La 5.16 è nota come teorema di Bernoulli che nel caso più generale si scrive: 
l 2 
zPY +p = cost. (5.17) 


dove u è la velocità della corrente fluida lungo una linea di corrente. La quantità 
5 pu? & nota come pressione dinamica. 


CAPITOLO 5. AERODINAMICA DEI CORPI TOZZI 49 


5.2 Flusso vorticoso 


Consideriamo una corrente fluida che si muove lungo un percorso curvilineo (Fig. 
5.3). Sia dr la distanza tra due linee di flusso, l'equazione di moto per il flusso sarà 
data da: 


2 
dpdA = pdrdA~ (5.18) 
r 


dove p è la pressione lungo su un elemento elementare di fluido, dA è l’area di tale 
elemento, u? /r è accelerazione del flusso dovuta alla curvatura di raggio, r. 


pdA + d(pdA) dr 
dr 


STREAMLINES 


pdA 


Figura 5.3: Corrente fluida lungo una linea curva 


La 5.18 indica che la variazione di pressione normale alle linee di corrente curvi- 
linee è data dalla relazione: 


gdr 


dp = pu (5.19) 


= 
Nel caso in cui il percorso è circolare si ha il flusso vorticoso. Per il teorema di 
Bernoulli si ottiene: 


d d 
SE (5.20) 
u r 
che integrata fornisce la seguente relazione del vortice libero: 
ur = cost (5.21) 


secondo la quale la velocità diminuisce all’ aumentare del raggio secondo una legge 
iperbolica. Tuttavia tale relazione è vera soltanto nelle zone più esterne del vortice, 
dal momento che in quelle più interne prevalgono gli sforzi di attrito (Fig. 5.4). 
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TRANSITION 


Figura 5.4: Distribuzione di pressione e velocita in un moto vorticoso 


Figura 5.5: Separazione del flusso in corrispondenza di uno spigolo 
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Figura 5.6: Separazione e formazione della scia di un profilo alare 


Figura 5.7: Separazione del flusso di una superficie circolare 
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5.3 Strato limite e separazione 


Le forze fluide si differenziano in forze fluide inerziali e forze fluide viscose. Il 
rapporto tra le due è un indice del tipo di flusso in esame. Tale rapporto è noto come 
numero di Reynolds, Re: 


Roe eee (5.22) 
u V 


dove v è nota come viscosità cinematica e L è una dimensione caratteristica del cor- 
po. Quando Re è grande predominano gli effetti inerziali, altrimenti prevalgono gli 
effetti viscosi. Generalmente, una corrente fluida che investe un ostacolo origina una 
varietà di numeri di Reynolds, a seconda della regione del corpo che si considera. 

La separazione del flusso lungo un ostacolo è dovuta alla presenza di un gradiente 
avverso di pressione, come ad esempio può accadere in corrispondenza dello spigolo 
di un corpo (Figg. 5.5, 5.6 e 5.7). Il processo di separazione dello strato limite su- 
perficiale genera dei vortici che producono forti suzioni nel punto di separazione, e 
che defluiscono a valle del corpo fomando la scia. La scia è una regione costituita 
da vari strati tra i quali vi è uno scambio di energia continuo. La scia può essere 
di tipo turbolento o laminare a seconda se lo scambio di energia avviene al livello 
macroscopico (scia turbolenta) o a livello molecolare (scia laminare). 

In base al numero di Reynolds, si possono distinguere tra flussi laminari o turbo- 
lenti. Molti flussi turbolenti, possono essere intesi come una scia di vortici prodotta 
da una serie di ostacoli. Infatti, la turbolenza atmosferica, come menzionato nei capi- 
toli precedenti, è il prodotto dell’ attrito di un flusso laminare con gli ostacoli presenti 
sul suolo ed intesi come rugosità superficiale equivalente. 


5.3.1 Distacco dei vortici e formazione della scia 


La separazione dello strato limite superficiale e la formazione della scia dipendono 
dalla geometria del corpo e dal numero di Reynolds. 

I profili aerodinamici (profili alari) sono studiati in modo tale da non presentare 
separazione lungo la sezione, ma solo in corrispondenza della coda dove i vortici che 
si generano sono così piccoli che la scia risulta essere molto sottile. Tuttavia, se il 
profilo alare ha un angolo d’attacco pronunciato rispetto alla direzione del flusso, la 
separazione e la formazione della scia vorticosa variano sensibilmente a tal punto da 
compromettere le prestazioni del profilo stesso (Fig. 5.6). 

I corpi aventi superfici curve (Fig. 5.7) presentano punti di separazione e scie 
vorticose che dipendono sensibilmente dal numero di Reynolds. 

In particolare, nel caso di un cilindro circolare (Fig. 5.8), per valori bassi di Re, 
circa Re < 5, le linee di flusso seguono il contorno della sezione (flusso non se- 
parato). Nell’intervallo 5 < Re < 40, il flusso si separa nella regione a valle del 
cilindro formando una coppia di vortici simmetrici. Se il il numero di Reynolds vie- 
ne incrementato ulteriormente, 40 < Re < 150, la scia diviene instabile e i vortici 
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si staccano alternativamente dalla regioni laterali del cilindro formando una scia di 
vortici di tipo laminare (regime laminare). Nell’intervallo 150 < Re < 300 la scia 
comincia a diventare turbolenta anche se lo strato limite superficiale sul cilindro é an- 
cora laminare. I numeri di Reynolds compresi tra 300 < Re < 3 - 10° caratterizzano 
una scia completamente turbolenta (regime subcritico). Il regime critico si manifesta 
nell’intervallo 3- 10° < Re < 3.5 - 10° ed è caratterizzato da una scia turbolenta e 
disorganizzata a tal punto che non è più distinguibile il distacco dei vortici. Infine, 
per Re > 3.5 - 10° si ha il regime supercritico dove riappare nuovamente il distacco 
dei vortici all’interno di una scia turbolenta. 


— pÀ Re < 5 REGIME OF UNSEPARATED FLOW 


E, 
—t Ne 5 TO 15 < Re < 40 A FIXED PAIR OF FOPPL 
NEGATO = VORTICES IN WAKE 


40 < Re < 90 AND 90 < Re < 150 
TWO REGIMES IN WHICH VORTEX 
STREET IS LAMINAR 


150 < Re < 300 TRANSITION RANGE TO TURBU- 
LENCE IN VORTEX 


4 O 300 < Re Z 3X 10° VORTEX STREET IS FULLY 
UL 


3x10° = Re < 35x 105 


LAMINAR BOUNDARY LAYER HAS UNDERGONE 
TURBULENT TRANSITION AND WAKE IS 
NARROWER AND DISORGANIZED 


3.5 X 105 < Re 


RE-ESTABLISHMENT OF TURBU- 
LENT VONTEX STREET 


Figura 5.8: Caratteristiche della scia di un cilindro circolare al variare di Re 
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Figura 5.11: Flusso intorno ad un corpo tozzo: Re = 250 
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Figura 5.12: Flusso intorno ad un corpo tozzo: Re > 1000 
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Figura 5.13: Confronto tra sezioni quadrate e rettangolari 
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I corpi con spigoli vivi (corpi tozzi) presentano dei punti di separazione che di- 
pendono dalla forma del corpo, mentre il processo di formazione della scia dipende 
dal numero di Reynolds. Consideriamo ad esempio il corpo tozzo in Fig. 5.9: per 
valori molto bassi del numero di Reynolds, Re = 0.3, il flusso non si separa ma 
segue il contorno della sezione; per Re = 10 (Fig. 5.10) il flusso si separa in corri- 
spondenza degli spigoli vivi e forma una scia caratterizzata da due vortici simmetrici. 
Aumentando ulteriormente il numero di Reynolds, per Re = 250 (Fig. 5.11) la scia è 
composta da due vortici che si staccano alternativamente dagli spigoli della sezione. 
Infine, per Re > 1000 (Fig. 5.12), le forze fluide inerziali prevalgono pertanto non è 
possibile assistere alla formazione di grandi vortici, anzi la scia è composta da vortici 
di piccole dimensioni. 

Inoltre, analizzando il flusso intorno a cilindri quadrati e rettangolari (Fig. 5.13), 
si può notare un diverso comportamento aerodinamico, ad alti numeri di Reynolds. 
La sezione quadrata induce una separazione seguita da una scia ampia a valle del 
corpo. La sezione rettangolare, più allungata, presenta una separazione in corrispon- 
denza degli spigoli a monte del corpo, seguita da un riattacco del flusso lungo la 
sezione, e successivamente una nuova separazione a valle della sezione. La scia che 
ne consegue è più stretta. Tutto ciò si traduce, come vedremo più avanti, in coeffi- 
cienti di forza differenti tra le due sezioni. 


Figura 5.14: Legge del numero di Strouhal 
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Sezioni trasversali st 


Vento 
02 


per tutti i numeri di Ra 


Vento 0,5<db<10 Figura 5.16 


Cra 
d 
db-1 0,1 
Vento n db=1,5 0,10 
db=2 0,14 


interpolazione lineare 


Vento db=1 0,13 
db=2 0,08 


interpolazione lineare 
db=1 0.16 
Vento = db-2 0,12 
interpolazione lineare 
db=1,3 0.11 
= db=20 007 


TT 
db<5 0,154 
(a) 
Pei $<db<10 1,1 -d/b+1/ 


dbz 10 0,083 


impalcali di pont 


Figura 5.15: Numeri di Strouhal per diverse sezioni (Eurocodice 1) 


5.3.2 Numero di Strouhal 


I vortici si staccano dalle regioni laterali di un corpo con una frequenza, fys, che 
segue la nota legge del numero di Strouhal: 


fous = i (5.23) 


dove D è una dimensione caratteristica della sezione, e S; è il numero di Strouhal che 
dipende da diversi parametri, quali ad esempio la forma della sezione (Fig. 5.15) ed 
il numero di Reynolds. Si può facilmente notare che la legge del numero di Strouhal 
è una legge lineare (Fig. 5.14). 
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Figura 5.16: Numeri di Strouhal per sezioni rettangolari (Eurocodice 1) 


5.3.3 Coefficienti di pressione e di forza 


Un corpo soggetto ad una corrente fluida, è soggetto a delle pressioni locali, p, che 
seguono il teorema di Bernoulli: 


1 
30U° +p= cost (5.24) 


ossia la somma della pressione dinamica calcolata nelle immediate vicinanze del cor- 
po, 5 PU”, e di quella statica, p, costante lungo una linea di flusso. 
Generalmente, si preferisce far riferimento al coefficiente di pressione: 


P — Po 
Cp I o? (5.25) 
dove p — po rappresenta la differenza tra la pressione locale e quella del flusso indi- 
sturbato, po. 

Analogamente, è possibile definire il coefficiente di drag, Ca, definito nella dire- 
zione del flusso, il coefficiente di lift, Cz, definito nella direzione normale a quella 
del flusso, e il coefficiente di momento, C'm: 


Fp 


Cp = ——— (5.26) 
5pU2B 
Fr 
CL = = (5.27) 
5pU2B 
M 
Cy = > (5.28) 
5pU2B? 


dove B è una dimensione caratteristica della sezione, Fp e Fr sono rispettivamente 
le forzanti di Drag e di Lift per unità di lunghezza, e M è il Momento. Alcuni esempi 
sono riportati in Fig. 5.18. 

I coefficienti di pressione e di forza dipendono, in generale, dal numero di Rey- 
nolds. Consideriamo ad esempio la distribuzione dei coefficienti di pressione di un 
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cilindro circolare (Fig. 5.17). Nella scia del cilindro circolare, il coefficiente di pres- 
sione risulta essere costante al variare di Re. Ciò e dovuto al fatto che la separazione 
dello strato limite superficiale è avvenuta a monte del corpo stesso. Tuttavia, all’au- 
mentare di Re il punto di separazione arretra, pertanto la scia si forma più a monte 
e risulta essere più sottile. Ne consegue una minore suzione nelle regioni laterali. 
Integrando le pressioni locali nella direzione del flusso si ottiene la forzante di Drag, 
e quindi il coefficiente di Drag. Un tipico andamento del Cp per un cilindro circo- 
lare è riportato in Fig. 5.19, dove si può notare che nel regime critico Cp diminuice 
drasticamente per effetto dell’arretramento del punto di separazione e dell’assotti- 
gliamento della scia. Nel regime supercritico e in quello transcritico, il coefficiente 
di Drag tende ad aumentare sebbene il suo valore sia comunque minore rispetto al 
regime subcritico. Inoltre, l'andamento di Cp al variare di Re dipende anche dalla 
rugosità superficiale (Fig. 5.20). In particolare, la rugosità superficiale anticipa la 
transizione tra un regime ed un altro. 

Sebbene, per una sezione circolare, il coefficiente di Drag varia al variare di Re, 
lo stesso non si può dire per un cilindro a sezione quadrata, per il quale Cp rimane 
quasi costante al variare del numero di Reynolds (Fig. 5.20). Ciò è dovuto al fatto 
che la separazione è dettata dallo spigolo a monte della sezione, e non vi è possibilità 
di un riattacco del flusso lungo il corpo. Smussando gli spigoli del cilindro quadrato 
(Fig. 5.20), si ottiene, per contro, una notevole dipendenza del coefficiente di Drag 
dal numero di Reynolds, a tal punto che si può addirittura notare una regione critica 
come accade nei cilindri circolari. 


-2.0 


8 (DEGREES) 


Figura 5.17: Distribuzione dei coefficienti di pressione su un cilindro circolare, al 
variare di Re 
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Figura 5.18: Valori di Cp e Cr, per varie sezioni 
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Figura 5.19: Andamento di C'p per un cilindro circolare, al variare di Re 
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Figura 5.20: a, b) Andamento di C'p, al variare di Re e della curvatura degli spigoli 
di una sezione quadrata; c) Variazione di Cp, per un cilindro circolare, al variare di 


Re e della rugosità superficiale ottenuta con della sabbia avente grani di dimensione 
k 


Variabili aleatorie e processi stocastici 


Per molti problemi fisici, le incertezze sono costituite da valori numerici, ossia dai 
valori che una stessa variabile può assumere. La soluzione del problema, pertanto, 
non è deterministica, ma deve essere necessariamente discussa in termini statistici 
(Davenport [3], Newland [11]). 

In questo capitolo sono discusse le variabili aleatorie continue e i processi sto- 
castici che possono essere intesi come un’insieme di variabili aleatorie continue. 


6.1 Variabili aleatorie continue 


Data una variabile aleatoria continua, X, la probabilità che essa sia maggiore o uguale 
ad un certo valore, x, è data dalla funzione di distribuzione di probabilità, Fx (x): 


Fx (x) = P(X < 2) (6.1) 


Generalmente. è preferibile far riferimento alla derivata di Fx(x), che prende il 
nome di funzione di densità di probabilità: 


_ OF x(x) 
px(z) = Əz (6.2) 
da cui chiaramente si ottiene che 
Fx@)=PA<2)= | px (Og (63) 
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con la condizione che 
+00 
f px@de=1 (64) 


Una variabile aleatoria continua è completamente definita quando è definita la sua 
funzione di densità di probabilità, px (x), o alternativamente quando se ne conoscono 
i momenti di ordine n: 


00 (6.5) 


dove il momento del primo ordine E[X] = ux è noto come media della variabile 
aleatoria continua X. 
Analogamente, è si definiscono i momenti centrali di una variabile aleatoria: 


EX] = [- xpx(x)dx = px 
+00 
E|(X — ux)]] = f (a? — wx)px(x)dx = Var[X] = 0% 
Lög (6.6) 


+00 
BUX = ux)9)= fe" - ux)px(e)de 


dove 03, è la varianza di X. 
L’indice di variazione di una variabile aleatoria è definito come: 


pes EE (6.7) 


HX 


dove ox è nota come deviazione standard di X. 
Nella pratica, i valori di una variabile aleatoria sono noti dai risultati di N espe- 
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rimenti. In tal caso le precedenti definizioni dei momenti divengono: 


1 N 
E[X] = 2 = ux 
i=1 


1 N 
E[X?] = a 5x 
i=1 


(6.8) 
iS 
E[X"]= x DISP 
j=l 
mentre i momenti centrali saranno dati dalle seguenti relazioni: 
1 N 
EX] = IX = bx 
si 
1 N 
EIX — ux)?]= 7 DOG — px)? = ox 
i=l (6.9) 


ret 
E[(X — ux)"]= N di — px)" 


6.1.1 Variabili aleatorie Gaussiane 


Una variabile aleatoria Gaussiana, X, è definita dalla seguente funzione di densità di 
probabilità: 


px(x) = ceti ua) (6.10) 


La caratteristica principale delle variabili aleatorie Gaussiane, è che sono com- 
pletamente definite quando sono definiti i primi due momenti o i primi due momenti 
centrali, poiché: 


E(X — px)"] =0 Vp>2 (6.11) 


6.1.2 Variabili aleatorie congiunte 


Nel caso di più variabili aleatorie si definisce la funzione di distribuzione di probabi- 
lita congiunta. Ad esempio, date due variabili aleatorie, X e Y si ha: 


x y 
Fxy(x,y)=P(X <anY <y)= f f pxy (€, n)dêdn (6.12) 
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con la condizione che: 
+00 +00 
ni ni pxy (x, y)dady = 1 (6.13) 
—oo —oo 


È possibile risalire alla funzione di densità di probabilità di una sola variabile 
aleatoria, ad esempio X, nel seguente modo: 


+00 
LE Ji badi (6.14) 


—co 


Se la variabile aleatoria Y assume un valore noto, y, allora si puo definire la 
funzione di densita di probabilita condizionale: 


vale = oc) SACD 


= 6.15 
JIZ pxy (a, y)de py (y) oT 


La dipendenza di due variabili aleatorie è espressa dal coefficiente di correlazio- 
ne: 
_ E(X — ux)(Y - uy)] _ Cov[X,Y] 


PXY (6.16) 
OXxXOY TXOY 


Se pxy = 0 le due variabili aleatorie sono statisticamente indipendenti pertanto 
pxy(« |y) = px(2) (6.17) 
da cui si ottiene: 


pxy(£,y) = px(x)py (y) (6.18) 


Se pxy = 1 le due variabili aleatorie sono statisticamente dipendenti, e dunque la 
densità di probabilità sarà espressa da pxy (x, y). 


6.1.3 Momenti di variabili aleatorie congiunte 


Date due variabili aleatorie congiunte, X ed Y, i loro momenti centrali del primo 
ordine (medie) saranno definiti dalle seguenti relazioni: 


+00 +00 
n] = EX] = px = | f epxy(e,y)dedy (6.19) 


+00 +00 
xi[Y]= E[Y] = uy = / f ypxy (z, y)dzdy (6.20) 


mentre i momenti centrali del secondo ordine definiscono la matrice di covarianza, 
T(X,Y): 


Tey = | E[(X — px)?| uri i 2 Ee ci 
AY [EIN — wy )(X — px)] E((Y — py)? oyx oY 
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(6.21) 


dove oxy = oy x è la covarianza tra X e Y. La definizione dei momenti centrali di 
ordine superiore è analoga a quelli appena visti. Inoltre, tali definizioni, in generale, 
possono essere estese al caso di n variabili aleatorie congiunte. 


6.1.4 Variabili aleatorie Gaussiane congiunte 


Date due variabili aleatorie continue Gaussiane, la densità di probabilità congiunta si 
esprime nel seguente modo: 


1 
pxy (2, y) E ii 
TOXOY 
1 (ux)? , w-uy)? (@-ux) (a-py) 
1 fark ) | a { a 2pxy =x sey }} 
x e xy (6.22) 


y 1 — Py 


dove pxy è il coefficiente di correlazione tra X e Y. 


6.2 Processi stocastici 


Un processo stocastico è una famiglia di variabili aleatorie che dipendono da un pa- 
rametro. Nella pratiche strutturali, tale parametro è il tempo, dal momento che i 
processi stocastici a cui si fa riferimento, sono delle forzanti che variano nel tempo 
(Fig 6.1). 
La struttura di probabilità di un processo stocastico, X (t), è completamente defi- 
nita quando sono definite tutte le funzioni di densità di probabilità ad esso associate: 
Px) (21) 


PX (t1) X (ts (£1, £2) 
(t1) X (t2) (6.23) 


PX (t1)X (t2)...X (tn) (T1 T2, +++, Ln) 


In alternativa, un processo stocastico è completamente definito se sono definiti i 
suoi momenti: 


(6.24) 


EX (t1)X(t2) ...X(tn)] 


oppure i suo momenti centrali 
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Figura 6.1: Esempio di processo stocastico definito nel tempo 


(6.25) 


E{TX (t1) — ux (tr) |[X (te) — ux (t2)] - - - [X (tn) — ux (tn)]} 


dove uy (t) è la funzione di media, Ry x (t1, t2) è la funzione di autocorrelazione, 
mentre 


kol X (t1) X (t2)] = F{[X (t1) — ux (t1) [X (t2) — px(t2)]} = 
Rxx(tı, t2) — ux(tı)ux(t2) (6.26) 


è la funzione di covarianza o cumulante del secondo ordine. Inoltre si può dimostrare 
che il cumulante del primo ordine è uguale alla funzione di media: 


K[X(t)] = E[X(ti)] = px(t;) i=1...n (6.27) 
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Se ciascuna variabile aleatoria, che compone il processo stocastico, è nota come 
risultato di N esperimenti, allora i momenti assumono la seguente forma: 


EIXO = 7 XO) = ax) 
= a 
EIX (t1) X (to)] =D) Xi(t)X;(t2) 
i=1 j=1 (6.28) 
1 N N N 
E[X(t1)X(t2)...X(t)]= 7 DVD Kilts) X;(t2)...Xk(tn) 


(6.29) 


6.2.1 Processi stocastici stazionari 


Un processo stocastico si dice completamente stazionario se i suoi momenti sono 
tempo-invarianti: 


E|X (t1)?] = E[X(t0)P] =... = E[X (tn)? (6.30) 
oppure stazionari di ordine m se 


E[X(t1)P] = E[X(t2)”] =... = E[X(tn)?] Yp< m (6.31) 


6.2.2 Processi stocastici stazionari e Gaussiani 


Un processo stocastico Gaussiano, è un processo stocastico le cui componenti sono 

caratterizzate da un funzione di densità di probabilità di tipo gaussiano (Eq. 6.10). 
Inoltre, dal momento che la densità di probabilità di una variabile aleatoria Gaus- 

siana è completamente definita quando sono note la sua media, wx, e la sua varianza 
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o3, (Eq. 6.10), allora un processo Gaussiano è completamente stazionario se i primi 
due momenti (o momenti centrali) sono stazionari: 


E[X(t)] =x Vt (6.32) 


kol X (t1)X (t2)] = F{[X (t1) — ux (t1)][X (t2) — ux (ta) = 
= Rxx(t1,t2) — ux(tı)ux(t2) = Rxx(T)— u% (6.33) 


In questo caso la funzione di autocorrelazione, Rxx(7) dipende da 7 = tı — t2 e 
non da t e to separatamente. 


Osservazione 


Se il processo è stazionario allora la media può essere sottratta dal processo stocastico 
e considerata a parte. In questo caso, il processo stocastico originario viene scompo- 
sto in una componente media ed una fluttuante (stocastica). Ad esempio, la turbolen- 
za atmosferica è un processo stocastico Gaussiano, che risulata essere stazionario su 
un intervallo di tempo compreso tra 10 minuti ed un’ora. Pertanto, mediando su tale 
intervallo il processo associato alla turbolenza atmosferica, avremo una componente 
media che induce una risposta strutturale statica, ed una componente stocastica che 
induce una risposta dinamica. 


6.2.3 Processi stocastici ergodici 


Un processo stocastico si dice ergodico se i suoi momenti coincidono con quel- 
li associati alla singola registrazione, X(t), intesa come unica variabile aleatoria 
rappresentativa di X (t) : 


T Jo N f 
T 
KR H] = ox = z f (X(t) — uxļ?dt = £ Yt) — px]? (6.34) 


Inoltre, la funzione di autocovarianza di un processo stocastico ergodico si scrive: 


1 /T 
ka[X(t)X(t+7)] = o f [X(t) — wx |[¥(t +7) — ux]dt = 


YO x(t) — wx] [X (Gi + 7) — ux] = Rxx(7) — u% (6.35) 


i 


z| = 
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si può subito notare che per T = 0 si ha K2[X?(t)] = 0% = Rxx(0) — u}. 

L'ipotesi di ergodicità è quasi sempre accettata nella sperimentazione. In partico- 
lare, acquisito un certo processo stocastico, X(t), durante un intervallo di tempo T, 
con una frequenza di campionamento fe, il numero di campioni monitorati risulterà 
N = T fe. Noti X(t), T, ed N si possono trovare le grandezze di interesse riportate 
nella 6.34 e nella 6.35. 

Se il processo ergodico fosse anche Gaussiano, sarebbe completamente descritto 
dai primi due cumulanti «K1 € K2. 


Osservazione 


La funzione di covarianza ha le seguenti propietà: 


ko[X(t)X(t+7)] = Ka(7) < K2(0) (6.36) 
K2(T) = K2(-7) (6.37) 
k2(T) 20 T —> +00 (6.38) 


6.3 Processi stocastici multivariati 


Un processo stocastico multivariato, X(t), è un vettore le cui componenti, X;(t), so- 
no dei processi stocastici. 

Un processo stocastico multivariato è completamente definito quando sono defi- 
niti i suoi momenti centrali o i momenti statistici: in particolare, avremo il vettore 
delle medie 


E[X(t)] 
E[X(t)] 
nO] = px(t) = (6.39) 
E[Xy(6)] 
la matrice di covarianza: 
Kxixi(ti,to)...  Rxixy(t1,t2) 
Txx(t1,t2) = : : = 
KXyX1(t1, ta)... KXyXy(t1, t2) 


= Rxx(ti, t2) — ux(ti)ux(te)’ (6.40) 


e tutti gli altri momenti centrali associati al processo stocastico multivariato X(t). 
La matrice Rxx(t1,t2) è la matrice di correlazione di X(t), ossia il momento 
del secondo ordine di X(t): 
E[X1(t1) X1(t2)] sal E|X1(t1) Xn (t2)] 
: : : (6.41) 


Rxx(t1,t2) = : : 
E|Xn (ti) X1(t2)] su E|Xn(t1) Xn (t2)] 
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6.3.1 Processi multivariati stazionari 


Nel caso stazionario i momenti centrali (ed i momenti statistici) del primo ordine 
non dipendono dalla variabile ¢, mentre tutti gli altri momenti centrali (e i momenti 
statistici) dipendono solo da T = ta — tı: in particolare il vettore delle medie sarà 
dato da 


na [X(8)] = px = vt (6.42) 
E[Xn(t)] 
mentre la matrice di covarianza si scriverà 


Kx X(T). 6X1 X(T) 
Lxx(7) = De 
KXNX1 (7) -+ RXNXN (7) 
= Rxx(7) — uxu% (643) 
dove la matrice di correlazione Rxx (T) si scrive: 
E|Xi(Xi(t+7)] ... EDGMAXN(t+7)] 
Rxx(T) = : : : (6.44) 
E|Xn(t)Xi(t+7)] ... EDXy()XN(t+ 7) 
La stazionarietà si ripercuote chiaramente su tutti i momenti centrali (e tutti i 
momenti stastici) di ordine superiore al secondo. 
Per un processo stocastico multivariato, Gaussiano, stazionario, sarà sufficiente 
conoscere i primi due momenti centrali per descriverlo completamente. Un esempio 


è dato dalla turbolenza atmosferica che generalmente è riconducibile ad un processo 
con le suddette caratteristiche. 


6.4 Densità spettrale 


La densità spettrale di potenza di un processo stocastico stazionario, X(t), è data 
dalla trasformata di Fourier della funzione di correlazione di X (t): 


I! se 
Sxx(w) = = | Rxx(7)e dr (6.45) 


Pertanto, la funzione di correlazione risulta essere la trasformata inversa di Fou- 
rier della densità spettrale, Sx x(w): 


+00 
Rxx(7) = f Sxx(w)e™ dw (6.46) 


—oo 
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Inoltre, noto lo spettro di X (t) si dimostra che 


+00 
ot = Rxx(0) = J Sx x(w)dw (6.47) 


—co 


Figura 6.2: Esempio di densita spettrale di potenza 


Poiché la densità spettrale è una funzione pari, vale la seguente relazione: 
+00 +00 
ox = Rxx(0) = f Sxx(w)dw = Gxx(w)dw (6.48) 
0 


—00 


dove Gx x(w) = 2Sx x(w) è detto spettro unilatero di X(t). 


Figura 6.3: Esempio di densità spettrale bilatera, Sx x(w), e unilatera, G x x (w) 


Nel caso di un processo multivariato stazionario, X(t), si definisce la matrice 
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delle densita spettrali: 


Sx,x,(w) Sx,x,(w) ... Sxxy(w) 

Sxx(w) = 2 : ay : (6.49) 
OX XG (w) SXNX2 (w) yah Okay (w) 

: i n | 
Sxx(w) = a si Rxx(7)e dr (6.50) 
dove 

E|Xi(t)Xi(t+ 7)] (RE E|Xi(t)Xn(t+7)] 

Rxx(r) = i (6.51) 
E|Xn(t)Xi(t+7)] ... EDXy(0)XN(t+ 7) 


è la matrice di correlazione di X(t). 
La matrice di densità spettrali è una matrice Hermitiana, ossia le componenenti 
extra-diagonali sono l’una il complesso coniugato dell’altra: 


Sx;x; (w) = SX;x; (w) iFj (6.52) 


La parte reale dei cross-spettri, Re[Sx,x;(w)], è una funzione pari, mentre la 
parte complessa, Im[Sx,x;(w)], è una funzione dispari. Gli auto-spettri, Sx,x;(w), 
sono funzioni reali e pari. 

Inoltre, è possibile definire la coerenza Cohx,x;(w), che è in generale una fun- 
zione costituita da una parte reale (co-spettro), Rx; x; (w), ed una complessa (quad- 
spettro) Tx,x, (w): 


Cohx,x; (w) = Rx;x; (w) pt Dx (w) (6.53) 


In generale si ha: 


s 
\/|Cohx,x,(w Pa xxl) (6.54) 
XX: 


w) Sx; x, (w) 


Nota la matrice delle densità spettrali di X(t) è possibile risalire alla matrice di 
correlazione ad esso associata: 


2 
ox, ses OX{XN 


Rxx(0) = : ic : (6.55) 


2 
OXnX, «+: TX, 
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dove 
o%, = / Sx,x,(w)dw (6.56) 
OX;X; = I OX OG (w)dw (6.57) 


sono rispettivamente la varianza di X;(t), e la covarianza tra X;(t) e X;(t). Inoltre, 
poiché la parte immaginaria del cross-spettro è una funzione dispari, il suo integrale 
è nullo. Pertanto, la covarianza è un numero reale ed inoltre 


OX;X; = OX;X; (6.58) 


6.5 Esempi di processi stocastici 


Si riportano le funzioni di correlazione e le densita spettrali di alcune tipologie di 
processi stocastici ricorrenti nelle pratiche ingegneristiche: 


6.5.1 Processo a banda stretta 


R | | 
2 SL 0) =2 5, S,() A 


| | 
7 So 


Figura 6.4: Processo a banda stretta 
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6.5.2 Processo a banda larga 


vX) 


l | | | | 


' 


| di 


el 


Figura 6.5: Processo a banda larga 


XY OM" 


a) 


R= 2 nS, Ht) 


S(O) 4 


Figura 6.6: Rumore bianco 


6.5.3 Rumore bianco 


Si definisce rumore bianco un processo stocastico le cui variabili aleatorie sono 
statisticamente indipendenti tra loro. 
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6.6 Valore estremo 


Si definisce valore estremo di un processo stocastico Gaussiano stazionario, X (t), la 
seguente quantità: 


Xe = Lx + 90x (6.59) 


dove g, noto come fattore di picco, è dato dalla seguente relazione: 


0.5772 
= /2ln(vT) + ————} 6.60 
g (vT) TT) (6.60) 


dove T è la lunghezza di X(t), e v è la frequenza di attraversamento di soglia: 


v= 


1 SÉ w2Sxx (w)dw ae! (Oe 


6.61 
27 ES Sxx(w)dw 27 oe ( ) 


Generazione di storie di vento 


In questo capito verranno proposte alcune delle più comuni tecniche per la gene- 
razione di storie di vento a partire dalla matrice di densità spettrale associata alla 
turbolenza atmosferica o, più in generale, alle forzanti acquisite in galleria del vento 
oppure mediante sperimentazione al vero. 

In questa sede, faremo riferimento alla componente fluttuante longitudinale della 
velocità del vento, u = u(t), intesa come un processo stocastico stazionario, gaus- 
siano, a media nulla. Chiaramente le tecniche di simulazione che verranno proposte, 
sono estendibili a qualsiasi processo stocastico con le stesse propietà statistiche. 


7.1 Simulazione di singoli processi stocastici 


La trattazione comincia con la simulazione di singoli processi stocastici autocorrelati. 


7.1.1 Constant Amplitude Wave Superposition (C.A.W.S.) 


Un processo stocastico stazionario Gaussiano a media nulla può essere inteso come la 
sovrapposizione di funzione armoniche. Pertanto, la velocità eolica u(t) potrà essere 
scritta nel seguente modo: 


N 


u(t) = HY” cos (2rnxt + ox) (7.1) 
k=1 
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dove H è un’ampiezza costante delle armoniche, N è il numero di bande di frequenza 
in cui si divide l’intero intervallo di frequenza!, ny sono le frequenze centrali delle 
varie bande di frequenza Anı ... Any, dx sono delle fasi casuali uniformemente 
distribuite tra 0 e 7. 

Il metodo si basa sul calcolo di H, in modo che lo spettro associato alla 7.1 sia 
quello desiderato. Dalla teoria dei processi stocastici sappiamo che la varianza, o°, 
di un processo stocastico stazionario, gaussiano, a media nulla, ed ergodico, è data 


dalla seguente relazione: 


+00 1 T 
oi =| Su(n)dn = F u?(t)dt (7.2) 
0 T 0 


dove T è la lunghezza della storia di vento, e S,,(n) è la densità spettrale (in questo 
caso unilatera) di u(t). Ricordando la 7.1, possiamo scrivere: 


2 
He cpl | N 
o? = ca f 3 cos(2rnxt +@x)| dt = Ho (7.3) 


Pertanto, l'ampiezza H sarà data dalla seguente relazione: 
H? = 03% (7.4) 


Per assicurare la convergenza del metodo, indicata con uy(t) la k — esima com- 
ponente di u(t), dovrà essere verificata la seguente relazione: 


H? 


1 T 
= Su(nx) Ang = = if uz(t)dt Yk (7.5) 
5 T Jo 


Tale condizione porta, infine, alla seguente relazione: 


H? oo 
— = += Gu (nk) Ang = cost. 7.6 
z = y Sua) Are (7.6) 
ossia, dato lo spettro dal quale si vuol generare una storia di vento, è necessario suddi- 


viderlo in N bande di frequenza tale che l’area sottesa dallo spettro in corrispondenza 
di queste bande sia costantemente pari a ci (Fig. 7.1). 

Il rapporto di convergenza di questo metodo è proporzionale a VIN. La genera- 
zione di N bande di frequenza, aventi frequenze centrali nı . . . ny può comportare un 
sostanziale onere computazionale. In ogni caso, alle alte frequenze si hanno intervalli 
troppi estesi (Fig. 7.1). 


! Generalmente la turbolenza atmosferica longitudinale ha una banda di frequenza compresa tra 0Hz 
e 4Hz. 
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Figura 7.1: Constant Amplitude Wave Superposition (C.A.W.S.) 


7.1.2 Weighted Amplitude Wave Superposition (W.A.W.S.) 


In questo metodo, le ampiezze delle armoniche che compongono il processo stoca- 
stico u(t), non sono costanti: 


N 
u(t) =) Hy cos (2rngt + dp) (7.7) 


k=1 


Per calcolare le ampiezze Hy si sfrutta la seguente uguaglianza: 


DE oa =1 f m Ordine ata 


da cui, infine, si ottiene che la k — esima ampiezza Hy, è data da: 
A; = 2S,(nx)Ank (7.9) 


Il rapporto di convergenza di questo metodo è proporzionale a N?. I vantaggi, 
rispetto al precedente metodo, risiedono nel fatto che il metodo può essere usato con- 
giuntamente con l’algoritmo FFT (Fast-Fourier-Transform), e le bande di frequenza 
hanno ampiezza costante (Fig. 7.2). 
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Figura 7.2: Weighted Amplitude Wave Superposition (W.A.W.S.) 


7.1.3 Filtri Autoregressivi 


Il processo stocastico, u(t), è scritto nel seguente modo: 


u(t) = X pu(t — kAt) + e(t) (7.10) 


k=1 


dove p è l’ordine del filtro, Yọ sono dei coefficienti da stimare, e e(t) è un rumore 
bianco a media nulla e varianza 02. In particolare, bisogna stimare i coefficienti yy. 


La funzione di auto-correlazione di un processo stocastico, ad esempio u(t), è data 
da: 


Rault) = E[u(t)u(t+7)] (7.11) 
che nel caso di un processo ergodico diviene: 
1 T 
R(t) = J u(t)u(t + 7)dr (7.12) 
0 


I coefficienti Yy si ottengono risolvendo il seguente sistema lineare: 


P 
Ra(0) = X` Yk Ru (kAt) + 02 


Kal (7.13) 
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cons=1...p. 

Il rapporto di convergenza del metodo é proporzionale al numero dei parametri da 
stimare. Tuttavia, non si conosce a priori l’ordine del filtro, p, che pertanto dovrà es- 
sere stimato iterativamente, ossia fino a quando non vi sia una buona corrispondenza 
tra lo spettro effettivo e quello simulato. 


7.2 Simulazione di processi multivariati 


Generalmente il campo di vento è un processo stocastico multidimensionale e multi- 
correlato nello spazio e nel tempo. Si procede, pertanto, all’estensione dei precedenti 
metodi al caso più generale e d’interesse per le strutture. 


7.2.1 Weighted Amplitude Wave Superposition (W.A.W.S.) 


Un vettore di velocita eoliche fluttuanti, definite nello spazio: 


u(t)= |: (7.14) 


puo essere simulato mediante la seguente relazione: 


N 
u(t) = 5 H(nx)cos(2rî, + Dx) (7.15) 
k=1 
dove: 
n, = nk + On, (7.16) 


sono i vettori delle frequenze centrali degli N intervalli di frequenza, con |n| << 
An/2; 


ð=] : (7.17) 
LIVORNO 


è il vettore delle fasi random coprese tra 0 e 7. 
Considerando, ad esempio, M = 2, si ha in forma esplicita: 


N Ai (ng) 0 0 cos (Qr7(1) kt + Pak) 
ug(t) = 5 Ho (nk) H22(nk) 0 COS (277 (2) pt + ®(2) x) 
ug(t)|  k=1(|H31(nx) Hs2(nk) H32(ng)] [cos (277i(3) gt + O(3),4) 
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(7.18) 
In maniera analoga a quanto visto per i singoli processi, si può dimostrare che: 
H(p) H (np)! = 2AnS, (nx) (7.19) 


ossia che la matrice H(ny) è il fattore di Cholesky di 2AnS,(nx) dove Su (ng) è 
la matrice di densità spettrale della turbolenza atmosferica, intesa ora come processo 
stocastico multivariato nello spazio e nel tempo. 

Per tenere in conto la correlazione del vento è possibile simulare l’i-esima com- 
ponente fluttuante del vento nel seguente modo: 


t) = DE D Hy; (nk) cos [ana (t + Ta) + D()) x (7.20) 
k=l j=l 
dove: 
Tij = AZ (7.21) 


OEO 

Pertanto, calcolato il profilo di vento per il sito in questione, si sceglie la dura- 
ta della storia di vento, T, il passo di campionamento, At, la frequenza massima 
dello spettro, il numero di intervalli di frequenza, N, e si calcola la matrice delle 
densità spettrali Sy (ng). Successivamente, si opera la fattorizzazione di Cholesky di 
2AnS,(nx), si generano i vettori fix e Py e si calcolano le storie di vento associate 
alle k-esime frequenze nz. A questo punto si itera per tutte le N bande di frequenza 
scelte, e di volta in volta si sommano le componenti armoniche ottenute. 


7.2.2 Filtri Autoregressivi 


Supponiamo di voler simulare M processi cross-correlati (come la velocità fluttuante 
del vento a diverse quote), allora scriveremo che l’i-esima componente fluttuante 
u;(t) è data da: 


SSL Juj(t — kAt)| + Ni(t) (7.22) 


k=1 j=i 


dove, al solito, p è l'ordine del filtro, e N(t) = [N;(t) ... Nag(t)] è il vettore compo- 
sto da M rumori cui corrisponde la matrice di cross-correlazione Rw) (7). 
Per determinare i coefficienti del filtro, Yj è necessario risolvere il seguente 
sistema di equazioni: 
p M 
Riml (lAt) > 5 Yizj(k ~~ I)At] 
a 1j=1 
p M 


ay valk) Rjm [kAt] + Rw) im(0 ) 


UE 


(7.23) 
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dovel=1...p. 

Per stimare N;(t) si opera la fattorizzazione di Cholesky della matrice di cross- 
correlazione di N;(t) e Nm(t), Rx N),im(0), ottenendo le matrice L detta fattore di 
Cholesky. Si puo dimostrare che 


N;(t) = 5 Lije;(t) (7.24) 
j=l 


dove €;(t) sono dei rumori bianchi a media nulla e varianza unitaria, completamente 
scorrelati tra loro. 

Pertanto, calcolato il profilo di vento per il sito in questione, si sceglie la durata 
della storia di vento, T, il passo di campionamento, At, e si calcolano i coefficienti 
del filtro 7,; dalla prima equazione del sistema 7.23. Successivamente, dalla seconda 
equazione del sistema 7.23 si calcola Rọ N);im(0), se ne opera la fattorizzazione di 
Cholesky, trovando L. Per ogni intervallo di tempo t;, si generano i processi e;(t;), 
si calcolano dunque i processi N(t;), ed infine il valore di u(t;). Il processo si ripete 
per ciascun istante t;. 


Dinamica delle strutture 


8.1 Sistemi ad un grado di libertà 


Le strutture possono essere assimilate a sistemi ad n gradi di libertà, riconducibili a 
n sistemi ad un solo grado di libertà. Pertanto, per affrontare lo studio delle strutture 
complesse, conviene introdurre la dinamica dell’oscillatore elementare. 


8.1.1 Oscillazioni libere 


Loscillatore semplice è caratterizzato da una massa m collegata ad una molla di 
rigidezza k che esercita su di essa una forza di richiamo elastico proporzionale allo 
spostamento della massa stessa (Fig. 8.1). 

Il moto delsistema è descritto dalla seguente equazione differenziale ordinaria 


x(t) 


m 


Figura 8.1: Oscillatore semplice 


omogenea di secondo ordine, lineare a coefficienti costanti: 


mà + kx =0 (8.1) 
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L'equazione 8.1 può essere riscritta nella forma: 
m+ wir =0 (8.2) 
previa introduzione della pulsazione naturale wo: 


k 

Wo = 4/ — (8.3) 
m 

La soluzione dell’equazione di moto (Eq.ni 8.1 e 8.2) fornisce la risposta dell’o- 


scillatore semplice, e può essere scritta come segue: 
x(t) = Csin(wot + ¢) (8.4) 


dove le costanti C e ¢ dipendono dalle condizioni iniziali. La soluzione dell’ oscilla- 
tore (Eq. 8.4) risalta il fatto che il moto dell’oscillatore e di tipo armonico caratteriz- 
zato da una pulsazione pari alla pulsazione naturale del sistema, ampiezza C, e fase 
Q. 


Date le condizioni iniziali: 


x(0) =x 
l ee (8.5) 
LO) = to 
si possono facilmente trovare le costanti d’integrazione: 
1.2 
C = | x2 + a 
wo (8.6) 
To 
tan(¢) = wo— 
To 


Nella pratica risultano essere interessanti i casi in cui £zo(0) = 0 oppure #,(0) = 


è la massa m inizialmente ferma e spostata dalla sua configurazione a riposo della 
quantità xo: 


i 7 o > x(t) = £o cos(wot) (8.7) 


e la massa m è spostata dalla sua configurazione a riposo da un impulso pari a 
Mio: 


so E i => x(t) = do sin(wot) (3.8) 
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Figura 8.2: Oscillazioni libere 


Le oscillazioni libere! sono oscillazioni armoniche, pertanto sono caratterizzate da 
un periodo proprio, To (Fig. 8.2): 


2 
T, = Z = m, 2 (8.9) 
Wo k 


ossia la risposta si ripete per multipli del periodo proprio di oscillazione. L’inverso 
del periodo proprio prende il nome di frequenza propria del sistema e sarà dunque 
data dalla seguente relazione: 


Le E (8.10) 
27 Vm 


8.1.2 Oscillazioni forzate 


Se l’oscillatore di Fig. 8.1 è soggetto ad una forzante armonica, F(t) = Fsin(wt), 
ľ equazione di moto assume la forma: 


F 
m + wr = —sin(wt) (8.11) 
m 


la cui soluzione é data dalla sovrapposizione della soluzione omogenea (Eq. 8.4) e di 
quella particolare. Quest’ultima si ottiene a partire dalla soluzione di tentativo 


xp(t) = Asin(wt) (8.12) 


ISi noti che l’oscillatore non è soggetto ad alcuna forza esterna, pertanto si parla di oscillazioni 
libere. 
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Sostituendo tale equazione nell’equazione di moto (Eq. 8.11) è facile trovare 
l’ampiezza A: 


F = F (8.13) 


mie) | f 3 S 


Pertanto è possibile accorpare le due soluzioni per ottenere la risposta di un oscilla- 
tore semplice forzato: 


x(t) = C sin(wot + ¢) 4 : z sin(wt) (8.14) 
1-(@ 


dove le costanti C' e ¢ si calcolano a partire dalle condizioni iniziali. 

In definitiva, la risposta è data dalla somma di due funzioni armoniche di cui una 
pulsa alla frequenza propria del sistema e l’altra alla frequenza della forzante F(t). 
Si può facilmente notare che se le due frequenze in questione sono uguali, la risposta 
dell’oscillatore risulta infinita. Tale situazione è nota come condizione di risonanza. 


8.1.3 Oscillazioni libere smorzate 


Nel caso in cui l’oscillatore è munito di uno smorzatore (Fig. 8.3), la massa sarà 
soggetta ad un’ulteriore forza di richiamo. Tale forza è di natura viscosa ed è propor- 
zionale alla velocità della massa stessa. L'equazione di moto sarà pertanto descritta 


Figura 8.3: Oscillatore smorzato 


dalla seguente equazione: 
m + ct + kz =0 (8.15) 


che, analogamente a quanto visto per l’oscillatore non smorzato, può essere riscritta 
come segue: 


ï + 2vwot + wea =0 (8.16) 
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con 
maree (8.17) 
m 


relazione che permette di definire il fattore di smorzamento, v: 


c c 
v= = (8.18) 
2MWo 2V km 
Le radici dell’ equazione caratteristica associata all’equazione di moto (Eq. 8.15) 
sono le seguenti: 


hiss (-v+ y2 1) Wo (8.19) 


Escludendo il caso in cui v > 1 per il quale la radici sono reali ed il moto è 
aperiodico, prendiamo in esame gli altri casi: 


e v = 1: le radici sono coincidenti, \j = Ag = —wo, e la soluzione risulta 
essere: 


a(t) = Are We! + Agte vot (8.20) 


In questo caso lo smorzamento del sistema corrisponde a quello critico pari a 
2Vkm 


e v < 1: Le radici sono complesse coniugate: 


A1,2 = —VWo F iwo V 1 — v? (8.21) 


e l'integrale dell’equazione di moto lo si può scrivere come: 


x(t) = Ce” cos (vov 1-r2t+ 6) (8.22) 


Il caso v < 1 è quello più frequente per le strutture, pertanto conviene studiarlo in 
dettaglio. Imponendo le condizioni iniziali: 


HO) =H (8.23) 
(0) =0 
si ottiene 


x(t) = ae = e Yet cos (wov 1-rv2t+ $) on 
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0.57 


Figura 8.4: Oscillazioni smorzate 


Se invece si pone: 


x(0) =0 
: : (8.25) 
“£(0) = do 
si ottiene 
x 
x(t) = —— el sin (woy 1 — vt 

ia ( ) (8.26) 

o=5 


Come si vede dalla Fig. 8.4, se non vi è carico applicato le oscillazioni decadono 
finché la massa non torna nella sua situazione iniziale di riposo. Il periodo di tempo 
che intercorre tra due massimi relativi successivi è pari a: 


27 
S woy 1 — v2 
Inoltre il rapporto tra due massimi relativi in Fig. 8.4 è pari a: 
2rv 
r= ex (8.28) 
È (- 1- 3) 
Il logaritmo di questo rapporto si definisce decremento logaritmico, 6: 
2 
ga (8.29) 


V1- p2 
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Se lo smorzamento è piccolo, v < 1, si ha: 
6 X 2Qrv (8.30) 


Infine se 6 è piccolo si ha: 


Tn 5 Tn — Tn+1 


e 1+ò (8.31) 


Tn+1 Tn+1 


8.1.4 Oscillazioni forzate con smorzamento 


x(t) 


m 


Figura 8.5: Oscillatore forzato con smorzamento 


Se applichiamo all’ oscillatore smorzato una forzante armonica (Fig. 8.5), avremo 
che il sistema sara descritto dalla seguente equazione differenziale ordinaria: 


F 
E+ 2vwot + wir = —sin(wt) (8.32) 
m 


All’integrale generale dell’omogenea (Eq. 8.26) andrà aggiunta la soluzione partico- 
lare che scriveremo: 


xp(t) = Asin(wt — y) (8.33) 
Sostituendo tale soluzione nell’equazione di moto si ottiene: 
1 F 
elt) = - n sin(wt — w) 
va — a2)? + 4y2a2 (8.34) 
2va 
tan(w) = ar) 


a Lodi 
dove si è posto a = &. 


Si può osservare dalla prima delle Eq.ni 8.36 che il fattore F'/k è il valore statico 
dell’ampiezza, ossia il valore che si avrebbe se il fenomeno fosse statico, mentre la 
quantità 


1 


(8.35) 
va — a2)? + 4y2o2 


94 LEZIONI DI INGEGNERIA DEL VENTO 


1 


(1-a 2)2+4 va 2 


Figura 8.6: Fattore di amplificazione 


assume il significato di fattore di amplificazione, mentre ) rappresenta lo sfasamento 
tra la risposta e la forzante. 

Per a = 1 si ha la risonanza ossia la massima risposta dell’oscillatore. In corri- 
spondenza della risonanza, se v = 0 l’amplificazione tende all’infinito, mentre per 
valori diversi la risposta tende a diminuire all’aumentare dello smorzamento (Fig. 
8.6). L'andamento della fase, +), è riportato in Fig. 8.7. 

In generale potremmo quindi scrivere la risposta totale del sistema come segue: 


x(t) = Ce? cos (wov 1—vt+ ) + 


1 F 
sin(wt — w) 
va — a2)? + 4y2a2 (39) 
2 
tan(w) = o> 


dove le costanti C' e ¢ dipendono dalle condizioni iniziali. 

Poiché l’oscillatore è smorzato, a parte un transitorio iniziale, le oscillazioni li- 
bere tendono a scomparire lasciando spazio alle sole oscillazioni forzate. Quando si 
perviene a tale situazione il sistema ha raggiunto la condizione di regime. 


8.1.5 Risposta ad una forzante qualsiasi 


Finora abbiamo trattato le forzanti armoniche, tuttavia la soluzione del problema di- 
namico può essere esteso al caso più generale di forzante che varia secondo una legge 
qualsiasi. 
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Figura 8.7: Andamento della fase 


A tal proposito è utile introdurre la risposta ad una forza F applicata all’oscil- 
latore all’istante to per la durata ôt. A tale forza corrisponde l’impulso elementare 
dI = Ft, mentre la risposta all’impulso elementare risulta essere: 


1 
h(t) = ——T == e sin (v 1- Pt) (8.37) 
( MWoV 1 — v2 s 
La risposta ad una forzante qualsiasi può essere vista come la somma delle ri- 
sposte corrispondente ai singoli impulsi elementari in cui si ipotizza di suddividere la 
forzante stessa. La sovrapposizione degli effetti è garantia dall’integrale di Duhamel: 


ole) — i F(7)h(t = r)dr (8.38) 


Pertanto la soluzione complessiva sarà: 


t 
x(t) = Ce” cos (wov 1— vt + ) + J F(r)h(t — 7)dr (8.39) 
0 


8.2 Sistemi a più gradi di libertà 


Le strutture reali sono sistemi aventi infiniti gradi di libertà, cosicché lo studio rigo- 
roso del loro comportamento dinamico richiederebbe un formulazione al continuo. 
Tuttavia, la teoria degli elementi finiti permette di discretizzare la struttura complessa 
in un sistema caratterizzato da un numero finito di gradi di libertà. La trattazione di 
tali sistemi è ulteriormente semplificata riconducendo il calcolo della risposta dina- 
mica alla sovrapposizione della risposta di singoli oscillatori elementari smorzati e 
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forzati. In definitiva, un sistema lineare a n gradi di libertà può essere scomposto in 
n sistemi ad un solo grado di libertà, per i quali vale la trattazione svolta finora. La 
suddetta scomposizione è nota come decomposizione modale e si basa sul calcolo dei 
modi propri di vibrare. 


8.2.1 Oscillazioni libere non smorzate 


Le oscillazioni libere di una struttura non smorzata, ricondotta ad un sistema a n gradi 
di libertà, si determinano risolvendo un sistema di equazioni differenziali ordinarie 
che potremo scrivere in forma matriciale come segue: 


MU +KU =0 (8.40) 


dove M e K sono rispettivamente la matrice delle masse e la matrice di rigidezza del 
sistema: 


Mii Mia ... Min ku k12 e kin 
M21 Mar ... Man ko ka ... kon 

M=] . E, . K=]| . i ; (8.41) 
Mni Mn2 +--+ Mnn kni kn2 ... knn 


mentre Ü e U sono rispettivamente i vettori delle accelerazioni e degli spostamenti 
della struttura: 


Ul Ul 

A U2 U2 

U= . U = f (8.42) 
Un Un 


Se la struttura non é labile si dimostra che M e K sono definite positive, per- 
tanto la soluzione del problema la si può scrivere, applicando il principio di sovrap- 
posizione degli effetti, come la somma delle risposte di n sistemi ad un grado di 
libertà: 


n 


U(t) = PQH) = $ [C cos(wjt + 4;)] 4y (8.43) 
j=l 
le costanti C} e j dipendono dalle condizioni iniziali, mentre ® = [d1, d2,..., dn] è 


la matrice delle forme modali o modi naturali di vibrare del sistema, ed è una matrice 
reale. Il vettore Q(t) = [q1.92,..., In” è il vettore delle singole risposte associate al 
j-esimo oscillatore elementare: analogamente a quanto visto precedentemente (§ 8.1) 
avremo: 


qj = C; cos(wjt + W;) (8.44) 
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Per la stessa analogia andremo a definire la j-esima pulsazione propria, w;, il 
j-esimo periodo proprio, T}, e la j-esima frequenza propria, fj, della struttura: 
o 27 1 


haz (8.45) 


wj =2rf; T; 
J J J wj T; 


Le quantità wj e @; si determinano trovando le soluzioni non banali del sistema di 
equazioni differenziali che governa il moto (Eq. 8.40), proiettato nello spazio del- 
le coordinate modali. In particolare il problema iniziale si riconduce al seguente: 
ricordando le relazioni 8.40 e 8.43 si ha: 


M®Q+K9Q=0 (8.46) 


Considerando una soluzione del tipo q; = C; cos(wjt + w;) è facile dimostrare che 
il problema diviene: 


—wiM¢,;q; + K6;g; = 0 (8.47) 
che può essere riscritta in forma compatta: 
[K — w? M] $j4;= 0 (8.48) 


Escludendo le soluzioni banali del tipo q; = 0, vogliamo trovare le quantità w; e 
Qj tali per cui: 


[K — ©M] ¢; =0 (8.49) 


Il problema si risolve mediante la teoria spettrale. In particolare, siamo alla ri- 
cerca delle radici del polinomio caratteristico associato al sistema. Tale polinomio è 
dato dalla seguente relazione: 


Det [K — \M]=0 (8.50) 


dove le radici Aj assumeno il significato di autovalori 
Otteniamo dunque n valori di A tali che 0 < Ay < Ag < ... < An. Le pulsazioni 
proprie del sistema saranno date dalla seguente uguaglianza: 


we = Àj (8.51) 


Le forme modali o autovettori, dj, si trovano applicando le equazioni definite 
nella 8.49. Fissato j si ha un sistema di equazioni lineare omogeneo, la cui soluzioni 
Qj, è definita a meno di una costante di proporzionalità. Per eliminare tale indetermi- 
nazione, basta porre una delle sue componenti pari all’unità, ad esempio la prima. 

Si può dimostrare che gli autovettori sono vettori ortogonali rispetto alla matrice 
delle masse e alla matrice di rigidezza della struttura, ossia sono tali da diagonaliz- 
zare tali matrici. Poiché ciascun autovettore viene determinato a meno di costante 
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arbitraria, generalmente si opera una normalizzazione in modo tale che la matrice 
degli autovettori, ®, sia anche ortonormale rispetto alla matrice delle masse: 


TM? =I (8.52) 
dove J è la matrice d’identita. Inoltre, in questo caso si ha che: 
&’Kd = (8.53) 
avendo posto: 
wt 
Q? = É 9 (8.54) 


2 
Wn 


Le propietà di ortogonalità e ortonormalità della matrice degli autoettori, ®, per- 
mettono di semplificare il calcolo della risposta dinamica di una struttura complessa. 
Infatti, premoltiplicando tutti i termini della 8.46 per T si ottiene: 


PTMI + PTKƏQ =0 (8.55) 
da cui per quanto detto si ottiene: 
IQ+9°2Q=0 (8.56) 


che corrisponde ad n equazioni differenziali lineari omogenee linearmente indipen- 
denti: 


di twin =0 


do + 0392 = 0 
(8.57) 


Gn = Wi Gn =0 


ossia la struttura viene scomposta in n oscillatori elementari le cui risposte sono date, 
come visto in precedenza, dalle seguenti soluzioni: 


qı = Ci cos(wit + 41) 
q2 = Ca cos(wst + Y2) 


(8.58) 


qn = Cn cos(wnt + Un) 
con le costanti C; e =; da determinare in base alle condizioni iniziali. La soluzione 
infine sarà data dalla sovrapposizione degli effetti: 
n 
U(t) = Y_[C; cos(w;t +%;)] bj (8.59) 


j=l 
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8.3 Oscillazioni libere smorzate 


Prima di entrare nella trattazione analitica delle oscillazioni libere smorzate conviene 
sottolineare che vi è un grande numero di incertezze nella valutazione dello smorza- 
mento delle strutture reali. Le cause principali di queste incertezze sono di varia natu- 
ra: l’attrito dei giunti, l’isteresi del materiale, la presenza di elementi non strutturali 
interagenti con la struttura, le resistenze del mezzo. Generalmente non è possibile 
tenere in conto tutte queste cause, peraltro non sempre facili da individuare, ma si 
preferisce valutarne l’effetto complessivo mediante indicazioni di natura sperimen- 
tale. In quest’ ottica risultano giustificate alcune semplificazioni che possono essere 
sintetizzate in due ipotesi fondamentali: 


e Le forze di richiamo viscoso, resistenze viscose, sono proporzionali alla velo- 
cità; 


e La matrice di smorzamento, C, è esprimibile come una combinazione lineare 
della matrice delle masse e delle rigidezze. 


Quest’ ultima ipotesi può essere riscritta con l’ausilio di due costanti, a e 8: 
C=aM+ 6K (8.60) 


La 8.60 costituisce una condizione sufficiente affinché la matrice di dissipazione sia 
diagonalizzabile. 

L'importanza della diagonalizzazione C è evidente. La dinamica di una struttura 
complessa smorzata è descritta in forma matriciale nel seguente modo: 


MU +CU +KU =0 (8.61) 


introducendo nella 8.61 la traformazione U(t) = ®Q(t) e premoltiplicando tutti i 
termini per T si ottiene: 


1Q+NQ+’Q=0 (8.62) 
dove 
2V1WI 0 0 
T 0 2V2w9 wee 0 
N = Cé= . . A f (8.63) 
0 0 0° 2UnWn 


è la matrice di dissipazione disaccoppiata. 
Le equazioni che descrivono il moto (Eq. 8.62) risultano dunque essere linear- 
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mente indipendenti: 


di + 2ywigi + wig =0 


do + 220292 + waq = 0 
(8.64) 


dn + 2UnWnln + Wi In =0 


ossia la struttura viene scomposta in n oscillatori smorzati le cui risposte sono date, 
come visto in precedenza, dalle seguenti soluzioni: 


a(t) = Cie 11° cos (on /l—vet+ v) 
th=C —vzw2t 1 — v2t + 
q2 (t) ve cos (cav V5 vs) (8.65) 


dn(t) = Ce E cos (@nvi — VIt+ Va) 


con le costanti C} e 4; da determinare in base alle condizioni iniziali. 
La soluzione sarà data dal principio di sovrapposizione degli effetti: 


U(t) = 5 [Cje 54" cos Gr fjl—- vit + vi)| dj (3.66) 


j=1 


8.3.1 Oscillazioni forzate con smorzamento 


Il problema delle oscillazioni forzate è descritto dalle seguenti equazioni date in 
forma matriciale: 


MU + CÙ + KU = F(t) (8.67) 
operando la decomposizione modale si ha: 

IÖ + NQ +Q = P(t) (8.68) 
dove 

P(t) = PT F(t) (8.69) 


è la proiezione del carico sulle forme modali. Nella precedente formulazione si 
intende che 


F(t) = [Fi (0), Fa(t),...,Fa()] PE = [P(t), Po(t),..., Pat] (8.70 
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La struttura viene quindi scomposta in n oscillatori smorzati e forzati: 


di + 2ywigi + win = p(t) 


do + 2v.woge + w392 = pa(t) 


(8.71) 


Gn + 2UnWndn + wen = Pn(t) 


La soluzione dipende dal tipo di forzante. 


8.3.2 Risposta a forzanti qualsiasi 


Si procede al calcolo della risposta per una forzante multidimensionale F(t) di tipo 
qualsiasi, essendo tale soluzione di validità generale. 

Le equazioni scritte nella 8.61 sono linearmente indipendenti, in quanto disac- 
coppiate. La risposta del j-esimo oscillatore sarà data dalla sovrapposizione della 
j-esima soluzione omogenea e della j-esima soluzione particolare. La soluzione omo- 
genea è data dalle equazioni 8.65, mentre, analogamente a quanto visto nel paragrafo 
8.1.5, la soluzione particolare sarà espressa mediante l’integrale di convoluzione, o 
di Duhamel. 


t 
q(t) = Cie 11" cos (è 1—vet+y ) f pı(T)hı(t— T)dT 
i(t) = Ci 1 itty F i(7)hi(t — 7) 
t 
qo(t) = Coe?" cos (è ber ) f po(T)ho(t — 7)dr 
2(t) = C2 2V gt + Y2 : 2(T)ha(t — 7) 


(8.72) 
t 
dn(t) = Ce?" cos (env 1—v2t+ Un) +f Pn (T)hn(t — T)dr 
0 
dove: 
hj(t) = a Sin e Vi%it sin (o; 1- vt) (8.73) 


wj 1—13 


è la riposta ad un impulso unitario dell j-esimo oscillatore semplice. Infine, la solu- 
zione sarà data dalla sovrapposizione degli effetti: 


n t 
U(t) = 5 jer cos (o; /1— vt + vi) +f p(T)hylt—7)dr| Qj 


j=l 
(8.74) 


ricordando che nella condizione di regime l’effetto delle soluzioni omogenee svanisce 
in quanto la struttura è smorzata. 


Dinamica aleatoria delle strutture 


Nell’ambito della dinamica deterministica di un oscillatore ad un grado di liberta, si 
assume che ciascun parametro definito nell’equazione di moto é deterministico. In 
tali condizioni, la risposta è unicamente determinata una volta imposte le condizioni 
iniziali. Se almeno uno di questi parametri è aleatorio, la risposta sarà a sua volta 
aleatoria. Ciò significa che la risposta dovrà essere necessariamente caratterizzata 
statisticamente. 

In questa trattazione si assume che il carico sia un processo stocastico, dal momento 
che generalmente il vento può essere inteso come un processo stocastico multidimen- 
sionale e multivariato. 


9.1 Risposta aleatoria dell’oscillatore smorzato 


Nel Cap. 8 abbiamo introdotto l’equazione di moto dell’oscillatore smorzato e forza- 
to: 


1 

# + Qvwot + wre = — f(t) (9.1) 
m 

In questo caso f(t) è un processo stocastico, applicato all’oscillatore nell’istante 


t=0. 
La soluzione del problema (cfr. Eq. 8.39) è data dalla seguente relazione: 


a(t) = Ce” cos (wov 1-v2t+ ) + [ f(7)h(t — 7)dr (9.2) 
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dove il primo termine, soluzione omogenea, dipende solo dalle condizioni iniziali e 
quindi é deterministico. Poiché vogliamo caratterizzare la risposta aleatoria, possia- 
mo tralasciare la parte deterministica. Prendiamo cosi in considerazione la seguente 


soluzione: 
W= fs f(r)h(t — 7)da (9.3) 


dove la risposta ad un impulso, h(t), è data da: 


1 
h(t) = “ep = e “ot sinwoV 1 — vt (9.4) 


Dalla teoria dei processi stocastici sappiamo che per poter caratterizzare la rispo- 
sta aleatoria dobbiamo calcolare i suoi momenti statistici. L’n-esimo momento è dato 
da: 


E [x( (t2)...2(tn)] = 


ltd ie h(t, — 71) e f(Tn)hltn — m)dTi dn = 


-f 1 di E[f(71) <<< f(tm)] halti — T) h(tn— Tr)dt1 +++ din 
(9.5) 


ottenuto ricordando la linearità del operatore valore atteso, E|. ..]. Analogamente si 
può scrivere l’espressione dell’n-esimo cumulante! 


Kn [e(t1)x(t2) (tn) = 


ti tn 
= J De | kn [f(71) +++ f(t) | h(a 1) htn — T)d71 +++ din 
0 0 
(9.6) 
Un caso particolare è quello di forzante Gaussiana che produce una risposta Gaus- 


siana. In questo caso i momenti di ordine superiore al primo sono nulli, mentre i 
primi due assumono rispettivamente il significato di funzione di media e funzione di 


'Si ricorda che il primo cumulante è uguale al primo momento, mentre il secondo e terzo cumulante 
sono coincidenti con i corrispondenti momenti centrali, ossia: 


1 [x(t)] = Ele] = pelt) 
rale(tı)z(t2)] = B{le(ts) — pe (t1)lfe(t2) — pe (t2)]} = Rea (t1, ta) — p(t) pte (t2) 
rale(t)e(ts)e(ts)] = E{[e(t) — wo(tr)|[e(te) — wo(ta)|[e(ts) — wo(ts)]} 


mentre i cumulanti di ordine superiore non sono facilmente calcolabili. 


[x 
[x 
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covarianza: 


E [x(t)] = pelt) = I E[f(r)] h(t — 7)dr 
Kea (t1, t2) = Reo(ti, t2) — Ue (ti) Ue(te) = (9.7) 


ty t2 
= f Kp (T1, T2)h(ti — Ti)h(t2 — 72)dr1dr2 
0 0 


dove Rx è la funzione di autocorrelazione della risposta x(t). 


9.2 Forzanti aleatorie debolmente stazionarie 


In diversi problemi di ordine pratico, la forzante può essere assunta debolmente 
stazionaria, condizione verificata quando: 


EFO] = ep (t) = ep Vt (9.8) 
E[f?(t)] = cost. Vt (9.9) 

Si può dimostrare che quando il processo è debolmente stazionario la funzio- 

ne di autocorrelazione, R;;, non dipende dai singoli instanti tı e to ma dalla loro 


differenza, tı — t2 = T. Pertanto: 


Elf) f(t2)] = Rep — to) = Rpy(7) (9.10) 


con T = tı — to. 
La funzione di autocorrelazione può essere calcolata anche operando la trasfor- 
mata inversa di Fourier della densità spettrale di potenza della forzante, S;y(w): 


+00 i 
Rya ta) = f Spy(w)el =) du (9.11) 


Utilizzando questa relazione è possibile trovare la funzione di autocorrelazione 
della risposta: 


Rra(t1,t2) = Ela(t1); x(t2)] = 


ti t2 

| {EIA I - ecm = 
ti n 

=| f Re p(t h(t1 — 71) A(t — T2)dt1dT2 = 


ty tg +00 . 
= ni J r Spy) MT h(t = T1)h(ta = T2) dr; dradw (9.12) 
0 0 - 0 
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Ponendo u = ty — T4 € v = t2 — 73 ed integrando per ue v si ottiene: 


+00 . 
Ree(t1, t2) = i Sip p(w) M(w, t1) MN" (w, tee) du (9.13) 
con: 
t è 
M(w, t) =} h(u)e “du (9.14) 
0 


La funzione h(.) è la risposta ad un impulso unitario, mentre M(w, t) è una funzione 
complessa, pertanto M* (w, t) è il suo complesso coniugato. Si può osservare che la 
risposta ad una forzante stazionaria sembrerebbe non essere stazionaria perché questa 
è funzione di tı e t2 piuttosto che di tı — t2. Tuttavia tenendo presente che la funzione 
di risposta ad un impulso unitario è nulla per t < 0 è possibile cambiare l’estremo 
d’integrazione inferiore della 9.14 e porlo pari a —oo: 


t 

M(w,t) = f h(uje™™"du (9.15) 
=00 

se poi consideriamo la condizione di regime possiamo porre t + +00, e otteniamo 

che la 9.15 diviene la trasformata di Fourier della funzione di risposta ad un impulso 

unitario, nota come ammettenza strutturale H(w): 


1 t —iwu 1 1 
H(w) = — h(u)e “du = — (9.16) 
"1 "ap gE] 
wi Wo 
Possiamo dunque scrivere la funzione di autocorrelazione della risposta stazionaria: 


+00 . 
Rie(t1 — t2) = f |H(w)|? Sjp(w)e 12) do (9.17) 


—co 


Si può notare che la funzione a secondo membro della 9.17 è la trasformata inver- 
sa di Fourier della quantita |H(w)|?.S;y(w). Ciò permette dunque di trovare la densità 
spettrale della risposta stazionaria che per definizione è la trasformata di Fourier della 
funzione di autocorrelazione Rx (t1 — t2): 


ii = | S Ras(ti — to)e 2 d(ty — ty) (9.18) 
da cui risulta infine: 

Sre(w) = |H(w)|? Sp;(w) (9.19) 
con: 

HWw)P = : (9.20) 


24 2 2 
mg) +02 (2) 


Risolto il problema dinamico, la valutazione della risposta dell’ oscillatore è carat- 
terizzata dai primi due momenti statistici, essendo la forzante debolmente stazionaria 
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9.2.1 Valutazione della risposta media 


In base a quanto visto in precedenza, la risposta dell’oscilatore in condizione di 
regime è data da: 


x(t) = ie f(r) h(t — r)dr (9.21) 


Pertanto il valore medio della risposta, juz, sara dato da: 


+00 
ua = Ele) = E | JOME- rar] = 
+00 
= DO) f h(t — r)dr = pyH(0) (9.22) 


Avendo sfruttando la linearità dell’ operatore valore medio, E |...], e ricordando 
che H(0) = [72° h(t — 7)dr. Infine, notando che 


1 1 


si perviene alla ben nota relazione, di valenza generale, tra spostamenti e forze medie: 


Gi 
eet 9.24 
Ha = 7 (9.24) 


9.2.2 Valutazione dei momenti del secondo ordine 


Operando la trasformata inversa della densità spettrale della risposta si ottiene la 
funzione di correlazione della risposta debolmente stazionaria: 
+00 ; 
Ryo (7) = Ela(t1)x(t2)] = f Sra (we? dr (9.25) 


—oo 


avendo posto tı — tg = T. 
A partire da Ryg è possibile trovare il momento centrale del secondo ordine, ossia 
la funzione di covarianza, kx, associata alla risposta debolmente stazionaria: 


Ko(T) = E {[x(t1) — pel[e(t2) — Mal} = Res(T) — pz (9.26) 
ponendo tı — to = T. Infine, ponendo 7 = 0 si ottiene la varianza della risposta: 
kzz(0) = o2 (9.27) 


La risposta aleatoria dell’ oscilaltore risulta così completamente caratterizzata. 
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9.3 Risposta dei sistemi a più gradi di libertà 


La valutazione della risposta aleatoria di una struttura complessa può essere ricon- 
dotta allo studio di oscillatori ad un solo grado di libertà mediante la decomposizione 
modale (Cap. 8). 

Consideriamo il problema dinamico di un sistema smorzato, a più gradi di libertà, 
soggetto ad una forzante aleatoria multidimensionale ($ 8.2): 


MU + CU + KU = F(t) (9.28) 


Risolto il problema agli autovalori, w;, e agli autovettori, $;, è possibile proiettare 
il problema nello spazio delle coordinate modali, qj, scrivendo la soluzione come 
combinazione lineare delle coordinate modali: 


U(t) = &Q(t) (9.29) 


Si ricorda che in generale gli autovettori costituiscono una base ortogonale, ossia 
sono tali da poter diagonalizzare la matrice delle masse e delle rigidezza. Inoltre, 
operando un’opportuna normalizzazione è possibile rendere gli autovettori anche or- 
tonormali rispetto alla matrice delle masse, che risulta così essere pari alla matrice 
d’identità. 

In questa trattazione, si preferisce considerare una base ortogonale, ma non orto- 
normale, di autovettori. Il significato sarà chiaro nel seguito. 

Sostituendo la 9.29 nella 9.28 e moltiplicando entrambi i membri per ®7 si 
ottiene la decomposizione modale del problema: 


MQ +C + KQ = P(t) (9.30) 
dove: 
P(t) = ®T F(t) (9.31) 
e 
o 0 0 & 0 0 
0 0 0 0 
M= ji du (9.32) 
0 0 On 0 0 0 & 
K1 0 0 
0 ar MEE T 0 
K=]. . , f (9.33) 
0 0 0 kr, 


Il problema si riconduce pertanto al calcolo di un sistema di equazioni differenziali 
linearmente indipendenti. 
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9.4 Forzanti aleatorie debolmente stazionarie 


Abbiamo più volte ricordato che per caratterizzare un processo stocastico debolmente 
stazionario è sufficiente calcolare i primi due momenti statistici. Tuttavia, è impor- 
tante ricordare che la risposta della struttura sarà stazionaria solo a regime (cfr. $ 9.2), 
condizione che si avvererà quando l’effetto della soluzione omogenea del problema 
svanisce. 

Consideriamo il sistema di equazioni che risolve la struttura: 


pui + 191 + K191 = pi(t) 
pago + odo + K2q1 = p2(t) 


(9.34) 
undn + É10n + kng = Prlt) 
che può essere riscritto come: 
x : 2 1 
Gt 2017191 + wiqi = aa 
i, : 2 1 
d2 + 2vqweqz + wiqo = —po(t) 
02 (9.35) 
3 ; 2 1 
Qn + 2UnWnAn + WnIn = — Pn (t) 


trascurando la soluzione omogenea, avremo che la j-esima risposta sarà data dalla 
sola soluzione particolare: 


q(t) = si hj(t— T)pj(7)dr (9.36) 


dove al solito h(t) è la funzione di risposta ad un impulso unitario: 


hj(t) = —_ eit sin (w;/1- 84) (9.37) 


0jW54/1— Vi 


Ricordando che h(t) è per definizione nulla per t < 0 possiamo cambiare l’estre- 
mo inferiore di integrazione della 9.36: 


q;(t) = f hj(t — T)p;(T)dT (9.38) 


_- 00 
Nell'ipotesi in cui la struttura ha raggiunto la situazione di regime, possiamo 
cambiare anche l’estremo superiore di integrazione della 9.36, t + +00: 


q(t) = f n hj(t — 7)p(7)dr (9.39) 


—co 
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Poiché le forzanti sono stazionarie e la struttura ha raggiunto la condizione di 
regime, la risposta sara a sua volta stazionaria. In particolare, essendo le forzanti 
debolmente stazionarie, per caratterizzare il processo stocastico multidimensionale e 
multivariato della risposta, sara necessario calcolare i suoi primi due momenti sta- 
tistici. Procediamo quindi con il calcolo della funzione di cross-correlazione della 
risposta totale, ricordando che w(t) = Y5_, $14; (t) avremo: 


Ruyum (t1 — t2) = E [wi(t1)um(ta)] = 


n n +00 +00 
PIO a 1; PkmE [pj (71) Pe (72)] X 


Pr Ee 

x hj(t1— 71) he (te — T2)dr1dr2 = 
Nn n +00 +00 

= >S Ji DrjPrmBo;p (Tı — T2)hj(t1— 71) he (ta — T2)d rd 
j=1 k=1“ 700 “790 


(9.40) 


La funzione di cross-correlazione si ottiene applicando la trasformata inversa di 
Fourier della densità spettrale della risposta: 


+00 . 
Rosy (Tı = 72) = f Soo EA dy (9.41) 
pertanto si ottiene: 
Tag se +oo +00 +00 
Ruyum(t1—t2) = f / | 1; Pkmhj (ta — T1) he (te — T2) x 
JELRE a a eo 


> Spip (we —™ dry drgdw (9.42) 


Ponendo u = ty — 7] € v = to — 73 ed integrando per u e v avremo: 


n n +00 . 
Ruun (ti —t2) = DR: $1; Pim Hi; (W) Spip, (0) H} (we) dw (9.43) 


j=l k=1" 7 


che si ottiene ricordando che l’ammettenza strutturale è la trasformata di Fourier della 
funzione di risposta ad un impulso unitario: 


+00 
Hj(w) -f h;(u)e “du (9.44) 


io la 


Si può facilmente notare dalla relazione 9.43 che il secondo membro è, a me- 
no delle forme modali, la trasformata inversa di Fourier di Hj(w)Sp,p,(w) H% (w). 
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Si può dunque calcolare la densità spettrale della risposta che per definizione è la 
trasformata di Fourier di R.,,,,,(t1 — ta): 


Surum (4) = DD by bem HH; lw) Spp, (7) HE (1) (9.45) 
j=1 k=1 


Complessivamente si ottiene la matrice delle densità spettrali della risposta strut- 
turale, Syy (w): 


Suu(w) = PH (w)Spp(w)H* (wo! (9.46) 
dove H(w) e Spp(w) sono rispettivamente la matrice delle risposte impulsive e la 


matrice delle densità spettrali dei carichi modali. In generale le matrici che compaio- 
no nella 9.46 sono definite come segue: 


Su,u,(w) Suu:(W) ... Suu, (w) 
Su. S cot ASU 
SRA UaU (w) vat, (w) vaun (w) cay 
SU„U (w) SU„U2 (w) sons SUnUn (w) 
Sp.p(w) Sp,p.(W) ... Sp.p,(w) 
S 211 S 212 STE S 21n 
Spp(w) = PoP (w) PaPa(w) PoP (w) (048) 
Sp, Pi (w) SP, P> (w) Fei SP, Pa (w) 
b. i 1 0 0 
SCOLO) 
H(w) Di 0 o da 0 
0 (E e ee 
en (3) +i (4 )| 
(9.49) 


Risolto il problema dinamico la valutazione della risposta strutturale si completa 
trovando i primi due momenti statistici della risposta”. 


*Le forzanti applicate sono debolmente stazionarie, pertanto soltanto i primi due momenti statistici 
sono stazionari. 
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9.4.1 Valutazione della risposta media 


In base a quanto visto in precedenza, proiettando gli spostamenti strutturali nello 
spazio delle coordinate modali, otteniamo: 


U(t) = DQ(1) (9.50) 


Pertanto il valore medio della j-esima risposta, uy sara dato da: 


uu = E[U(t)] = E[Q(t)| = DE | L h(t — 7)P(r)dr (9.51) 
dove: 
hat) 0... 0 
h(t) = . a ta (9.52) 
È oe 


Sfruttando la linearita della funzione valore medio, possiamo riscrivere la 9.51 
come segue: 


py = ®E if. h(t — PPU) =6 ie h(t — 7) E[P(t)]dr = 


—co —oo 


=o i h(t —7)E [ÐT F(t)| dr = up® / 


—co —co 


+00 
h(t — nar] dI = 


= up®PH(0)DT (9.53) 


ottenuta ricordando che P(t) = DT F(t) e H(0) = [72° h(t — r)dr. 
Infine, si può dimostrare che 


6H(0)6? = K7! (9.54) 


pertanto si perviene alla ben nota relazione, di valenza generale, tra spostamenti e 
forze medie: 


ig Sie (9.55) 


9.4.2 Valutazione dei momenti del secondo ordine 


Operando la trasformata inversa della matrice delle densita spettrali della risposta 
si ottiene la matrice delle funzioni di cross-correlazione della risposta debolmente 
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stazionaria: 


+00 . 
Ruu(7)= E[U(t)U (ta)] = f Suu(w)eTdr = 


Ruu,(7) Ruvz(7) --- Ruun(7) 
Ry. Ruy. ioe GR. 

_ ni (7) va» (7) f 03, (7) 050 
Ru,U, (7) Ru,us(7) eS Ru,U, (7) 


avendo posto tı — tg = T. 

A partire da Ruy è possibile trovare il momento centrale del secondo ordi- 
ne, ossia la matrice delle funzioni cross-covarianza, kuyuy, associata alla risposta 
debolmente stazionaria: 


kuult) = E{[U(t1) — wo] [U (t2) — uul} = Ruu(7) — uug (9.57) 


dove si è posto tı — ta = T. Infine, ponendo 7 = 0 si ottiene la matrice di covarianza 
della risposta: 


2 
TUU YUz see OU Un 
OU2U1 OUUU + ~FU2Un 
xuu(0) = i f i (9.58) 
2 
OU,U, TUnUg «+: TU,Un 


che nel caso di processo stazionario coincide con Ryy(0), essendo la media stazio- 
naria e quindi potendola considerare a parte. Tuttavia, generalmente, è preferibile 
considerare solo le componenti fluttuanti delle forzanti, in modo da avere dei proces- 
si stocastici a media nulla. In tal caso l’identità xyy(0) = Ryu (0) è rigorosa. 

A questo punto la risposta aleatoria della struttura è completamente caratterizzata. 


9.5 Formulazione generale della dinamica dei continui 


L'equazione del moto di un continuo deformabile è data dalla seguente relazione 
definita nel generico punto individuato dal vettore r = x -i1+y-j+z-k: 


m(r)U(r, t) + e(r)U(r, t) + LIU(1)] = F(x, t) (9.59) 


con il consueto significato dei simboli fatta eccezione per il funzionale £ che dipende 
dal tipo di struttura in esame. Ad esempio, per una piastra si ha 


Eh3 04U(x,y,1) Ulz yt) OU (x,y, t) 


are di dog > 


(9.60) 
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dove h è lo spessore della piastra E è il modulo di Young, e vp è il coefficiente di 
Poisson. Per una struttura snella, che si sviluppa prevalentemente in direzione z, si 
ha invece: 


o 


aa (9.61) 


LU (r, t] = 


È Rea 


Oz? 


In ogni caso, con un’opportuna scelta di coordinate è possibile decomporre il 
problema analogamente a quanto visto per i sistemi discreti. Proiettando il problema 
nello spazio delle coordinate modali, @;(r) con j = 1,2, .. . , +00, si ottiene: 


+00 
=) 4y) (9.62) 


Le coordinate modali sono tali per cui lungo l’intera struttura, S, si ha: 


fot gilr (r)dS=0 if#j (9.63) 
ee gilr (r)dS=0 ij (9.64) 
x gilr \dS=0 i#j (9.65) 

ed inoltre: 
Li[di(r)] = m(x)w7di(r) (9.66) 


dove w; = 27n, è l’i-esima frequenza angolare della struttura. 
Pertanto, sostituendo la 9.62 nella 9.59, pre-moltiplicando l’equazione di moto 
per $;(r) ed integrando su tutta la struttura S, si ottiene: 


Ís ®i(r) F(r, t)dS _ mlt) 
Js Blaym(a)dS — mi 


= ni di (1) F(r, t)dS (9.68) 
Ss 


è l’i-esimo carico proiettato sulla forma modale ¢;(r), e 


Qi + 2Viwidi + wi di = i=1,2,...,+00 (9.67) 


dove 


mi = f p? (r)m(r)dS (9.69) 
S 


è la massa associata all’ i-esimo modo di vibrare. 
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9.5.1 Risposta ad una forzante Gaussiana stazionaria 


Come più volte menzionato, la risposta totale indotta da una forzante Gaussiana sta- 
zionaria è data dalla componente media e da quella fluttuante. Inoltre, nell’ipotesi di 
Gaussianità, i primi due momenti statistici caratterizzano completamente la risposta 
strutturale. 


Risposta media 


In questo caso la 9.59 si riduce a: 
LU (r)|] = ur(r) (9.70) 


ossia la risposta media, uy (r), si riconduce all’integrazione dell’equazione differen- 
ziale della linea elastica della struttura soggetta al carico medio ur(r). 


Risposta dinamica 


La risposta dinamica la si può trovare nel dominio del tempo integrando le equazioni 

di moto (Eq. 9.67) e poi sommarle in base alla 9.62, oppure nel dominio delle fre- 

quenze. In ogni caso bisogna caratterizzare il momento statistico del secondo ordine. 
Abbiamo visto che l’i-esima risposta associata alla struttura è data da: 


TE f "aci 071) 


dove al solito h;(t) e la funzione di risposta ad un impulso unitario: 


1 
hit) = —— e sin (w;/1- vi) 9.72 
( (siii (0.72) 


midi Jl- vi 


Ricordando che h(t) è per definizione nulla per t < 0 possiamo cambiare l’estre- 
mo inferiore di integrazione della 9.71: 


ale J de (0.73) 


Nell’ ipotesi in cui la struttura abbia raggiunto la condizione di regime possiamo 
cambiare anche l’estremo superiore di integrazione della 9.71, t + +00: 


att) = | Aes 0.74) 


—co 
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Analogamente a quanto visto per i sistemi discreti si può calcolare la funzione di 
cross-correlazione: 


Ryu (r1,4e, ti — t2) = E [U (r1, t1)U (r2, t2)] = 


n n 


TOE 


Jaiki YS 
x halts — Ti)hg(t2 — T2)dTidrə = 


n n 


+oo too 
=. a bj (1) be (2) Ro, pa (71 — T9)X 


PI 
x hj(t1 — T1)he(te2 — T2)drıdr2 (9.75) 


La funzione di cross-correlazione si ottiene applicando la trasformata inversa di 
Fourier della densità spettrale della risposta: 


+00 . 
Rpip: (Tı — T2) = / Spin (w) dw (9.76) 


—00 


pertanto si ottiene: 


Ruu(r1,r2,t1- te) = 


n n 


E DI ~ Poni dj(11)dx(r2)hj(t1 — T1) he (ta — T2) x 


galkal 


X Sp;px (w)e MP) dry drodw (9.77) 


Ponendo u = ty — 7] € v = to — 73 ed integrando per u e v avremo: 


Ruu(ri,r2,t1—t2) = 


n n 


+00 | 
=> f M $5 (t1) be (12) Hj (w) Spp (V) Hk lwe dv (9.78) 


j=l k=1" 7 


che si ottiene ricordando che l’ammettenza strutturale è la trasformata di Fourier della 
funzione di risposta ad un impulso unitario: 


+00 
Hj(w) F hyj(u)e“™"du (9.79) 


= oF bai 


Ricordando la 9.68 si ottiene: 


Spin) = | | da(eryoules)See(e1.22,w)dSas = jk(w) (9.80) 
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dove I;x(w) è il cross-spettro delle forze generalizzate proiettate sui modi j e k. 

Si può facilmente notare dalla relazione 9.78 che il secondo membro è, a me- 
no delle forme modali, la trasformata inversa di Fourier di H;(w)Ijx(V)H{(w). Si 
può dunque calcolare la densità spettrale della risposta che, considerando le for- 
zanti medie a parte, è per definizione la trasformata di Fourier della funzione di 
cross-correlazione Ruu(ri,ra,t1 — te): 


Suu(r1, Yo, w) = 5 5 $3 (01) Ok (t2) Hj;(w)Ijx(W)Hi(w) (9.81) 


j=l k=1 


A questo punto, la risposta risulta essere completamente caratterizzata. 


Risposta strutturale along-wind 


La turbolenza atmosferica produce un carico eolico fluttuante agente sulle strutture, 
a tal punto che queste iniziano a vibrare. Se le oscillazioni strutturali sono rilevan- 
ti allora la risposta media, ossia indotta dalla sola componente media della velocita 
del vento, é insufficiente per descrivere la risposta strutturale che di per sé evolve in 
campo dinamico. 

Dal momento che la turbolenza é descrivibile come un processo stocastico allora 
anche la risposta strutturale andra caratterizzata in termini statistici. In particolare, 
si assume che sia la turbolenza sia la risposta siano processi stazionari Gaussiani, 
caratterizzabili dunque dalla conoscenza dei primi due momenti statistici, media e 
funzione di correlazione. Si ricorda inoltre che grazie al gap spettrale la media è 
stazionaria su un intervallo di tempo compreso tra 10 minuti ed un’ora, lo stesso in- 
tervallo di tempo considerato dalle normative nella definizione del profilo di velocità 
media nel sito d’interesse. Pertanto, il campo di vento agente su una struttura lo si 
può scomporre in una componente media ed una fluttuante, dove la prima è responsa- 
bile della risposta statica e la seconda della risposta dinamica. I due effetti vengono 
successivamente sommati per individuare la risposta strutturale totale. 

In questo capitolo si descrive la procedura di validità generale per la valutazio- 
ne della risposta strutturale nella direzione del vento (along-wind), distinguendo tra 
strutture assimilabili a sistemi ad un grado di libertà e strutture complesse che invece 
richiedono lo studio dell’intero campo di pressione eolica e degli effetti che ne con- 
seguono. In quest’ultimo caso la risposta sarà dunque un processo multicorrelato e 
multidimensionale. 

Nel Cap. 9 abbiamo considerato la dinamica aleatoria dei sistemi a più gradi 


120 LEZIONI DI INGEGNERIA DEL VENTO 


di libertà. Tale trattazione è sempre valida soprattutto per strutture complesse il cui 
studio è condotto operando una discretizzazione mediante elementi finiti. Tuttavia 
in taluni casi è possibile studiare la struttura nel continuo mediante una trattazione 
analoga al caso discreto. 


10.1 Strutture riconducibili a sistemi ad un grado di libertà 


In questa sezione faremo riferimento alla trattazione di Davenport [2]. 

Quando la dimensione dei vortici del vento turbolento è molto maggiore della 
dimensione caratteristica della struttura in esame, allora quest’ultima può essere in- 
tesa come una struttura puntiforme. Entro certi limiti, la trattazione delle strutture 
puntiformi può essere estesa alle strutture con dimensioni maggiori. Ne deriva la 
possibilità di poter studiare tali strutture come sistemi ad un grado di libertà (Fig. 
10.1). 

La risposta dinamica di tali strutture può essere di tipo quasi-statico oppure ri- 


U+u 


>) 


we 


77, 77, 


Figura 10.1: Struttura intesa come sistema ad un solo grado di libertà 


sonante, in base alla rigidezza nella direzione del vento. Nel primo caso, la struttura 
è così rigida che le sue oscillazioni non interagiscono con la forzante applicata, men- 
tre nel secondo caso il moto della struttura induce un’amplificazione della risposta 
(risonanza), rispetto al caso precedente. 

Nelle suddette ipotesi, la struttura può essere intesa come una mensola con tutta 
la massa applicata sulla sua sommità ossia su un area A, pari all’intera area investita 
dal vento, posta ad un’altezza h rispetto al suolo (Fig. 10.1). Per quanto visto nel 
Cap. 9 avremo che la risposta dinamica della struttura in Fig. 10.1 sarà data da: 


See(w) = |H(w) |? Spr(w) (10.1) 
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con: 
1 1 


= mus (1 DI | 42 O) 
wi Wo 


con il consueto significato dei simboli, a parte la forzante eolica che viene definita 
come segue. 


(10.2) 


10.1.1 Definizione della forzante eolica 


x 


La forzante nella direzione del vento è espressa a partire dalla definizione della 
velocità del vento 


U(r,t)=[uu(r) + u(r, t)] -i+ v(r,t)-j+w(r,t)-k (10.3) 


dove y(r) è la velocità media e u(r, t), v(r,t) e w(r,t) sono le componenti flut- 
tuanti (turbolenza), agenti lungo le tre direzioni individuate da un sistema di rife- 
rimento Cartesiano x — y — z, nel generico punto individuato dal raggio vettore 
r=ax-i+ty-jtz-k. 

Definito con Cp(r) il coefficiente di pressione locale per unità di superficie, nella 
direzione del vento avremo: 


F(r,t) = 50C p(X) [uu(r) + u(r, t) + vlr, t) + wlr, t)? (10.4) 


dove A = 1 è l’area investita dal vento. 
Generalmente le componenti u?(r, t), v?(r, t) e w?(r, t) sono molto più piccole 
della velocità media e quindi possono essere trascurate, pertanto si ottiene: 


F(1,t) = 30Op(7) [uu (1) + ulr, t) + vlr, t) + wlr, OP ~ 


= 5eCp(X) [uo (1) + 2uy(r)u(r,t)] (10.5) 


L'equazione 10.5 può essere scomposta in una componente media ed una flut- 
tuante, 4. ;(r) e F}(r, t), rispettivamente: 


py(r) = suis ()Cp(r) (10.6) 
Fy (r,t) = pCp(r)uu(r)u(r,t) (10.7) 


Nel caso di una struttura riconducibile ad un sistema puntiforme possiamo ipo- 
tizzare che la forzante totale agisca su una massa posta ad una quota pari all’ altezza 
h della struttura, pertanto: 


1 
up = nsh) = 50 | wb )Ce)dA (10.8) 


Fi(h,t)= p | Car aa (10.9) 
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Inoltre, ipotizzando che il coefficiente di pressione non vari nello spazio, ossia Cp(r) = 
Cp(h) = Cp, si ottiene: 


1 
uf = 5 CDA (10.10) 
Fi(h,t) = puvCo | u(r,t)dA (10.11) 
A 


dove il coefficiente di pressione locale Cp assume il significato di coefficiente di 
trascinamento (drag), C'p, per unità di superficie, nella direzione del vento. 


10.1.2 Ammettenza aerodinamica 


Nel caso della struttura in esame la funzione di autocorrelazione si scrive: 
Rre,(h, T) = (puuCp)? | f E [u(rı, t)u(r2, t+ T)] dAdA = 
AJA 


(onCo)? | f Rus(risra add (10.12) 


Operando la trasformata di Fourier e sfruttandone la propietà di linearità, possiamo 
scrivere la densita di potenza spettrale della forzante: 


Sh) = (uCp)? | f Suler w)dAdA (10.13) 
AJA 


Inoltre ricordando la 4.29, possiamo scrivere: 


Sia (Et; r2, w) = V Suu(ti, w)Suu(12, w)Puulr1, r2, w) (10.14) 


che nel caso di una struttura puntiforme si riduce ulteriormente alla seguente espres- 
sione: 


Suu(11,12;0) = Suul h, w)Yuulr1, r2, w) (10.15) 


La 10.15 deriva dall’aver ipotizzato che la forzante turbolenta non vari nello spazio, 

ma solo nel tempo, pertanto considerando la forzante concentrata sulla sommita della 

struttura ad una certa quota h, avremo Syu(r1,0) = Suu(12,0) = Suulh, w). 
Pertanto la densità spettrale della forzante eolica espressa dalla 10.13 diviene: 


1 
Sr; (h, w) = (pCpp A)” Suho) f / WPuulr1, r2, w)dAdA (10.16) 
AJA 
che possiamo scrivere: 


Sr,r,(h,w) = (pCppy AY Suulh, w) x? (w) (10.17) 
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dove 
1 
x7 (w) = zÍ f Vuu(r1, r2, w)dAdA (10.18) 
A? AJA 


assume il significato di ammettenza aerodinamica. La funzione di coerenza è espres- 
sa dalla relazione 4.30 che viene riscritta di seguito: 


2 V (Cyry)? +(Cerz) 


Puulri, r2,w) =e FU (10.19) 


dove ry = |y1—y2| e rz = |z1 — 22] sono le distanze lungo y e z tra due punti sull’ area 
A, mentre Cy e Cz sono due costanti che garantiscono l’appropriato decadimento 
della funzione. 

Generalmente, l’ammettenza aerodinamica è espressa in funzione del rapporto 
tra la dimensione caratteristica della struttura, /, e la dimensione dei vortici, U/n, 
dove n = 5 è la frequenza. Fisicamente xL) può essere intesa come una misura 
della correlazione spaziale dei vortici: maggiore è la dimensione dei vortici rispetto 
a quella della struttura, maggiore è l’ammettenza aerodinamica, perché maggiore è 
la correlazione spaziale degli stessi vortici lungo la superficie investita dal vento; 
a parità di dimensione della struttura, poiché i vortici piccoli sono meno correlati 
nello spazio, l’ammettenza aerodinamica diminuisce (Fig. 10.2). Possiamo dunque 
affermare che per valori piccoli di nl il carico eolico, che in generale è una forzante 
multidimensionale e multivariata, diviene una forzante monodimensionale, puntuale, 
quasi completamente correlata, e data dalla risultante di tutte le pressioni locali nella 
direzione del vento agenti ortogonalmente alla faccia A, infatti LA) — 1. Per 
valori sempre maggiori del parametro ni la diminuzione della correlazione spaziale 
delle varie pressioni locali indotte dal vento diviene sempre più accentuata, pertanto 
xe (7) tende a valori sempre più piccoli dell’unità. 


10.1.3 Calcolo della risposta strutturale 


Dal momento che la forzante eolica è un processo stocastico, anche la risposta strut- 
turale avrà la stessa natura. Nell’ipotesi di forzante e risposta stazionarie e Gaussiane, 
vogliamo cercare di definire lo spostamento massimo della struttura. Si ricorda, che 
per valore massimo si intende il massimo valore atteso (valore estremo) della risposta, 
e non il massimo assoluto. Vogliamo dunque definire la quantità: 


Ee(h) = uelh) + geoe(h) (10.20) 


dove ue(h) è la risposta media, o¢(h) è la deviazione standard, e ge è il fattore di 
picco della risposta. 
La risposta media è data da: 


Lf Il 2 
= = — -pu CpA 10.21 
Hg ie ke HU D ( ) 
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nl 
U 
Figura 10.2: Ammettenza aerodinamica x? 


essendo ke la rigidezza alla traslazione nella direzione del vento. 
Per la valutazione della risposta dinamica, si consideri nuovamente la relazione 
9.81 che esprime la densità di potenza spettrale della risposta di un sistema continuo: 


See(h,w) = |H (w) Sr r (h, w) (10.22) 
da cui considerando la 10.17 si ottiene nel caso in esame: 
See(h,w) = (oC py AY |H (w)? Suu lh, wxw) (10.23) 


La varianza della risposta si ottiene integrando la 10.23, e nell’ipotesi di spettro 
unilatero avremo: 


+00 
oz(h) = (Conway? | | (w)|? Sulh, w)x?(w)dw (10.24) 


mentre il fattore di picco della risposta è definito secondo la seguente formula di 
validità generale per qualsiasi processo stazionario Gaussiano: 


2 
= 4/2 ln(vT) + Ucr (10.25) 
2ln(veT) (eT) 


dove T è la lunghezza della registrazione, mentre vg è la frequenza di attraversamento 
di soglia definita da: 


l Jo" ISS (10.26) 
27 N Se 


Ve = 
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10.1.4 Risposta quasi-statica 


Nel caso in cui la struttura è molto rigida l’ammettenza strutturale la si può approssi- 
mare come segue: 


Dy 2S = 
Hw)!" = |HO) = — = (10.27) 
pertanto la 10.24 diviene: 
CpuuA) rte i 
o2(h) = a J Suu(h, w)x2(w)dw (10.28) 
0 


che può essere riscritta, utilizzando la 10.10 e la varianza della turbolenza o?(h), nel 
seguente modo: 


Aus o2(h) [T9 Suu(h,w) 
ML Put w) 2 
o¢(h) = P f F(R) x (w)dw (10.29) 


Strutture puntiformi 


Come accennato precedentemente, se la dimensione dei vortici della componente 
fluttuante della velocità del vento è molto maggiore della struttura in esame, allora 
x?(w) — 1 pertanto la 10.29 diviene: 


Aut 02(h) (T° Suu(h 
Ruga PET f SII (10.30) 
k HU 0 o2(h) 
da cui introducendo l'intensità di turbolenza L,(h) = polti si ottiene: 
43 02 (h Api 
oe(h) = 3 ae 2 IBM) (10.31) 
by 


E dunque possibile trovare la deviazione standard della risposta: 
oe = “fy (h) (10.32) 


che sommata alla risposta media fornisce l’espressione della risposta totale (Eq. 
10.20): 


u 
&(h)= all + ge2Iv(h)] (10.33) 
La relazione 10.33 suggerisce la definizione di un carico statico equivalente: 


F.(h) = pp[1+gy,2Iu(h)] (10.34) 
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dove ge = gf, nell’ipotesi di risposta quasi statica. Inoltre si può introdurre il fattore 
di raffica, p, definito come segue: 


Fe(h 
o= FO a1 ap 2th (10.35) 


Possiamo quindi riscrivere la risposta totale: 


Ee(h) = Ty (10.36) 


da cui si evince che il fattore di raffica assume il significato di fattore di amplificazio- 
ne della risposta media. 


Strutture grandi 


Nel caso di strutture di dimensioni maggiori l’ammettenza, per quanto detto, diminui- 
sce dal momento che 1 vortici risultano meno correlati lungo la superficie strutturale 
investita dal vento. In questo caso si ottiene: 


so Suu h, 
oe = Lota (h) y [AP ewe = uva (0037) 


avendo posto: 


Salo) 
ke= Ni sy wdw (10.38) 
0 ou(h) 
In questo caso la risposta totale sara data dalla seguente relazione: 
u 
Ge(h) = 5 [1 + ge2Iu (h) Va] (10.39) 
Pertanto il carico statico equivalente assume la forma: 
Fe(h) = p;[1+ gf, 21u (h) v ke] (10.40) 
ed introducendo il fattore di raffica, y, 
Fe(h 
pe ” 1g dT (10.41) 


é possibile infine definire la risposta totale della struttura analogamente a quanto visto 
nella 10.36: 


Ee(h) = Ly (10.42) 


dove il fattore di raffica in questo caso tiene conto anche della correlazione spaziale 
della turbolenza. 
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10.1.5 Risposta in risonanza 


Nei paragrafi precedenti abbiamo calcolato la risposta di strutture rigide, per le qua- 
li il moto della struttura non influisce sulla risposta stessa (da cui appunto risposta 
quasi-statica). Nel caso di una struttura flessibile, dovremmo tenere in conto anche 
l'amplificazione in risonanza. Inoltre, se il moto della struttura è piuttosto pronun- 
ciato si generano ulteriori forze aerodinamiche. 

Conviene riscrivere l equazione di moto della struttura ricondotta ad un sistema 
ad un grado di libertà (Fig. 10.1): 


mé(h,t) + cé(h,t) + KÉ(h,t) = F(h,t) (10.43) 


In questo caso la forzante è data da: 


F(r,t) = 3POp(0)A [yur(n) + u(x,t) + vle, t) + (r,t) — é(e,)] © 
~ SeC(r)A [u8 (0) + 24u æu, t) — 2yv(e)É(1,1)] 10.44) 


che possiamo scomporre in una componente media jf, una fluttuante, F;(h, t), ed 
una componente aerodinamica, F,(h,t), applicata sulla sommità della struttura, alla 
quota h: 


1 

Lf = 5PHUCDA (10.45) 

F,(h,t) = puun | u(r, t)dA (10.46) 
A 

Fa(h,t) = puvé(h,t)CpA= caé(h,t) (10.47) 


La forza F,(h,t) è di natura viscosa, come la forza cÈ (h,t) nella 10.43. Inoltre, 
Fa(h,t) è anche di natura aerodinamica, indotta cioè dal flusso di vento inciden- 
te, pertanto Ca prende il nome di smorzamento aerodinamico. Possiamo dunque 
riscrivere l’equazione di moto nel seguente modo: 


mé(h,t) + (c+ Ca)E(h, t) F k£(h, t) = uf “e Fi(h, t) (10.48) 


Il problema consiste sempre nella valutazione del valore estremo della risposta 
(Eq. 10.20): 


Ee(h) = uelh) + geoe(h) (10.49) 


Procediamo quindi al calcolo dei vari contributi. 
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Risposta media 
La risposta media è data da: 


Mp _11 


= A 1 . 
a a Put OD (10.50) 


essendo ke la rigidezza alla traslazione nella direzione del vento. 


Risposta in risonanza 


Per quanto visto nel Cap. 9 la risposta della struttura, ricondotta ad un solo grado di 
libertà, è data da: 


See(w) = |H(w)|? Sr, r, (w) (10.51) 


dove però l’ ammettenza assume la seguente forma: 


1 1 
Hw)? = — 5 ; (10.52) 
meo (1 A 4 ay (2) | 
essendo: 
GF: Ca 
4 = an (10.53) 


il rapporto di smorzamento totale, comprensivo del contributo strutturale e di quello 
aerodinamico. Il coefficiente di smorzamento aerodinamico c, dipende dalla forma 
della struttura, ed è a favore di sicurezza se positivo (contribuisce ad aumentare lo 
smorzamento strutturale), mentre per valori negativi è a sfavore di sicurezza (contri- 
buisce a diminuire lo smorzamento strutturale). 

Nel caso di strutture flessibili, avremo: 


+00 
og (h) = (pC A)” f |H (w)? Suulh, w)x (7) dw (10.54) 


che può essere riscritta, ricorrendo alla 10.10 e alla varianza della turbolenza 02 (h), 
nel seguente modo: 


4u? oè(h) ft? Suulh, w) 
2 rS Tü 2 2 Puuh, 2 
oe(h) = 33 iz I k?|H(w)| 0) X (w) dw (10.55) 


da cui la quantità che cerchiamo, o¢(h), è data da: 


oelh) = Tf atu(h ff k2|H(w pee Bull) aes (10.56) 
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La risposta dinamica può essere calcolata come la somma della componente quasi 
statica, che si ottiene per |H(w)|? ~ |H(0)|? = 1/k?, e di quella dinamica: 


a Suu(h, w) 
2 2 ’ 2 = 


dove: 


too Suulh,w t Suulh,w 
b= fori eda A) 


alh) oalh) 
(10.58) 
TE Suu(h, wo) WoSuu(h, wo) T 
kr a 2 H 2 Fuu 3 (0) 2 g d = on uu bj oO 2 i 
(RR. ede I oo) 
(10.59) 
Pertanto, la risposta totale della struttura sarà data da: 
&(h) = a E + g¢2Iy(h) ky + kr (10.60) 
ed il carico statico equivalente assume la forma: 
Fe(h) = well + gf 2Iu(h)v ko + kr] (10.61) 
ed introducendo il fattore di raffica, y, 
Fe(h 
= A =1+ gp,2Iu(h) WV kp + kr (10.62) 
f 


é possibile infine definire la risposta totale della struttura analogamente a quanto visto 
nel caso delle strutture rigide: 
u 
ée(h) = Fo (10.63) 
con © che tiene conto della correlazione spaziale della turbolenza, della componente 
quasi-statica e di quella dinamica. 


10.1.6 Strutture puntiformi 


Nel caso di strutture puntiformi si ha x?(w) — 1 e quindi si ottiene: 


+00 EOS 
ky = I p (0) Pee) -1 dw = I Sul hw) | l dw=1 
0 0 


o%(h) oCh) 
(10.64) 
too Suulh, wo) WoSuu(h, Wo) T 
kp = kH 2 UUN OOD o g a OR UMN O E, 10. 
j PROP e do e (069 
pertanto il fattore di raffica assume la forma: 
Felh 
p= Fe(h) = 1+912Iu(h)V14+kr (10.66) 


Mf 
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10.2 Strutture allungate nella direzione del vento 


Nei precedenti paragrafi abbiamo trattato il caso di strutture a prevalente sviluppo in 
una direzione (e.g. lungo z per la struttura in Fig. 10.1). Con riferimento alla Fig. 
10.3, consideriamo ora il caso tridimensionale, ossia anche la dimensione lungo la 
direzione del vento. In questo caso la forzante di trascinamento sarà data da due con- 


Z 


X 0.0 25.0 50.0 75.0 100.0 125.0 150.0 
Time {s} 


Figura 10.3: Struttura estesa nella direzione del vento 


tributi, la componente agente sulla faccia direttamente investita dal vento (windward 
wall) e da quella che agisce in scia (leeward wall). Possiamo sempre scomporre il 
carico eolico in una forza statica ed una dinamica: 


1 1 
Hf = bp, + Ufi = 5 PHU CDy Aw + srt CDA (10.67) 


Filh, t) = Fw(h,t) + Filh, t) = 
= puvC. f u(r, t)dAy + puuCD, si u(r,t)dA; (10.68) 
Aw Aj 


Il calcolo della risposta media é analogo a quanto visto precedentemente, mentre 
per il calcolo della risposta dinamica conviene fare alcune osservazioni. In parti- 
colare, siamo interessati al calcolo della densita spettrale della forzante dinamica 
F,(h,t) = Fw(h,t) + Fi(h, t) che sappiamo essere la trasformata di Fourier della 
sua funzione di autocorrelazione: 


Rap (ht) = E [F;(h, t)F;(h,t + 7)] = 
= E{[F,(h,t) + Fi(h,t)] [Folh,t +7) + Fi(h,t+7)]}} = 
= E|F,,(h,t)Fo(h,t+7)] + E [Fi(h, t)Fi(h,t + 7)] + 
+ E[F,(h,t)F(h,t+7)] + E[Fy(h,t+7)F(h,t)] = 
= Rep, (ht) + Rag (h, T) + Rew (h, T) + Rrr (h,-T) (10.69) 
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tenendo presente che Rrr „(h,T) = Rr, w(h,—T). Operando la trasformata di 
Fourier si ottiene dunque la densita spettrale di potenza del carico eolico: 


Srp (h, w) =S pP (h, w) + Srr(h, w) + Sr, F;(h, w) + Sr r (h, w) (10.70) 


che sarà data dunque dalla somma degli auto-spettri e dei cross-spettri delle singole 
componenti. 

Generalmente, come si vede dalla Fig. 10.3, le due componenti della forzante so- 
no poco correlate tra loro, pertanto Rp „r (h, +7) — 0 con buona approssimazione. 
Conseguentemente Sr, r (h,w) > 0 e Sh r (h,w) — 0, quindi la densità spettrale 
di potenza sarà data da: 


Sr,r,(h,w) = Spr, (h,w) + Srp (h, w) (10.71) 


Noto lo spettro della forzante eolica si può procedere al calcolo della risposta 
come visto nei precedenti paragrafi, distinguendo tra struttura rigida o flessibile, 
puntiforme o estesa. 


10.3 Trattazione generale della risposta strutturale along- 
wind 

Nei precedenti paragrafi abbiamo considerato il calcolo semplificato della risposta 
along-wind, considerando la struttura come un sistema ad un grado di libertà. In que- 
sta sezione viene descritta la trattazione, di carattere più generale, per il calcolo della 
risposta strutturale nella direzione del vento, con particolare riguardo alle strutture 
line-like e plate-like (Dyrbye and Hansen [4]). 

Tale procedura si riferisce a strutture aventi forme semplici, riconducibili alle 
seguenti tipologie: 


e edifici a torre; 
e ciminiere; 
e ponti a mensola; 


e linee di trasmissione elettrica; 


opere per l'illuminazione. 

Inoltre, si assumono le seguenti ipotesi: 
e Ilcarico eolico è determinato dal campo di vento indisturbato; 
e La struttura è elastica-lineare con smorzamento viscoso; 


e I modi propri di vibrare nella direzione del vento sono disaccoppiati tra loro. 
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10.3.1 Equazione di moto per le strutture che vibrano nella direzione 
del vento 


L'equazione di moto della struttura che vibra nella direzione del vento è la seguente: 


m(r)€(r, t) + c(r)£(r,t) + L[E(r)] = F(x, t) (10.72) 


con il consueto significato dei simboli. In analogia a quanto visto nel paragrafo 9.5 è 
possibile decomporre il problema proiettandolo nello spazio delle coordinate modali, 


w(x): 

= Ís w(x) F(x, t)dS _ pit) 
Js U2 (2) m(r)dS Mi 

Se più modi di vibrare nella direzione del vento contribuiscono in maniera signifi- 


cativa alla risposta strutturale, allora avremo che la risposta dinamica totale, R, sarà 
data da: 


oh=) oR: (10.74) 


di + Wiwids + w? qi i=1,2,...,+00 (10.73) 


10.3.2 Funzioni d’influenza 


Finora abbiamo considerato come risposta strutturale lo spostamento estremo, tutta- 
via conviene formulare il problema in modo da poter estendere la valutazione del- 
la risposta in funzione di qualsivoglia quantità (spostamenti, deformazioni, tensioni 
e caratteristiche di sollecitazione). Peranto, introduciamo le funzioni d’influenza, 
Ig(r), della generica risposta R. 

La funzione di influenza è la generica risposta che si ottiene in un generico pun- 
to della struttura per effetto di un singolo carico unitario applicato in un altro punto 
della struttura stessa. 

Calcoliamo, ad esempio, il valore medio del taglio alla base di una struttura ver- 
ticale assimilabile ad una mensola. Per ogni singolo carico unitario applicato alla 
quota z, ossia F(z) = 1, avremo che il taglio alla base è uguale, ma di verso con- 
trario, al carico applicato, T(0) = F(z) = 1. Pertanto, in questo caso la funzione 
d'influenza del taglio, Ir(z), è costantemente pari ad 1. Infatti, se si applica una 
forza generica alla quota z, sia essa F(z), avremo T (0) = 1- F(z), con I7(z) = 1. 

Nel caso del momento flettente, avremo M (0) = z - F(z), per qualsiasi carico 
F(z) applicato alla generica quota z. Conseguentemente, essendo la funzione d’in- 
fluenza la risposta generica dovuta ad un carico unitario, avremo Iyy(z) = z- 1 = z. 

In generale, per un carico non uniformemente distribuito, il momento flettente 
alla base sarà dato dalla seguente relazione: 


h h 
mo) = | Iu(2)F(2\dz f ZF (z)dz (10.75) 
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E facile dimostrare che nel caso di un carico uniformemente distribuito si ha: 


2 


(10.76) 


h h 
mo) = | uF = F | de = PE 


ossia il ben noto momento d’incastro di una mensola caricata uniformemente. 


10.3.3 Valore estremo della risposta 


Lo scopo di questa trattazione è quello di valutare il valore estremo della risposta: 


Re = UR + 9ROR (10.77) 


or= Jo, + OR, (10.78) 


essendo oR, € oR, la deviazione standard della risposta quasi-statica e in risonanza, 
rispettivamente. 


dove: 


10.3.4 Strutture Line-Like 


Per strutture Line-Like s’intendono strutture con prevalente sviluppo in una direzio- 
ne. Assumendo che il vento agisce in direzione x, possiamo distinguere tra strutture 
che si estendono prevalentemente in direzione y (e.g., ponti) ed altre che si estendono 
in direzione z (e.g., ciminiere). 

In questa sezione consideriamo strutture che si estendono in altezza, dunque lun- 
go z, ma la trattazione è del tutto analoga per le strutture allungate nella direzione 


y. 


Equazione di moto 


L'equazione di moto, associata all’i-esimo modo di vibrare, sarà dunque data dalla 
seguente relazione: 


h 
è aa Sh Vi(2)F(2,t)dz _ pilt) 
Git 2Viwidi + w; q IEZIENTOTE mi 


i=1,2,...,+00 (10.79) 


Definizione della forzante eolica 


La forzante eolica è data da: 


F(2,t) = 5eC(2)D(2) [pu (2) + (2,8) + olet) + wlt) -Eld © 


= 50c(2)D(2) (A) + Quy (z)u(z,t) — 2uu (2 Elz, 5) (10.80) 
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che possiamo scomporre in una componente media ju p(z), una fluttuante, F;(z, t), ed 
una componente aerodinamica, F’,(z, t), 


up(2) = put; (2)C(2)D(2) (10.81) 
z)u(z,t) (10.82) 
2)é(2,t) = pC(z)D(z) uu (z)divi(z) (10.83) 
dove nella 10.83 si è tenuto conto della decomposizione modale che porta alla re- 


lazione 10.79. La quantita D(z) è la dimensione della struttura nella direzione 
perpendicolare alla direzione del vento, mentre C'(2) è un coefficiente di forma. 


Smorzamento aerodinamico 


La forzante definita nella 10.83 è di natura aerodinamica, in quanto indotta dal vento 
incidente, ma anche di natura viscosa perché è proporzionale a q. Pertanto, è possibile 
definire uno smorzamento aerodinamico: consideriamo la 10.79 dove la forzante sia 
solo quella di natura aerodinamica: 


So UR (2)pC (2) D(z) uy (2) dz 


Gi + 2viwidi + w2G = di (10.84) 
Se ee) m(z)dz 
ponendo 
oy, wo — te VieC@P@)mo(2)dz (10.85) 
Di fo v2(2)m(a)dz 
avremo: 
Gi + 2(0; + Va iwigi + wi gi = 0 (10.86) 


pertanto, l’ammettenza strutturale associata all’i-esimo modo di vibrare sarà: 


1 1 
4 


mu 3\2 2 
0-3) +42 (2)) 


con Vii = (vi + Vai). 
E possibile scrivere il rapporto di smorzamento aerodinamico 1,,; nel seguente 
modo: 


|H (wi)? = (10.87) 


Vag = DOPOD (A)de _ 1 oy Uri 
n FE AmAde ae Mae 


(10.88) 
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dove: 

ee 

Urid.i = niD(2*) (10.89) 
— _ Mgalh 

Mrid,i and pD2(z*) (10.90) 

h 2 
Mgi = I m(z) ante (10.91) 
1 (È C(2)D(2)uu(2)b?(2) 

a,i 7 1 d 10.92 
de ii CE) DE) pv (EVO) ae 


dove z* è una quota di riferimento arbitraria, n; = w;/(27) è l’i-esima frequenza 
naturale del sistema, e h è l’altezza della struttura. 

Complessivamente, possiamo scrivere l’equazione di moto associata al i-esimo 
modo di vibrare nel seguente modo: 


Role) + F(z bd 


(10.93) 
JÈ w?(z)m(z)dz 


Gi + 24, wid; + Ww? Gi 


Risposta media 


Detta Ip(z) la funzione d’influenza della generica risposta possiamo calcolare la 
risposta media come segue: 


h 
m=) Ig(z)up(z)dz (10.94) 


che possiamo scrivere, facendo riferimento ad una generica quota z*, come segue: 


uR = ZPhD(2*)O (4) pt (2 IR") Ym (10.95) 
dove: 
1 h 
Yn =F Im(z)dz (10.96) 
h Jo 
DACIA) 
Im) = DECENTI) eo 


La funzione gm(z) è una funzione adimensionale che esprime la variazione del 
carico medio e della funzione d’influenza della generica risposta. Pertanto, Ym rap- 
presenta l’effetto complessivo di gm(z) sull’intera struttura. 
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Risposta quasi-statica 
La risposta quasi-statica é data da: 
ch, = Rrr (T = 0) (10.98) 


dove Rr, R, (7) è la funzione di autocorrelazione della generica risposta, Ry(t), defi- 
nita come segue: 


h 
rap) = f Ip(z)Fi(z, thdz (10.99) 
0 


La funzione di autocorrelazione di Ry(t) è dunque data dalla seguente relazione: 
RryBy(T =f fa E [Fi (a1, t)Fi(22, t + 7)] Ir(21)Ir(22)d21d22 = 


= [ ic PO) AE [oC (22) D (22) tur 22)] 
x E [u (21, t)u (22, t+ T)| Ip(21)Ir(z2)dzidze = 


= [ i [pC (21) D(z1) uu (41)] [eC (22) D(z2) wu (z2) à 
x Ruu(21; 22, DIR(z1)Ir(22)dz1dz9 (10.100) 


Dal momento che stiamo considerando la sola componente fluttuante, allora oh, = 
RRR, (T = 0) ossia: 


h ph 
ab = ff POED PCI 
0 JO 
x Ruul(z, 0) LR (41) Ir (22) dz1dz2 (10.101) 


dove si è supposta che la funzione di cross-correlazione della componente fluttuante 
della velocità del vento dipenda dalla distanza r, = |z1 — 22]. Inoltre, moltiplicando 
e dividendo per 0,,(21) (22), e considerando una quota di riferimento z* possiamo 
scrivere: 


* 


2°) D(2°) pu (2*)ou(2*)Ir(2*))? x 
«Lf [eC (z1).D (21) uu (21)] [eC (z2).D (22) uu (22)] 
[pC (2*) D(z*) pu (2*)ou(2*)In(2*)]? ou(21)ou(z2) 
x Ruu(7::0)Ir(21)Ir(22)0u(21)0u(22)dz1dz9 = 
= [POD wu (TREE x 
ee [pC (21) D (21) mu (21) [eC (22) D (22) mu (22)] puulrz) 
[pC (2*) D(z*) nu (2*)ou(2*)In(2*))? 
x Ir(z1)Ir(22)0u(71)0u(22)dz1dz2 (10.102) 


X 
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dove p(rz) = —Buu(r=.9)_ è Ja funzione di cross-correlazione normalizzata della com- 
Tu (21)0u(22) 


ponente fluttuante del vento. 
L'equazione può essere riscritta nel seguente modo: 


o, = [POD wu (2)u(2) PREI (10.103) 
dove: 
1 h h 
Js DI f Pulrz)g(21)%(z2)dz1dz2 (10.104) 
C(2)D(z)o(z)uu(z)Ir(2) 


(10.105) 


9) = ADE E Ee 


La funzione 9p() esprime la variazione del carico fluttuante lungo la struttura, 
mentre Je è la varianza adimensionale della risposta e può essere scomposta in due 
integrali: 


1 h 
=> | klrz)o(r2)are (10.106) 
0 
2 h—-rz 
k(r.) = | 9o(2) 96(2 + r2)dz (10.107) 
0 
I valori asintotici di J? sono i seguenti: 
h 
JÈ > k(0)/6- per $: = 77 > +00 (10.108) 
h 
BR per be = 77 70 (10.109) 
dove 
1 h 
Yb = a g(z)dz (10.110) 
h 0 


fornisce la deviazione standard adimensionale della risposta, Jp, nel caso di turbolen- 
za completamente correlata, p,(r;) = 1. 

Nel caso in cui le funzioni gm € gẹ abbiano segno costante la varianza della 
risposta quasi-statica può essere scritta nella seguente forma: 


oh, = [ur2I,(2*)]? OF ky (10.111) 


dove [,(z*) è l’intensità di turbolenza alla quota di riferimento z*, mentre 


0, = (10.112) 
Ym 
J2 

ky = > (10.113) 
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Una buona approssimazione di ky è data da: 


1 


CSI 


(10.114) 


dove I° è una costante che dipende da dalla funzione g,(z). In Fig. 10.4 è rappresen- 
tata la J? per alcune funzioni g,(z). 


0.1 


Accurale expression 


Approximale expression 


0.01 


Non-dimensional response varianca, J? 


0.01 0.1 1 10 100 


Figura 10.4: J? in funzione di ¢, nel caso di g(z) = 1 Œ = 1/2); go(z) = z/h 
C = 3/8); a(z) = sin(rz/h) Œ = 4/7?) 


Nel caso in cui gm € gẹ cambiano di segno lungo la struttura, kẹ deve essere calco- 
lata caso per caso. In Fig. 10.5 è riportato l'andamento della varianza adimensionale 
della risposta per gp(z) = 2z/h — 1. 


Risposta in risonanza 


In questo caso non è possibile procedere direttamente con l’ausilio delle funzioni 
d’influenza, dal momento che queste sono modificate dal moto della struttura che di 
per sé è l’incognita del problema. Pertanto, bisogna dapprima risolvere il problema 
dinamico, e successivamente definire un carico dinamico inerziale, ossia proporzio- 
nale all’ accelerazione della struttura, che permette l’utilizzo delle funzioni d’influen- 
za. 

Siamo alla ricerca di una forzante F(,t) definita nel seguente modo: 


F4(z,t) = m(z)Q(z,t) = m(z) De pil ült) (10.115) 


dove m(z) è la massa per unità di lunghezza, Q(z, t), è l'accelerazione totale della 
struttura che può essere decomposta, come più volte menzionato, nella somma delle 
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2 
d 


J, 


0.1 


Non-dimensional response variance, 


Figura 10.5: J? in funzione di œ, nel caso di gy(2) = 2z/h — 1. Una buona 
approssimazione è data da .J}($.) = HI 


accelerazioni associate ai sistemi elementari in cui è possibile scomporre la struttura. 
Le funzioni +;(2) sono dunque le forme modali della struttura. 

Consideriamo l’i-esimo modo di vibrare della struttura, allora l’i-esimo carico 
dinamico sarà!: 


Faalz, t) = mepi ült) = m(z)w? di(z)ai(t) (10.116) 


che consente la valutazione della generica risposta in risonanza Ry; con l’ausilio 
della funzione d’influenza Ip(z): 


h 
Rei) = f Ip(2)Fai(2,t)dz (10.117) 


La risposta in risonanza si riconduce al calcolo della varianza della generica risposta 
in risonanza: 


+00 
Oh, = f SR,.iRri(7)dw (10.118) 


dove Sr,.R,;(w) è lo spettro della generica risposta in risonanza, R, (t), definito 
come: 


SR,iRri (w) = 
h rh 
= / . m(z1)m(z2) Wi (21) vi(22)wh In(21) Ir (22)dz1dz2| Sqgi(6) 
0 JO 
(10.119) 
'Il carico dinamico è definito a meno del segno negativo che deriva dal fatto che ġ; (w) = —w? q(t). 


Questo perché il segno è ininfluente sulla valutazione della risposta. 
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Dobbiamo pertanto calcolare 5,4; (w) risolvendo il problema dinamico associato 
all’i-esimo modo di vibrare della struttura: 
h 
Gi + 2 wii + wi gi = hi i (10.120) 
J we(z)m(z)dz 


dove 
Filz, t) = pC(z) D(z) pu (z)u(z, t) (10.121) 


Nell’ ipotesi di modi disaccoppiati avremo che, per l’i-esimo modo di vibrare, la 
densità spettrale della risposta è data da: 

Saqi (0) = |Hi(w)]?Ki() (10.122) 

dove l’ammettenza strutturale | H;(w)|? tiene conto del coefficiente di smorzamen- 


to totale, v: i, comprensivo del coefficiente di smorzamento aerodinamico, va,i (Eq. 
10.88). 


1 1 
m2w4 32 2 
“Gae rea 


Il cross-spettro delle forze modali associato all’i-esimo modo di vibrare è dato dalla 
seguente relazione: 


|H(w)? = (10.123) 


h h 
Is) = f | Wi (21) Ui (22) Spr(41, 22,w)dzdzo (10.124) 


ottenuta calcolando la funzione di autocorrelazione del carico i Wi(z) Fi(z, tdze 
successivamente operando la trasformata di Fourier. 

Tenendo presente la 10.121 il cross-spettro del carico può essere espresso nel 
seguente modo: 


Srr(21, 22,7) = 
= [pC(a1)D(a1) Mu (21)] [eC (Z2) D(z2) uu (22)| Suu (21, 22,60) = 
= [pC(21) D(a) mu (41) [PC (22) D(z2) Hu (22)] x 
X VJ Suu(Z1,W)Suu(Z2,W)Puu(rz,w) (10.125) 


dove r, = |z1 — za], € Yuulrz,w) è il co-spettro normalizzato della componente 
longitudinale della velocità del vento. 

Considerando una quota di riferimento z* possiamo riscrivere la densità spettrale 
della risposta associata all’i-esimo modo di vibrare nel seguente modo: 


Sgaile) = [hpC(z*) D(z") uv (A di (291° |H(wi) |? IL (0) Suu (2*) (10.126) 
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dove J? ;(w) assume il significato di ammettenza congiunta e può essere scritta nel 
seguente modo: 


1 sh ph 
Pie) = a f Ni Gr,i(21, 0) Gri( 22, W) Puu(Tz, w)dzıdz2 (10.127) 


essendo da w) la funzione che tiene conto della distribuzione del carico eolico in 
condizioni di risonanza: 


C(2)D(2)uu(2)V Suu(2,0)wi(2) 
C(2*)D(2*)uu(2*)V Suu(2*w)wi(2*) 


L’ammettenza congiunta J 2 (w) può essere calcolata scomponendo l integrale 
nel seguente modo: 


Gr,i(Z,W) = (10.128) 


ir 
J} iw) = ral kz ilrz, w) puulrz,w)drz (10.129) 


2 h-rz 
kz i(rz,w) = a Or,i(Z, W) Gr,i(Z + Tz, w)dz (10.130) 


ce 

= 

2 

Ss 0 

® 

o 

= 

w 

a 

0.01 

= ; : 

2 —— Exponential function 

Modified exponential function 

0.001 


Figura 10.6: J2(w) in funzione di ¢, nel caso di (z) = 1, variando il co-spettro 
normalizzato. (—) J?(w) ottenuta utilizzando wu 


Varianza della risposta in risonanza 


Risolto il problema dinamico, ossia nota la densità spettrale S7,g; (7), è possibile 
calcolare la varianza della generica risposta in risonanza associata all’i-esimo modo 
di vibrare, oi , secondo la 10.118. 

Si può dimostrare che 


Ong = Reta") 0 ins (10.131) 
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Nel caso di modi propri di vibrare che non cambiano segno si ha: 


; h . 
Ym Mg, JO wi(z*) In(2*) 
T? Wi Syy(2*, wi) 
kri = li ee Is i) I 10.1 
ina Re) E on 
con 
il h 
ri = > si Gri(Z, O)dz (10.134) 
h 0 
h 2 
pi (z) 
Mg i = m(z)—4 dz (10.135) 


Quando in modi propri di vibrare cambiano segno, la varianza della risposta deve 
essere calcolata caso per caso. Ciò è dovuto alla diversa interazione tra la turbolenza 
e i modi di vibrare della struttura. 

Nel caso in cui la turbolenza è costituita prevalentemente da vortici piccoli, l’ am- 
mettenza congiunta tende al suo valore asintotico, 


iL 
Te). (10.136) 
A ) Td. 
dove J, è una costante che dipende dalla funzione g,(z), dz = Cambi, essendo C, la 


costante che garantisce il decadimento di Yuu(rz, w) lungo z. 

Quando invece la turbolenza è costituita da vortici di dimensioni paragonabili alla 
struttura, se non maggiori, allora l’interazione tra la turbolenza e la struttura è così 
diversa che il co-spettro normalizzato %,u(rz,w) non tende più all’unità alle basse 
frequenze. Questo fenomeno è maggiormente pronunciato per strutture che vibrano 
con modi non aventi segno costante. 

Un esempio di come varia l’ammettenza congiunta relativamente ad un modo di 
vibrare che cambia segno, è dato dalla Fig. 10.7 


Risposta totale 


Generalmente la risposta in risonanza è principalmente data dal primo modo di vi- 
brare. Pertanto, tenendo presente la 10.77, la 10.94, la 10.111 e la 10.131, possiamo 
scrivere la risposta totale come: 


Re = UR h + gr2Tu(2*)4/02kp + DI (10.137) 


dove la risposta dinamica è quella associata al solo primo modo di vibrare. 
Infine, possiamo quindi scrivere il coefficiente di amplificazione della risposta 
media, coefficiente di raffica, nel seguente modo: 


p= € £1 + g,21,(2)/028,+ 02%, (10.138) 
LR 
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Joint acceptance lunclion. 
o 
= 


Exponeniial function 
nn Modilied exponential function 
TTT a 


a a | Tm] 


0.01 0.1 1 10 100 


Figura 10.7: .J?(w) in funzione di $. nel caso di (z) con segno non costante, 
variando il co-spettro normalizzato. (-) J2(w) ottenuta utilizzando uu 


Strutture line-like orizzontali 


Nel caso di strutture Line-Like orizzontali (e.g., ponti), che si estendono cioè in 
direzione y per una lunghezza b, avremo: 


1 fè 
Til) = al ky i(ry, W)PuulTy, w)dry (10.139) 
2 b-ry 
kyi w) = 3 | Iri(Y:W)gyi(Y + ry, w)dy (10.140) 


dove in questo caso 


C(y)D(ly)uu (y) V/Suuly,w)vily) _  Cly)D(y)vily) 
C(y*) D(y*) uu (y*) Suly w)vily*)  Cly*)D(y*)vily*) 
(10.141) 


Gr ily, w) = 


essendo py (y) = pu (y*) e Suu(y, w) = Suu(y*, w). La trattazione è uguale a quella 
delle strutture Line-Like verticali. 


10.3.5 Strutture plate-like 


Le strutture Plate-Like sono strutture che si estendono sia in direzione verticale, z, 
sia in direzione orizontale, y. La trattazione é analoga a quella delle strutture Line- 
Like, pertanto si riportano direttamente le quantita d’interesse per la valutazione della 
risposta estrema: 


Re = UR h + Gp2Tu(2*) / 02k + DI (10.142) 
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dove la risposta dinamica è quella associata al solo primo modo di vibrare, che ge- 
neralmente prevale sugli altri. Il coefficiente di raffica è definito secondo la nota 
relazione: 


ga 2149212) 0218 + Oe, (10.143) 
HR 


Vediamo in dettaglio i vari contributi della risposta, ossia la risposta media, quasi- 
statica e in risonanza. Quest’ultima è calcolata tenendo in conto soltanto il primo 
modo di vibrare. 


Risposta media 


1 
LR = SPORCARE Y") Ym (10.144) 


con: 


1 h prb 
Im = af Ni Im(y, 2)dz (10.145) 
bh Jo Jo 


Cz, Yu (2) IR(z, y) 
C(2*, y*) ut (2*)Ir(2*, y*) 
dove h e b sono rispettivamente le dimensioni lungo z e y della struttura, mentre 


Ir(x, y) è la funzione (superficie) d’influenza della generica risposta R. Il punto di 
coordinate (z*, y*) è un punto di riferimento della struttura. 


(10.146) 


Im(£, y) = 


Risposta quasi-statica 
La varianza della risposta quasi-statica è data dalla seguente relazione: 

oh, = [pbhC(z*, y*) uy (2*)ou(z*))? 2 (2*, y*) JP (10.147) 
dove: 


1 h h b b 
ioe. I f i J Pul" z, y)go(21, Y1)9(22; v2)dzıdz2dyıdyz (10.148) 
bh? Jo Jo Jo Jo 


C(z,y)o(z)uu(z)Ir(z,y) 


glz, y) = — i oe (10.149) 
ED = TEMA MEET 
Se le funzioni gm € gp non cambiano segno, avremo: 
g=® (10.150) 
Ym 
2 
ky = eb (10.151) 


N 
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dove: 


h b 
W= aa | f mlz, y)dydz (10.152) 
Dh Jo Jo ( 


che sarebbe il valore di J? quando p,,(rz,7ry) = 1. Assumendo go(z,y) = go(z) go(y) 
il fattore di risposta k, è ben approssimato dalla seguente formula: 


1 


b= 5 3 z (10.153) 
1+ Cydy)? FEAS + CGn 
h b 
oz = Tz dy = DI (10.154) 


dove I è una costante che dipende dalla funzione gp. Alcuni valori sono riportati in 
Fig. 10.4. 


Risposta in risonanza 


Per la risposta in risonanza l’ammettenza congiunta è la seguente: 


1 
2 
Jey) = papa 
e in) nG Ab 
dl f | forti) ge (castra te) erro 
o Jo Jo Jo 
(10.155) 
dove: 
9r(2,4;6) = 239) 12) y Suul% yw) P(% y) (10.156) 


ae „Y 2*) y Suu(z*, y*, w)v(2*, y*) 


Per il calcolo della risposta in risonanza si hanno le seguenti quantita associate al 
primo modo di vibrare con segno costante: 


h pb ie Ir, 
ZIE Jo Jo i Tae pe dude 
0, a (10.157) 
Ym f o (2: y) gagy Ide 
2 
(E: T w1 Suul?” LAA TA (10.158) 


Any, =a 2(2*) 


dove w = 27n, è la frequenza angolare associata al primo modo di vibrare, m(x, y) 
è la massa per unità di superficie, mentre yr è data da: 


1 h pb 
-5 [accu 20)duas (10.159) 
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Infine, K,(w) è una funzione di riduzione che tiene conto dell’ interazione della tur- 
bolenza con il modo proprio di vibrare della struttura. Nel caso di modi di vibrare con 
segno costante si ha che K(w) tende all’ unita alle basse frequenze in corrispondenza 
delle quali la turbolenza è maggiormente correlata, mentre decresce progressivamen- 
te alle alte frequenza dal momento che i vortici sempre più piccoli sono sempre meno 
correlati. Un’approssimazione della funzione di riduzione Ks(w) è così data: 


Ks(w)= 2 : 2 2 me) 
1+ y (Gupp)? + (G.0. + (2Gy6yG-42) 
A Chw B Cybw 
= 7e) #7 TE) ale 


dove Cy e C, sono le costanti della 4.30, mentre G; e Gy sono delle costanti che 
dipendono dalla funzione gr. 

Nel caso di modi di vibrare che cambiano segno, la risposta in risonanza deve 
essere calcolata caso per caso, visto che l’interazione tra turbolenza e modo di vibra- 
re induce un cambiamento della funzione di ammettenza congiunta come discusso 
precedentemente. 


yz.) =1 


o 


o 
o 


Size reduction function, Ke 


2001 


95% 


D{y.2.0) = zm 


o 


ini 
ici 


Size reduction function, Kg 


0.001 


Figura 10.8: Funzione di riduzione K, al variare di g(z,y,w), dy e dz 


Risposta strutturale along-wind secondo 
l’Eurocodice 1 


11.1 Scopo 


L’ECI fornisce regole e metodi per il calcolo delle azioni del vento agenti su strutture 
che non superino un’altezza di 200m, e sui loro componenti strutturali; su ponti, 
viadotti stradali e ferroviari che si estendono fino a 200m di lunghezza; su passerelle 
e ponti pedonali fino a 30m. 

I carichi del vento devono essere calcolati per l’intera struttura, o su parti di essa 
(componenti, rivestimenti ed elementi di fissaggio, etc.). Le ciminiere e le strutture 
a mensola sono contemplate, ma non vi sono requisiti fondamentali per le torri a 
traliccio. Ponti strallati e ponti sospesi non rientrano nella normativa, e dovrebbero 
essere oggetto di studi specialistici. Non vi sono regole per le antenne strallate e per 
le strutture offshore. 


11.2 Modellazione dell’azioni del vento 


Le azioni del vento sono rappresentate in termini di pressione o forza. Le pressioni 
indotte dal vento s’intendo agenti perpendicolarmente alla superfici investite salvo 
diverse indicazioni come nel caso di forze tangenziali dovute al vento radente. 

Si fa riferimento ai seguenti parametri: 


© Qref intesa come la pressione eolica di riferimento; 
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e ce(z) ossia il coefficiente di esposizione che tiene conto dell’ effetto della ru- 
gosita superficiale del sito sulla velocita del vento alla quota z. In particola- 
re ce(z) modifica la pressione media in una pressione di picco (estrema) per 
tenere conto della turbolenza atmosferica. 


e cq è il coefficiente dinamico che tiene conto della correlazione del vento tur- 
bolento agente sulla struttura, e dell’amplificazione della risposta in risonanza 
dovuta al moto della struttura. 


11.3 Pressione dinamica di riferimento 


La pressione dinamica di riferimento è definita come segue: 


1 


dref = 5PU res (11.1) 


dove U,ef è la velocità eolica di riferimento che dipende dal sito in questione, e p è 
la densità dell’aria pari 1.25K g/m. 

La velocità eolica di riferimento è ottenuta operando una media della componente 
longitudinale della velocità del vento su un’intervallo di 10 minuti (all’interno del 
gap spettrale), ad un’altezza di 10 metri sopra un suolo appartenente alla seconda 
categoria. Tale velocità di riferimento è intesa come la velocità che ha una probabilità 
di eccedenza pari a 0.02, corrispondente ad un tempo di ritorno di 50 anni. 

In particolare si ha: 


Uref = CdirCtemCaltUref,0 (11.2) 


dove cq;r è il fattore di direzione posto pari a 1 se non altrimenti specificato; cem è il 
fattore temporaneo pari ad 1 se non altrimenti specificato; cj; è il fattore di altitudine 
posto pari a 1 se non altrimenti specificato, U,.e;,0 è il valore di base della velocità di 
riferimento del vento. 


11.4 Profilo della velocità del vento 


La velocità media del vento ad una generica quota z al di sopra del suolo dipende 
dalla rugosità del terreno e dalla topografia: 


U(z) = cr(2)4(2)Ure; (11.3) 


dove cr (z) è il coefficiente di rugosità che tiene conto della variazione della velocità 
del vento con la quota per effetto della rugosità del terreno; c;(z) è il coefficiente di 
topografia che tiene conto della variazione della velocità eolica dovuta ad eventuali 
ostacoli naturali (e.g., colline, scarpate). 

Il coefficiente di topografia è generalmente assunto pari ad 1, mentre l’influenza 
del coefficiente di rugosità sulla velocità eolica è dato dalla Fig. 11.1. 
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15 20 25 30 35 40 45 


Um (m/s) 


Figura 11.1: Profili di vento per differenti rugosità secondo  EC1 


11.4.1 Coefficiente di rugosità 


Il coefficiente di rugosità è definito nel seguente modo: 


cr(2) = krln (=) ZS Bye (11.4) 
20 


cr(z) = Cr(Zmin) Z < Zmin (11.5) 
con kr fattore di terreno, zo lunghezza di rugosita del sito d’interesse, Zmin quota 


minima sopra il terreno. Questi coefficienti dipendono dalla categoria del terreno e 
sono riportati in Tabella 11.1. 


11.5 Coefficiente di raffica 


Nel Cap. 10 abbiamo definito il coefficiente di raffica nel seguente modo: 


plz) = 1 + 2glu(z)V ko + kr (11.6) 


dove g è il fattore di picco, I,,(z) è l’intensità di turbolenza alla generica quota z, 
mentre ky e kr sono due costanti relative rispettivamente alla risposta quasi-statica e 
in risonanza, e che tengono in conto anche la correlazione spaziale della turbolenza 
attraverso l’ammettenza aerodinamica x? (w). 

Moltiplicando e dividendo la 11.6 per la quantità 1 + 29,,(£), si ottiene: 


1+ 2gIu(z) Vki + kr 


(11.7) 
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Tabella 11.1: Categorie del terreno e relativi parametri 


Categoria del terreno kr zo (m) 2min(m) 
I Mare aperto, laghi con minimo 5 km di 0.17 | 0.01 2 
estensione sopravento e campagna senza 
ostacoli 
II Terreno agricolo con recinzioni, piccole fattorie | 0.19 | 0.05 4 
isolate, case o alberi 
III Aree suburbane o industriali o foreste 0.22 | 0.3 8 
permanenti 
IV Area urbana in cui almeno il 15% della 0.24 1 16 
superficie sia coperta da edifici la cui altezza 
media superi i 15m 


Inoltre esprimendo l’intensità di turbolenza come: 


MO Dale) ue 
si ottiene: 
alie, kr 1 + 2gl,(z)V kp + kr 
y(z) = |! 29a | 1+ YLA (11.9) 


11.6 Forzante eolica 


Tenendo presente le relazioni appena ottenute, la forzante eolica indotta dal vento è 
espressa nel seguente modo: 


1 1 
F(2) = 5 PU" (z)elz)er Ares = SPUfesG (A) ple) ep Ares (11.10) 


dove y assume il significato di coefficiente di amplificazione della pressione me- 
dia, cf è un coefficiente aerodinamico che varia in base alla struttura considerata, e 
Aref è un’area di riferimento. La relazione 11.10 può essere riscritta considerando la 
definizione del coefficiente di raffica data nella 11.9. In particolare, si ottiene: 


_1 2 20/2 kr 1+2g1L,(2)vko + kr 
F(z) = 5 PU rep r (zen (2) |! EEO FALA cp Aref 
(11.11) 


che puo essere riscritta nel seguente modo: 


F(z) = dref Ce(Z)Ca(z)cp Aref (11.12) 
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dove: 
ce(2) = elez) |1+ i (11.13) 
" cr(2)c(2) 
è il coefficiente di esposizione, mentre 
1 + g2Iu(z) Vko +k 
spe EE (11.14) 


1+ 921,(z) 


è il coefficiente dinamico. 


11.7 Azioni eoliche agenti sulle strutture 


La relazione 11.12 costituisce la formula generale con cui 1’ ECI definisce le azioni 
di calcolo agenti sulle strutture. In particolare, l' EC1 distingue tra pressioni (forze 


per unità di superficie) e forze globali. 


11.7.1 Pressioni 


Se le superfici strutturali sono sufficientemente rigide da poter trascurare il contributo 
delle loro vibrazioni risonanti, le pressioni interne ed esterne possono essere calcolate 


come è di seguito riportato. 


Pressione 
interno 
negativa 


` Pressione 
interno 
positiva 

WES 


ra 


= e 3 
SS SSA 
SS aasq»qww 


c) 


RX 


% 
Z 
Hd 
H 


Figura 11.2: Convenzione dei segni per le pressioni esterne ed interne 
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Pressioni esterne ed interne 

Le pressioni sono date dalle seguenti relazioni: 
We = re fe (Ze) Cpe (11.15) 
Wi = Gre fCi( Zi) Cpi (11.16) 


ripettivamente per le pressioni esterne ed interne. Le quote Ze e z; sono le quote di ri- 
ferimento per la struttura in esame, mentre Cpe € Cp; sono rispettivamente i coefficienti 
di pressione esterna ed interna e dipendono dalla forma della struttura. 


b 


Figura 11.3: Forzante applicata con un’eccentricità e = a 


Forze globali 


Le forzanti globali agenti sulla struttura, o su parte di essa, può essere calcolata in 
base alla relazione 11.12. In particolare, ?EC1 definisce tali forzanti in funzione di 
una quota di riferimento che dipende dal tipo di struttura in esame, z = Ze: 


Fy = dref Ce(Ze)CdCf Ares (11.17) 


con il consueto significato dei simboli. Alcuni valori del coefficiente dinamico sono 
riportati nelle Fig. 11.4, 11.5 e 11.6. 

Gli effetti torsionali dovuti ad un vento inclinato oppure alla non completa cor- 
relazione del vento turbolento, possono essere tenuti in conto applicando la forza Fw 
con un’eccentricita pari ad un decimo della dimensione della struttura normale alla 
direzione del vento, e = }; (Fig. 11.3). 

Inoltre L’ECI1 fornisce anche la forzante d’attrito dovuta al vento radente, valida 
soprattutto per quelle strutture che presentano grandi aree investite dal vento. Tale 
forza è espressa dalla seguente relazione: 


Ffr = drefCelZe)CfrAref (11.18) 


dove cr, è il coefficiente d’attrito. In questo caso si può notare che non rientra il 
coefficiente dinamico nella definizione della forzante. 
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Legenda 
1 Vento 


x Larghezza binm 
Y Altezza hinm 


200 


100 


Non di 
interesse pralico 


Figura 11.4: cg per strutture in calcestruzzo e muratura 


11.8 Procedura per il calcolo della risposta along-wind 


L’ECI prevede due tipologie di analisi per la valutazione della risposta along-wind 
delle strutture: una procedura semplificata ed una dettagliata. 

Le due procedure non sono applicabili per ponti a travata continua, ponti ad ar- 
co, ponti sospesi, per strutture snelle che non presentino un comportamento elastico 
lineare oppure che non abbiano il primo modo di vibrare nella direzione del vento 
sufficientemente disaccoppiato dai modi superiori o che non siano riconducibili agli 
schemi semplificati (verticale, orizzontale e puntiforme) in Fig. 11.7. 


11.8.1 Procedura semplificata 


La procedura semplificata fa riferimento alle strutture non suscettibili all’ amplifica- 
zione dinamica. In particolare si applica nel caso di edifici e ciminiere con altezza 
minore di 200 metri, per viadotti stradali e ferroviari la cui campata massima non 
superi i 200 metri. In ogni caso cg < 1.2. 

Per procedere alla valutazione della forzante globale agente sulla struttura si 
considera un fattore di picco pari a 3.5 (g = 3.5), pertanto avremo: 


Fw = drefCe(Zeg)cacfAref (11.19) 
dove 


Tkr 


— 11.20 
Cr(Zeq) Ct(Zeq) ( ) 


Ce(Zeq) = Cr (ei) Ge (Zeg) Lee 
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Legenda 
1 Vento 
X Larghezza binm 
Y Altezzahinm 


200 


100} 
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interesse pratico 


Figura 11.5: cg per strutture in acciaio 


con Zeg è l’altezza equivalente che dipende dal tipo di struttura in esame (Fig. 11.7), 
mentre il cg è dato dalle tabelle (e.g., Fig. 11.4, 11.5 e 11.6). 

L'andamento del coefficiente di esposizione al variare della quota è riportato in 
Fig. 11.8 


11.8.2 Procedura dettagliata 


La procedura dettagliata si applica a quelle strutture per le quali la componente dina- 
mica della risposta non è trascurabile. In particolare si fa riferimento a questo tipo di 
analisi ogni qualvolta cg > 1.2. Tuttavia, è consigliata anche quando 1 < cg < 1.2. 

Nell ambito della procedura dettagliata il coefficiente dinamico viene definito 
come segue: 


=: 1 + 29Iu(zeg) V Q3 + R2 (11.21) 
1+ 2gIu(Zeq) i 


Cd 


dove Zeq è l'altezza equivalente definita nella Fig. 11.7, Q? = kp e R2 = kr secondo 
la simbologia utilizzata nel Cap. 10. 
L'intensità di turbolenza viene valutata secondo la seguente espressione: 


1 


11.22 
cilea) In (Ze9/20) oa 


Tu(Zeq) = 
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Legenda 
1 Vento 
X Larghezza binm 
Y Altezzahinm 


200, 


b 100} 
IR | 
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È 


Figura 11.6: cg per strutture miste acciaio-calcestruzzo 


mentre il fattore di picco è dato dalla seguente formula: 


0.6 
= y 2ln(vT) + ———: (11.23) 
4 7) 2In(vT) 
dove T = 600 secondi è il tempo su cui è mediata la velocità di riferimento del vento, 
mentre v è definito come segue: 


202 2 p2 
Vo Qo aa Mali 
QR? 
dove nı x è la frequenza di vibrazione della struttura nella direzione del vento, os- 


sia lungo x, e vo è la frequenza attesa delle azioni di raffica su strutture rigide. 
L'andamento del fattore di picco è riportato in Fig. 11.9. 


(11.24) 


Componente quasi-statica della risposta 


La risposta quasi-statica viene quantificata attraverso la seguente formula: 
9 1 
Qi = OES (11.25) 


b+h 
1+0.9 ky 


dove L;(Zeq) è la scala integrale di turbolenza alla quota di riferimento (Fig. 11.11), 
b e h sono rispettivamente la larghezza e l’altezza della struttura in esame. In Fig. 
11.10 è riportato l'andamento di Q3. 
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Strutture verticali 
(per esempio: edifici, ciminiere, torri, ecc.) 
Zogu = 98 A+ Zamin 


Strutture orizzontali 
(per esempio: ponti, cavi, ecc.) 
h=b.d=beb/h, 205 


Zequ= hy + h!2> Zanin 


Strutture puntiformi 
{per esempio: lampioni, serbatoi sospesi, ecc.) 
b/h, 0.5; h/h, -0.5 


Zogu =" +h! 22 Zmin 


5 
Ce (2) 


Figura 11.8: ce al variare della quota, nelle quattro categoria di terreno. c = 1 
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0,ol___| 
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Figura 11.10: Andamento della componente quasi-statica della risposta 
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Linm 
zinm 
£ = 0,13 - categoria 1 
e = 0,26 - categoria 2 
e = 0,37 - categoria 3 
e = 0,46 - categoria 4 


IS Io (ha < x 


ee, RDZ 


CO nn 
EHH 
X 


Figura 11.11: Scala integrale di turbolenza L;(z), per diverse categorie di terreno 


III 
TE 


Figura 11.12: Densità spettrale adimensionale della risposta, Ry 
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Componente in risonanza della risposta 


La risposta in risonanza è data dalla seguente relazione: 


2 
R= 33 EN Rn Re (11.26) 
dove 6 = 27vs è il decremento logaritmico dello smorzamento, vs è il rappor- 


to di smorzamento strutturale. Ry è la densità spettrale adimensionale del vento 
turbolento (Fig. 11.12): 


Ni x Owl te) 6.8 Nz 
Ry = AE = (11.27) 
o2 (1 + 10.2N,) 
con 
N, = Mnbilzea) (11.28) 


Um (Zeq) 


L’ammettenza aerodinamica è definita lungo la larghezza e l’altezza della strut- 
tura, rispettivamente Ry ed Rp: 


1 1 


= -2 
oa n) per n >0 (11.29) 
R,=1 per 7 =0 (11.30) 
con 
4.6Ny2h 
Ryn = Ri ponendo n= 1, (11.31) 
Li(Zeq) 
4.6N; zb 
Ry= Ri ponendo n= ——— (11.32) 
Li(Zeq) 


Nel caso in cui le forme modali nella direzione del vento cambiano segno (pre- 
sentano cioè dei punti nodali) l’ammettenza aerodinamica dovrà essere valutata con 
procedure più dettagliate. 
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TROIA 
tit] | {tly 


Figura 11.13: Ammettenza aerodinamica, R; (l=h,b) 


a) 


> 


b) e 


Figura 11.14: a) Edificio in acciaio: h = 10 + 200m, b = 20m, d = 10m; b) Torre 


= h/d. 


in calcestruzzo con snellezza À 
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11.9 Esempi 


Consideriamo alcuni esempi che mettono a confronto la procedura semplificata e 
dettagliata. 


11.9.1 Edificio in acciaio 


La struttura è situata in un’area urbana (categoria 2), la velocità di riferimento del 
vento è U,ef = 27m/s. Il coefficiente di forza è cf = 1.3. Secondo l’ECI, la 
procedura semplificata la si può applicare per un’altezza che non superi i 50m (Fig. 
11.15). 


11.9.2 Edificio basso in acciaio 


La struttura è situata in un’area urbana (categoria 2), la velocità di riferimento del 
vento è Ure = 27m/s. Il coefficiente di forza è cf = 1.3. In queste ipotesi, si 
ottengono i seguenti valori dell’azione eolica: 


e Procedura semplificata: Fu /Aw = 1.40K N/m? 
e Procedura dettagliata: F,,/A, = 1.19K N/m? 


La forzante eolica calcolata con la procedura semplificata è circa il 18% più grande 
di quella calcolata con la procedura dettagliata. Ciò è dovuto al fatto che la procedura 
semplificata non tiene conto degli effetti dovuti alla dimensione della struttura. La 
dimensione della struttura è tale che l’ammettenza aerodinamica diminuisce. 


11.9.3 Torre in calcestruzzo poco snella (À < 12) 


La struttura è situata in un’area urbana (categoria 2), la velocità di riferimento del 
vento è Ure fp = 27m/s. 

Per i criteri sulle torri, questa struttura può essere calcolata con la procedura 
semplificata. Sia la procedura semplificata sia quella dettagliata forniscono lo stesso 
valore della forzante eolica visto che gli effetti dovuti alla dimensione della struttura, 
e la sua risposta risonante sono trascurabili: 


e Foy/Awy = 1.47 KN/m? 


11.9.4 Torre in calcestruzzo molto snella (A > 12) 


La struttura è situata in un’area urbana (categoria 2), la velocità di riferimento del 
vento è Ure = 27m/s. 

Per i criteri sulle torri questa struttura deve essere calcolata con la procedura 
dettagliata. 
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Figura 11.15: Confronto tra la procedura semplificata e quella dettagliata 


e Fw/Aw = 2.23 KN/m? 


Con la procedura semplificata si otterrebbe una forza per unità di superficie pari al- 
la metà di quella ottenuta con la procedura dettagliata. La snellezza della struttura 
induce una rilevante risposta in risonanza. 


Risposta strutturale across-wind 


In questo capitolo viene discussa la risposta strutturale nella direzione ortogonale a 
quella del vento (across-wind), dovuta principalmente al fenomeno del distacco dei 
vortici. 

Prima di discutere il problema strutturale, conviene introdurre alcune nozioni 
relative al fenomeno del distacco dei vortici e alla formazione della scia. 


12.1 Fenomeno del distacco dei vortici 


Un corpo immerso in una corrente fluida è delimitato da uno strato limite superficiale 
che si separa alternativamente dalle regioni laterali dando vita ad una scia di vortici 
che defluisce a valle del corpo. La natura della scia dipende sensibilmente dal numero 
di Reynolds: 


par (12.1) 


dove v è la viscosità cinematica del fluido, U è la velocità media del flusso inci- 
dente, e D è la dimensione caratteristica del corpo nella direzione normale a quella 
del flusso. Nella Fig. 12.1 è riportata schematicamente l’effetto del numero di Rey- 
nolds sul processo di formazione della scia. Per valori bassi di Re, circa Re < 5, 
le linee di flusso seguono il contorno della sezione (flusso non separato). Nell inter- 
vallo 5 < Re < 40, il flusso si separa nella regione a valle del cilindro formando 
una coppia di vortici simmetrici. Se il il numero di Reynolds viene incrementato 
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Figura 12.1: Caratteristiche della scia di un cilindro circolare al variare di Re 


ulteriormente, 40 < Re < 150, la scia diviene instabile e i vortici si staccano al- 
ternativamente dalla regioni laterali del cilindro formando una scia di vortici di tipo 
laminare (regime laminare) Nell’intervallo 150 < Re < 300 la scia comincia a diven- 
tare turbolenta anche se lo strato limite superficiale sul cilindro è ancora laminare. I 
numeri di Reynolds compresi tra 300 < Re < 3- 10° caratterizzano una scia comple- 
tamente turbolenta (regime subcritico). Il regime critico si manifesta nell’intervallo 
3-105 < Re < 3.5 - 10° ed è caratterizzato da una scia turbolenta e disorganizzata a 
tal punto che non è più distinguibile il distacco dei vortici. Infine, per Re > 3.5 - 10° 
si ha il regime supercritico dove riappare nuovamente il distacco dei vortici all’inter- 
no di una scia turbolenta. 

I vortici si staccano dalle regioni laterali del corpo con una frequenza, fus, che 
segue la nota legge del numero di Strouhal: 


U 
fos = Si (12.2) 


dove S; è il numero di Strouhal che dipende da diversi parametri, quali ad esempio la 
forma della sezione ed il numero di Reynolds. Si può facilmente notare che la legge 
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Figura 12.2: Legge del numero di Strouhal 


del numero di Strouhal è una legge lineare (Fig. 12.2). 


12.2 Sincronizzazione del distacco dei vortici 


Quando la frequenza del distacco dei vortici uguaglia la frequenza di oscillazione 
trasversale del cilindro elasticamente sospeso o ad elasticità diffusa, si instaura una 
complessa interazione fluido-struttura che si traduce in notevoli ampiezze di oscil- 
lazione, pari addirittura al diametro del cilindro per cilindri poco smorzati. La sin- 
cronizzazione è un fenomeno aeroelastico molto complesso di cui, in questa sede, si 
ricordano le caratteristiche principali d’interesse strutturale. 

Considerando la 12.2, e ponendo fy; = fo, dove fo è la frequenza di oscillazione 
trasversale del cilindro, si ottiene la velocità critica, ossia la velocità per la quale si 
ha la piena sincronizzazione del distacco dei vortici: 

Ver = fD (12.3) 

St 
Tuttavia, il fenomeno si estende su un’intervallo di velocità, intorno alla velocità cri- 
tica, all’interno del quale la legge del numero di Strouhal risulta essere violata (Fig. 
12.3 e 12.4). La risposta diminuisce all aumentare dello smorzamento strutturale e si 
assesta su valori della velocità critica decrescenti (Fig. 12.3). Inoltre, all aumentare 
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dello smorzamento, l’intervallo di sincronizzazione diminuisce ed arretra leggermen- 
te (Fig. 12.4). 
La massima risposta risulta dipendere fortemente dal parametro massa smorza- 
mento: 
MVs 


K, = SD? (12.4) 


dove m è la massa per unità di lunghezza, p è la densità dell’aria, vs è il rapporto 
di smorzamento meccanico, e D è la dimensione trasversale della sezione del corpo. 
Generalmente si preferisce fare riferimento al numero di Scruton: 


Atm s 


Se = 41K, 2D? 


(12.5) 


ed esistono alcune relazioni matematiche che esprimono la massima ampiezza di 
oscillazione in funzione del numero di Scruton: si riportano la relazione di Griffin 
and Ramberg [8]: 


1.2 
(4) = 2a (12.6) 
d/max [1 + 0.43 (2782 Se)” 
e quella di Sarpkaya [14]: 
CE Sra (12.7) 


dove Sg è il numero di Skop-Griffin ed è un ulteriore parametro di massa smorza- 
mento. In Fig. 12.5 sono riportati alcuni risultati sperimentali comparati con le curve 
teoriche nelle relazioni 12.6 e 12.7. 
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Figura 12.3: Ampiezza adimensionale di oscillazione , $/D, al variare dello smor- 
zamento vs e del numero di Scruton Se: o vs = 0.243%, 5. = 3.1632; + 
v = 0.388%, Se = 5.05; * v = 0.459%, Se = 5.98; X v = 0.502%, Se = 6.54; 
Sı = 0.2074. Pastò [12] 
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Figura 12.4: Intervallo di sincronizzazione al variare dello smorzamento vs e del 
numero di Scruton Se: o vs = 0.243%, Se = 3.1632; + v = 0.388%, Se = 5.05; * 
v = 0.459%, S = 5.98; X v = 0.502%, Se = 6.54; Sı = 0.2074. Pastò [12] 
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Figura 12.5: Smooth flow: (3) max in funzione di Se: — Griffin and Ramberg [8] (Eq. 
12.6); - — - Sarpkaya [14] (Eq. 12.7); O Goswami et al. [7]; X Feng [6]; +, 4 Pastò 
[12] 


Un’ altra caratteristica del fenomeno della sincronizzazione è l’aumento della cor- 
relazione della forzante di Lift lungo il cilindro (Fig. 12.6). Per ampiezze di oscilla- 
zioni molto elevate la forza di Lift risulta quasi completamente correlata, mentre la 
risposta diviene isteretica, ossia raggiungendo la sincronizzazione per velocità cre- 
scenti o decrescenti si ottengono due curve differenti (Fig. 12.8). 

Da un punto di vista strutturale è importante ricordare che sebbene la riposta 
possa essere notevolmente amplificata dall’interazione fluido-struttura, le oscillazioni 
strutturali non divergono ossia presentano un ciclo limite per una data velocità com- 
presa nell’intervallo di sincronizzazione (Fig. 12.7). Ciò permette di poter mitigare 
la risposta agendo, ad esempio, sullo smorzamento strutturale. 


12.3 Risposta across-wind delle strutture 


In questa sezione viene proposto il modello di Vickery and Basu [18], valido per 
strutture snelle verticali a sezione circolare (e.g., ciminiere, torri di telecomunica- 
zioni etc.), che può essere esteso a qualsiasi tipo si sezione, previa conoscenza dei 
parametri necessari per l’analisi. 

I parametri che ricorrono nel modello di Vickery and Basu [18] hanno origine da 
dati sperimentali su cilindri circolari elasticamente sospesi (Fig. 12.9). 

L'equazione di moto associata al generico modo di vibrare è la seguente: 
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Figura 12.6: Correlazione della forzante di Lift lungo il cilindro in condizione di 
sincronizzazione 


0.75 


-0.75 


Figura 12.7: Ciclo Limite della risposta in condizioni di sincronizzazione del distacco 
dei vortici 


So VOF(, dz _ 
Jo ¥? (@)m(2)de 
fe b(2 [Fus(z,t) + Folz, t) + Fae(2,1)] dz 
Jo ¥?(2)m(2)dz 
dove wo = 27 fo è la frequenza angolare della struttura, associata al generico modo di 
vibrare. L'equazione 12.8 è ottenuta mediante una decomposizione modale del pro- 


blema dinamico della struttura ($ 9.5). La risposta strutturale across-wind associata 
a tale modo di vibrare sarà dunque data da: 


Y(2,t) = (z)y(t) (12.9) 


Y + 2Qewoy + wy = 


(12.8) 


Vediamo in dettaglio le vari componenti della forzante F(z, t), supposte scorre- 
late tra loro: 
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Figura 12.8: Isteresi dell’ampiezza adimensionale di oscillazione, $/D. 4 velocità 
crescente; + velocità decrescente. Pastò [12] 


12.3.1 Forzante dovuta al distacco dei vortici 
Il distacco dei vortici induce una forzante normale alla direzione del vento (forza di 
Lift) pseudo-armonica: 


1 
Fys(z,t) = 5 Pe D(2)Cr(2,1) (12.10) 


dove D(z) è la dimensione normale al vento (nella presente trattazione è il diametro 
della struttura). La quantità C7,(z,t) è il coefficiente di Lift per unità di lunghezza, 
e cè, è la sua varianza. La densità spettrale del coefficiente di Lift è generalmente 


espressa come segue: 


(12.11) 


ossia come una sorta di campana di Gauss (Fig. 12.10) la cui larghezza è determinata 
dal parametro B che dipende dall’intensità di turbolenza: 


B? = BẸ + 21? (12.12) 


dove Bo è il valore di B in flusso laminare, generalmente compreso tra 0.05 e 0.1. 
Tenendo presente la 12.8, la correlazione della forzante indotta dal distacco dei 
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Figura 12.9: Modello di un cilindro circolare elasticamente sospeso 


vortici è data dalla seguente relazione: 


Rie: =f ia G out) D(z1)D(22)E [Cr (21,t)Cr(22,t + 7)] x 


x p(z y (z2 \dzydzq = 


=f S (ot) D(21)D(z2) Ro, cy, (21, 22,7) (21) p(z2)dz1dze 
(12.13) 


Operando la traformata di Fourier della 12.13 si ottiene lo spettro della forzante 
di Lift: 


h h 2 
Sant = ff (Sout) Mbit 
X w( 21 \(z2)dz1dz2 = 


al: [ È ait.) D(21)D(22)\/ Serer (21, f) Serer (22, f)x 


x Cohys(21, 22, f)Y(a1)P(22)dzidzq (12.14) 


dove la funzione di coerenza è indipendente dalla frequenza 


lz1 -272| 
Cohys(21, 22, f) = cos( 11-22), al ( 3 ) | (12.15) 
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Figura 12.10: Spettro di Cz al variare della turbolenza 


12.3.2 Forzante dovuta alla turbolenza laterale 


Nell'ambito della teoria quasi-stazionaria linearizzata, la turbolenza trasversale as- 
sociata alla velocità v (lungo y) ha l’effetto di cambiare l’angolo d’attacco a (a ~ 
v/uy) della velocità media uy. Ne consegue una forzante del tipo: 


1 
F,(2,t) = 5PP(2)Cn(z) emu (12.16) 


dove Cp è il coefficiente di trascinamento per unità di lunghezza. Lo spettro della 
forzante indotta dalla turbolenza trasversale è dato da: 


h rh 2 
Snath = ff (Zoru) Co(a)Cr(a)D(a)Dla)x 
x V Solz, PS, f)x2(f)Cohy(z1, 225 Pelz (49) dz dzo (12.17) 


dove x?(f) è l’ammettenza della componente turbolenta trasversale, generalmen- 
te espressa a partire dall’ammettenza della turbolenza longitudinale u (Vickery and 
Basu [18]). 


12.3.3 Forzanti aeroelastiche 


L'interazione aeroelstica tra il flusso e la struttura si traduce in un’amplificazione 
della risposta che talvolta può essere anche notevole e insostenibile per la struttura. 
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Figura 12.11: Massimo valore di Kao al variare di Re 


Per riprodurre numericamente tale amplificazione si ricorre ad una forzante viscosa 
che induce uno smorzamento aerodinamico: minore è lo smorzamento della struttu- 
ra maggiore è la rispsosta strutturale (8.6). Tuttavia, si ricorda che la risposta della 
struttura è comunque limitata (presenta un ciclo limite) pertanto lo smorzamento to- 
tale dovrà essere sempre maggiore di zero. Sperimentalmente si osserva che, per 
ampiezze di oscillazioni molto grandi, forzanti viscose lineari non garantiscono un 
ciclo limite, pertanto si ricorre ad uno smorzamento non-lineare. 
In particolare, la forzante aeroelastica di natura viscosa è definita come segue: 


fae,a(2,t) = ka C - GY?) (12.18) 


dove ka e G sono delle costanti, e Y (z, t) = W(2)y(t). 
Tuttavia, si preferisce utilizzare comunque una forzante viscosa linearizzata, per- 
tanto avremo: 


wee 


2 
fae,a(2,1) = 47 fopD?(z)Kao 1-( at] w(z) y(t) (12.19) 


dove il termine (0y/(AD(z))? garantisce la presenza di un ciclo limite per oscilla- 
zioni molto elevate. Il parametro Kao è definito come segue: 


__ 6 {al} 
IPSE 


Kao (12.20) 
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Figura 12.12: Kao/Ka 


al variare di py /W,,. 


Omax 


dove 8 = uu/Wu., 
La forzante modale (secondo la 12.8) sara: 


h x 2 
Faea = Y | 4r fopD(z)Kao | — es Lou) | w?(z)dz (12.21) 


Dal momento che la Fae è proporzionale a y allora è possibile definire il coefficiente 
di smorzamento modale Cae: 


h T D(z) 
[ | Dref | Kao De; ao ( dd aa) 


avendo moltiplicato e diviso la forzante fae q per un diametro di riferimento, D2, f 
Si introduce la massa equivalente me: 


Cae = 4r fopD?o5 


h 2 
z)m(z 
Me = h A ude (12.23) 
Jo p? (z) 
e possiamo dunque scrivere (sempre con riferimento alla 12.8): 
ea E (12.24) 
Me 


da cui si ottiene infine il rapporto di smorzamento aeroelastico, Vae 


pD? 2 
ie ref ©» ( oy ) | (12.25) 
Me Dref 
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essendo 
fe Kao [BO] sota 
ci = dai (12.26) 
INEZIONE 
h 4 
selle ek (12.27) 
So (dz 


Inoltre, sperimentalmente si osserva che entro l’intervallo si sincronizzazione, il 
flusso e la struttura si sincronizzano non solo per valori di fys/ fp = 1, ma anche per 
valori che si discostano dall’unità. Ciò giustifica l’introduzione di una forzante che 
induce una rigidezza aerodinamica: 


fae(2,1) = 87° f20D° (z)Hat(z)y(t) (12.28) 


dove l’andamento di H, è riportato in Fig. 12.13. Sebbene per strutture immerse in 
un campo di vento si può assumere f,s/fo = 1 (che equivale a porre Ha = 0) la 
componente in fase della forzante aeroelastica, fae,x, gioca un ruolo fondamentale 
da un punto di vista fisico (cfr., e.g., Pastò [12]) per l’innesco del fenomeno della 
sincronizzazione. In questa sede si preferisce non appesantire la trattazione della 
sincronizzazione del distacco dei vortici, e assumeremo trascurabile la componente 
in fase della forzante aeroelastica, pertanto Fue(2,t) = Fae,d- 


Figura 12.13: Haal variare di uu/4u,, 
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12.3.4 Valutazione della risposta strutturale sincronizzata 


Come più volte menzionato nei capitoli precedenti, la risposta del sistema (decompo- 
sto) ad un grado di libertà (Eq. 12.8) può essere calcolata nel dominio delle frequenze 
secondo la formula: 


Syy(f) = H(A) Sreo) (12.29) 


dove l’ammettenza strutturale è definita come al solito nel seguente modo: 


2 1 1 
H(p) = Tor: I mE (12.30) 
(1 B) Pal (4) | 
con 
h 
M= I w?(z)m(z)dz (12.31) 
0 
Y= Vs + Vae (12.32) 
Lo spettro della forzante invece é dato dalle relazioni 12.14 e 12.17: 
Srr(f) = SF, Fos (f) + SF FCP) (12.33) 


nell’ ipotesi che la forzante indotta dal distacco dei vortici e quella dovuta alla turbo- 
lenza laterale siano scorrelate. 

Siamo alla ricerca della varianza della risposta che nel caso di spettri unilateri 
viene calcolata come segue: 


oo +00 
a= fl SuD = IMPE iosa 


In condizioni di sincronizzazione avremo fy; = fo, pertanto la varianza della risposta 
potrà essere scritta nel seguente modo: 


+00 
o2 = Spr(fo) | H(A) Raf (12.35) 


da cui si ricava: 


ES SFos Fos (f) + Spr (f) 
” (4m fo)? M? (Vs + Vae) 


(12.36) 


La massima risposta sarà data da: 


Ye = Hy + kpoy = kpOy (12.37) 
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Figura 12.14: o,/D al variare di Ks 


dove uy = 0 dal momento che le fluttazioni indotte dal distacco dei vortici, sia in 
regime stazionario sia in condizioni di sincronizzazione, hanno media nulla. Con 
kp si indica il fattore di picco che può essere calcolato con le consuete formule. La 
risposta massima modale sarà data da: 


Ye(z) = W(z)Ye (12.38) 


L'andamento della risposta adimensionalizzata in funzione del parametro Ks (de- 
finito nell’ Eq. 12.4) è riportato in Fig. 12.14, dove si vede che le ampiezze in con- 
dizione di sincronizzazione (Locl-In) sono decisamente maggiori rispetto a quelle in 
regime forzato dove l’interazione aeroelastica svanisce. 


Risposta strutturale across-wind secondo 
L’Eurocodice 1 


La risposta in condizioni di sincronizzazione contemplata dall’Eurocodice-1 [5] si 
basa sul modello di Ruscheweyh [13] che viene qui riproposto. 


13.1 Modello di Ruscheweyh 


L'equazione di moto associata al generico modo di vibrare è data dalla seguente 
relazione: 


fev F(z tide 
So @m(2)dz 


dove wo = 27f, è la frequenza angolare della struttura, associata al generico modo 
di vibrare. L’Eq. 13.1 è ottenuta mediante una decomposizione modale del problema 
dinamico della struttura ($ 9.5). La risposta strutturale across-wind associata a tale 
modo di vibrare sarà dunque data da: 


Y(z,t) = Y(z)y(t) (13.2) 


+ Qvewoy + wi y (13.1) 


Ruscheweyh introduce la seguente forzante dovuta al distacco dei vortici in con- 
dizioni di sincronizzazione: 


f(z, t) = a(2)D(2)c;(2)cos(27 fot +9(2)7) (13.3) 
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dove q e D sono rispettivamente la pressione dinamica ed il diametro della struttura 
alla generica quota z, e c;(z) è un coefficiente aerodinamico. Nel caso in cui il modo 
di vibrare non presenti punti nodali, possiamo porre y(z) = 0, e la forza modale sarà 
data da: 


h 
F(z,t)= cos(2n fot) | yp(z)a(z)D(z)cf(z)dz (13.4) 


In condizione di sincronizzazione, fvs = fo, si ottiene la seguente risposta massima 
(ricordando che uy = 0): 


Ymax = WmaxYmax Gori WmaxKpOy (13.5) 


dove kp è il fattore di picco e gy è data da: 


Oy = [Tua PACF — fo af w(z D(2)c;(2)dz (13.6) 


ottenuta tenendo presente che il secondo integrale è una costante. La funzione ô( f — 
fo) è il 6 di Dirac, e corrisponde allo spettro di cos(27 fot). La funzione sotto radice, 
E |H(f)|?5(f — fo)df, è definita solo per f = fo, pertanto otteniamo: 


1 1 Sè (a) 


oy = (OTF) Me Ws fe (2) (13.7) 
ricordando che: 
1 1 
|H(A)I? ia 
M2w4 2\ 2 2 
[( L) +a? (£) | 
h 
vo de ma) 2) (13.9) 
M 
Messe (13.10) 
SE P(2) 


Definita con z* la quota in cui avviene la sincronizzazione possiamo scrivere la 
massima risposta in questo modo: 


fou DES, Da c;(2)dz 


Ymax 1 E 
(27 fo)? Me 2Vs 2 JÈ W2(z 


D(z *) = Umax 


(13.11) 


dove la pressione dinamica g(z*) = $pU?(z*)D(z*). Inoltre, ricordando la defini- 
zione del numero di Scruton: 
ANV Me 


e= pD? 


(13.12) 
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e del numero di Strouhal 


fos D 


Sy = U 


(13.13) 


considerando una struttura a sezione costante D(z) = D, e un profilo di velocita 
uniforme, q(z) = q, si ottiene: 


Ymar =y) „ht z)dz 
D(2*) 5 2 Kp = 


(13.14) 


A questo punto si introduce la lunghezza di correlazione, L* definita come il dop- 
pio dell’integrale della funzione di cross-correlazione della forzante. In particolare 
Ruscheweyh pone: 


h 
by | c;(2)v(2)dz = ciak f w(z)dz (13.15) 
0 L* 
dove ciat è la deviazione standard della forzante applicata. 


Successivamente, Ruscheweyh introduce l lunghezza effettiva di correlazione, 
Le, tale che: 


L* Le 
kp W(z)dz = / wz)dz (13.16) 


Infine la massima risposta in regime di sincronizzazione viene scritta nel seguente 
modo: 


Ymaz 1 1 
Dx) Mv eag g (13.17) 
dove 
h 
mae 22 ud (13.18) 
4T da yp? 
Le 
i o (13.19) 
Jo p(z) 
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Tabella 13.1: Valori del coefficiente aerodinamico cat 


Rapporto critico di velocità del vento Clat,0 
Veri st 
Vai < 0.83 Clat = Clat,0 
0.83 < feriti < 1.25 ciat = (3 — 2.442%) ciato 
Vi sa 
` Vink Clat = 
aoe i 
1,0 
Gat, 0 
0 


Il Il 
i LI [UT {Do DI 
a inte 
Rai al ae 
i IE AH 


0 
104 2 345 7 10° 2 345 7 10° 2 345 7 107 2 3 


Figura 13.1: Valori di cjg¢,. per cilindri circolari 
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Figura 13.2: Esempi di lunghezza effettiva di correlazione 
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Strut Forma modale 
Piz) Ky K 
0,13 
0,10 
vedere figura C.9 L/b 4 L/b 
n=1;m=1 el sen] -(1-2)| 0,11 
m 
X flew de 
analisi modale eri 
n=3 ~ 1 0,10 
m=3 
L fiw de 


intl 


Figura 13.3: Valori di Ky, 


Risposta aeroelastica delle strutture 


Le strutture aventi una sezione non circolare sono soggette ad una forzante fluida che 
cambia con l’orientamento del corpo rispetto al flusso. Se la struttura oscilla, la sua 
posizione cambia e la forzante fluida oscilla a sua volta. Se la forzante fluida tende 
ad aumentare le oscillazioni strutturali, la struttura diviene instabile e si possono ma- 
nifestare grandi ampiezze di oscillazione. 

In questo capitolo studieremo la risposta alle oscillazioni galoppanti e in regime 
di Flutter. 


14.1 Galloping traslazionale 


L'analisi delle oscillazioni galoppanti si basa sull’ipotesi secondo la quale la forzante 
fluida agente sulla struttura dipende soltanto dalla velocità eolica istantanea (teoria 
quasi-statica), pertanto il fenomeno può essere studiato in galleria del vento su mo- 
delli stazionari, misurando la forza fluida al variare dell’angolo d’attacco. 

È doveroso ricordare che la teoria quasi-statica può essere ritenuta valida soltanto 
se la frequenza di distacco dei vortici, fys, è molto più grande della frequenza propria 
di oscillazione della struttura, fo. La condizione fys >> fo è verificata se: 


U 
> 20 (14.1) 
foD 
dove FD è la velocità ridotta, U = puy è la velocità media del flusso, e D è la di- 


mensione trasversale della sezione strutturale. 
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Fp 


Figura 14.1: Modello ad un grado di libertà per le oscillazioni galoppanti 


Tuttavia, alcune strutture presentano oscillazioni galoppanti per valori della ve- 
locità ridotta compresi nell’intervallo 1 + 20 (1 < 5 < 20), intervallo in cui 
l’applicazione della teroria quasi-statica è molto discutibile. 

In Fig. 14.1 è schematizzato il modello di un sistema ad un solo grado di li- 
bertà che si muove con velocità y, e soggetto ad una forzante di Lift ed una di Drag 


espresse come segue: 
1 
F, = 5 PV; DCL (14.2) 
1 
Fp = 3PVr DCD (14.3) 


e dirette rispettivamente ortogonalmente e parallelamente alla direzione della velocità 
relativa V, = VU? + ġ2. I coefficienti Cz e Cp sono chiaramente i coefficienti di 
Lift e di Drag. 

Il moto del sistema cambia la sua posizione rispetto alla velocità del flusso che di 
per sé risulta inclinata di un angolo a rispetto alla direzione del flusso indisturbato: 


1/9 
= tan | 14.4 
a = tan (4) (14.4) 
L'equazione di moto del sistema sarà dunque: 
$ f 2 1 
Ü + Wewyy + wy = —Fy(a) (14.5) 
m 
dove: 


F,(0) = 5eV2DCy(a) (14.6) 
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CL 


1.2 
Stable 
0.8 ` 


Unstable 


-129/-89 12°a 


Figura 14.2: Esempio di curve stabili ed instabili per Cz 


Cy(a) = — [Cr cos(a) + Cp sin a] (14.7) 

Per piccoli angoli d’attacco (a — 0), linearizzando le quantità d’interesse, si ottiene: 
y 

Ri 14.8 

aX (14.8) 

VU (14.9) 


Cy(a) = Cy(0) + Cy,a(0) a = -C1 (0) — [CL a (0) + Cp(0)] a (14.10) 
con C(.),a = OC(.)/Oa. 
Possiamo quindi scrivere l equazione di moto nel seguente modo: 
: RARO, 1 2 
Ü + 2Vswy) + wyY = ImbV DC,(a) = 


= 5 pU?D {-0,0 — [CL a(0) + Cp(0)] si (14.11) 


ossia: 
pUD 
Y 


È ; 1 
Ü+2wy fu, Pr [Ct a(0) + coro} y+wiy = -zp PU” DC1(0) (14.12) 


Dalla 14.12 è facile constatare che il moto della struttura induce una forza fluida 
viscosa che può essere tenuta in conto definendo uno smorzamento totale dato dalla 
somma dello smorzamento strutturale e di quello aerodinamico: 


UD 
vi = vs + n. [Cr a(0) + Cp(0)] (14.13) 


da cui si evince subito che il sistema è stabile se: 


[CL,a(0) + Cp(0)] > 0 (14.14) 
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aC,/da" 
Section Smooth flow Turbulent flow’ Reynolds number 
BEL 3.0 3.5 10° 
rt 
[Jie 0. -0.7 10° 
+ go 
EZ -0.5 02 10° 
>j je D/2 
D 
LI -0.15 0. 10° 
>| e-0/4 
+. 
30 C] 
Ti n 1.3 12 66 000 
D 
4. 
a E 2.8 -2.0 33.000 
Lira 
4. 
Dili -10. - 2 000-20 000 
| > 
| SR PTE | 
ji -6.3 -6.3 >10 
D 
doi 
es] -6.3 -63 >10 
D 
Dio -0.1 0 66 000 
Tua -0.5 29 51 000 
0.66 si 75 000 


Figura 14.3: CL a sezioni. Per Cz aœ < 0 ovvero CL a + Cp > 0 l'equilibrio risulta 


stabile 
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— 


RÖ SIN y/Rd 


Figura 14.4: Modello di Galloping torsionale 


Infatti, la soluzione dell’ equazione di moto può essere scritta: 


1 2 
—+pU?2DC_,(0 
y= 2P E (0) t Aye "©! sin (wyy1- 2+ 8) (14.15) 


da cui si evince che se lo smorzamento totale, v, diviene negativo, per effetto di 
[CL a(0) + C'p(0)] < 0, la struttura diviene instabile in quanto la sua risposta diverge 
(si veda l’Eq. 14.15). In Fig. 14.2 sono riportate due curve di esempio per il Cz. 

Pertanto, ponendo 1 = 0, è possibile trovare la velocità eolica critica che induce 
oscillazioni galoppanti: 


8am fyvs (IC, 
Uer = y's y 14.16 
pD ( Oa ) i ) 
In Fig. 14.3 sono riportati i valori di Cya = — [Cr,a(0) + Cp(0)] per alcune 


sezioni (Blevins [1]). 


14.2 Divergenza e Galloping torsionale 


Nel Galloping torsionale, l’ angolo d’attacco cambia con la posizione angolare 0 e con 
la velocità angolare 6 = 00 /Ot. La velocità angolare, 6, induce un angolo d’attacco 
che varia lungo la sezione: una velocità d’attacco positiva induce un angolo d’attacco 
negativo e viceversa. 

Dato un generico punto sulla sezione, l’angolo d’attacco è espresso come segue: 


a = 0 — tan”! _Rôsiny (14.17) 
U — Ro cos y 


mentre la velocità relativa è data dalla seguente relazione: 


V? = (Rô sin 1). + (u — RÊ cos DI (14.18) 
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dove U = py è la velocità del flusso, e R = Resin y è il raggio caratteristico della 

sezione. Nel caso di un rettangolo, R è pari a metà larghezza del rettangolo stesso. 
In generale, R dipende dalla sezione considerata. 

Per piccoli angoli di attacco (a — 0) si ottengono i seguenti valori dell’angolo di 

attacco e della velocità relativa: 
RO 
a=t—- T (14.19) 
vV? = U? (14.20) 


La sezione sarà pertanto soggetta ad un momento torsionale, per unità di lunghez- 
za, dato dalla seguente relazione: 


Mg = PU? D?Cu (14.21) 
che può essere espansa in serie di Taylor fornendo la seguente forzante: 

Mp = 5eU?D? [Cx (0) + Cua(O)a+...] (14.22) 
dove: 

Cma(0)a = Cm,a(0) ( = 7) (14.23) 


L'equazione di moto del problema dinamico viene scritta nel seguente modo: 


3 , 1 1 RÒ 
Joð+2Jgvowgð +kgð = 3PU’D?’Cm(0)+5pU?”D’Cm,a(0) ( = =) (14.24) 


Conviene riscrivere l’equazione di moto nel seguente modo: 


Job + (2mo + 30URD'Cx10(0)) 64 (i n 300 D°Co(0)) 0 = 
= È pU?D°Cy(0) (14.25) 
da cui si possono notare due tipi di instabilità, una dinamica (Gallopping torsionale), 
ed una statica (Divergenza) 
14.2.1 Divergenza 


La Divergenza è un’instabilità statica perché si verifica quando la somma della ri- 
gidezza rotazionale e della rigidezza aerodinamica è nulla, ossia quando (cfr. Eq. 
14.25): 


1 
ko — 3PU?DC1,0(0) =0 (14.26) 
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che corrisponde alla seguente velocita eolica critica: 


2ko 
po 14.27 
cr pD? CMa (0) ( ) 
Section 7 ICu/da Reynolds number 
Eie —0.18 104-105 
2D 
le y 
10 —0.64 5 x 10°-5 x 10° 
F 
C + Jw -18. 2x 10-2 x 10° 
F 
| 5D p 


-26, 2x 102 104 


>10 


olf 


wie a fe 


Figura 14.5: Cva (0) per alcune sezioni 


14.2.2 Galloping torsionale 


Nel caso in cui la somma dello smorzamento strutturale e di quello aerodinamico è 
nulla, si ha il Galloping torsionale la cui velocità critica si calcola analogamente a 
quanto visto precedentemente per il Galloping traslazionale: 


1 
2Jovgwe + 3PURD?C1,0(0) =0 (14.28) 


da cui si ottiene la corrispondente velocita critica: 


AJovewe 1 


14.2 
pD?R Cuia(0) aon 


Ucr = 


In Fig. 14.5 sono riportati i valori di C'M,a(0) per alcune sezioni. 


14.3 Galloping di sistemi a due gradi di libertà 


Nella realtà, le strutture possono sia ruotare che traslare. Pertanto, in questa sezione 
viene estesa la trattazione precedentemente illustrata, valida quando i due modi di 


192 LEZIONI DI INGEGNERIA DEL VENTO 


Figura 14.6: Sezione generica a due gradi di liberta 


vibrare possono essere considerati indipendenti l’ uno dall’ altro. 

La rotazione e la traslazione di una generica sezione immersa in un campo di 
vento (si veda la Fig. 14.6), sono aerodinamicamente accoppiate attraverso l’angolo 
d’attacco che in questo caso si scrive: 


Y 
a=0 uta (14.30) 
Inoltre, la rotazione e la traslazione sono accoppiate inerzialmente quando il centro 
elastico della sezione non coincide con il centro di massa della sezione stessa. Il 
centro elastico è quel punto della sezione dove qualsiasi carico che vi è applicato non 
produce effetti torsionali. 

Con riferimento alla Fig. 14.6, l’equazioni di moto del sistema a due gradi di 


libertà si scrivono come segue: 


is 1 
mij + 2mvywyy + Sz0 + kyy = zPU’DO;(a) (14.31) 
1 1 
Jo + 2Javgwoð + Sri + ko = 5 pU?’ D?Cumla) (14.32) 
dove: 
m= i Lmdédn (14.33) 
A 
Jo= f (+1? )umdédy (14.34) 
A 


Sy = f AT (14.35) 
A 
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sono rispettivamente la massa per unità di lunghezza, il momento di inerzia polare 
ed il momento statico della sezione. Inoltre, um, è la densità di massa, (€, n) sono le 
coordinate di un generico punto sulla sezione. 

L'angolo d’attacco a accoppia aerodinamicamente i due gradi di libertà della 
sezione, 0 e y (Eq. 14.30), mentre il momento statico Sz li accoppia inerzialmente. 

L’equazioni di moto del sistema risultano pertanto dipendenti tra loro. Tuttavia, 
attraverso una decomposizione modale è possibile disaccoppiare il sistema che può 
essere così studiato come la somma degli effetti di due sistemi ad un solo grado di 
libertà. Tralasciando i passaggi intermedi, si arriva alla definizione della velocità 
critica che comporta l’instabilità dinamica del sistema: 


—4MVyWy — 4c3 Jovowe 
pD (Ree + 1) [caDCm,a(0) — Cy,a(0)] 


—Acimvywy — 4Jolgwg 
pD(R+c1)[-c1Cy,a(0) + DCu,a(0)] 


dove cı e c2 sono due costanti, che derivano dalla decomposizione modale, date dalla 
seguenti relazioni: 


Sn? 
= EL (14.37) 
m (12 am wi) 
2 
es Silva (14.38) 


~ Jo (we — w) 


con wj e wa frequenze angolari proprie del sistema. Si può subito notare che se S} = 
0, allora cy = c2 = 0 ed il sistema risulta disaccoppiato in partenza. In quest’ultimo 
caso la velocità critica è la minima tra la velocità di Galloping traslazionale e quella 
di Galloping rotazionale, viste precedentemente. 


14.4 Flutter 


Il Flutter è un’instabilità dinamica originariamente studiata per i profili alari, e suc- 
cessivamente estesa alle strutture ed in particolar modo ai ponti di grande luce. Il 
termine Flutter, in realtà, riassume diverse tipologie d’instabilità dinamica: Flutter 
classico, Stall-Flutter, Flutter ad un grado di libertà. 

Il Flutter classico è un fenomeno cui corrisponde l’accoppiamento di due gradi 
di libertà della struttura, generalmente quello traslazionale e rotazionale, che innesca 
delle oscillazioni instabili. 

Lo Stall-Flutter è un’instabilità legata al moto torsionale della sezione, e indotta 
da una forzante di lift non lineare in prossimità dello stallo (perdita di portanza). 

Il Flutter ad un grado di libertà, include lo Stall-Flutter, ma in generale è associato 
a quelle sezioni la cui forma è tale da indurre una pronunciata separazione del flusso 


194 LEZIONI DI INGEGNERIA DEL VENTO 


in corrispondenza di essa. Sezioni tozze, poco aerodinamiche sono tipici esempi di 
sezioni strutturali candidate a questo tipo di Flutter. 

E doveroso ricordare che, in generale, il Flutter, in tutte le sue forme, coinvol- 
ge fenomeni aerodinamici non lineari. Tuttavia, sotto opportune ipotesi, è possibile 
linearizzare le forzanti indotte dall’interazione fluido-struttura e dunque studiare la 
soglia di stabilità della struttura. Lo studio dell’evoluzione dell’aerodinamica non 
stazionaria da cui trae origine il fenomeno, è rimandata a studi specifici generalmen- 
te più complessi. 

Le ipotesi che permettono lo studio del Flutter delle strutture sono fondamen- 
talmente due: la struttura è elastica lineare smorzata linearmente, e dunque presenta 
delle oscillazioni armoniche smorzate esponenzialmente; le oscillazioni sono piccole, 
e si ritiene che la soglia di stabilità sia raggiunta quando tali oscillazioni raggiungono 
un valore critico corrispondente alla velocità di Flutter. Pertanto, come per il Gal- 
loping e la Divergenza, l’analisi si riduce allo studio della velocità critica che porta 
all’instabilità della struttura, in questo caso dovuta al Flutter. La progettazione ae- 
rodinamica della struttura sarà tanto più affidabile quanto più la velocità critica di 
Flutter sarà alta, o comunque maggiore della velocità eolica di progetto, avente una 
certa probabilità di eccedenza nel sito dove sorgerà la struttura. 


Figura 14.7: Convenzione dei segni per una generica sezione 


14.4.1 Equazioni di moto su un impalcato da ponte 


In Fig. 14.7 è schematizzata la sezione di un ponte, avente due gradi di libertà: la 
traslazione verticale, h, e la torsione, a. L'equazione di moto sarà dunque scritta nel 
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seguente modo: 


myh + Sd + cnh +Crh = Lr 


* (14.39) 

Jo + Srh + Cad + Coa = Ma 
dove S, è il momento statico della sezione calcolato rispetto all’ asse x rivolto nella 
direzione della velocità media wy = U, Ja è il momento polare della sezione, my 
è la massa della sezione, cp, e Ca sono i coefficienti di smorzamento per i rispettivi 
gradi di libertà, mentre Ch e Ca sono rispettivamente le rigidezze alla traslazione e 
alla rotazione della sezione. 


14.4.2 Forzanti aeroelastiche L, e Ma per i profili alari sottili 


Come accennato precedentemente, il Flutter deriva dalle pratiche aeronautiche. Per- 
tanto, si accenna brevemente l’origine del problema facendo riferimento alla non 
stazionarietà delle forze aerodinamiche agenti su un profilo alare sottile in moto tra- 
slazionale e rotazionale. 

Un sezione di forma generica è tale da modificare le linee di flusso del corrente 
fluida incidente, creando dei disturbi distribuiti lungo la sezione stessa, e addirittura 
inducendo la separazione del flusso in corrispondenza di essa fino ad innescare la for- 
mazione di vortici che defluiscono a valle dando vita ad una scia vorticosa. Pertanto, 
quello che all’origine era un flusso irrotazionale, diviene un flusso rotazionale in cor- 
rispondenza della sezione, e conseguentemente la velocità delle particelle fluide, in 
prossimità della sezione, non è più tangente alla sezione stessa in quanto si genera 
una componente verticale della velocità. 

Nel caso di un profilo alare stazionario, i disturbi indotti dalla sezione stessa, so- 
no così contenuti che la velocità delle particelle fluide segue il contorno della sezione 
(sezione streamline), e dunque la portanza che si sviluppa è essenzialmente una forza 
media con una deviazione standard molto contenuta dal momento che i disturbi del 
flusso prossimo alla superficie, sono a loro volta molto contenuti. 

Se il profilo alare oscilla traslazionalmente e rotazionalmente, le particelle flui- 
de, in corrispondenza della sezione, presentano anche una componente normale alla 
superficie, sinonimo del fatto che la non stazionarietà del flusso, in questo caso, gio- 
ca un ruolo fondamentale inducendo delle forzanti aeroelastiche (dipendenti cioè dal 
moto della sezione) che allo stato limite ultimo inducono un’instabilità dinamica, no- 
ta appunto come Flutter. 

Theodorsen [17], utilizzando la teoria del flusso a potenziale, dimostrò che le 
forzanti Ly e Mp sono lineari in he ae nelle loro rispettive derivate del primo e 
secondo ordine. In particolare, Theodorsen arrivò alla definizione di una funzione, 
nota appunto come funzione di Theodorsen, che tiene in conto della non stazionarietà 
del flusso in prossimità del profilo alare, indotta dal moto della sezione stessa. Tale 
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LU. tr 
nb ok 


Figura 14.8: Funzione di Theodorsen [17], C(k) = F(k) + iG(k), in funzione della 


velocità ridotta U,;a = #2 


funzione, indicata con C'(k), è una funzione complessa: 
C(k) = F(k) + iG(k) (14.40) 
espressa in funzione della frequenza ridotta 


k= and (14.41) 


con b = È, essendo B la larghezza della sezione alare. In particolare, Ly e Mp, per 
un profilo alare, si esprimono nel seguente modo: 


a : 1 
Ln = —pb? (Urà +rh-— hati) ~2npC(k) [vo +h+b G z a) al (14.42) 


Ma = — pb? f G — a) Ubà + Tb? G l è) a arbh} | 


+ 2pUb? x G + a) C(k) [va +h+b G — a) al (14.43) 


dove ab è la distanza di metà sezione dal centro di rotazione (Fig. 14.7). 


14.4.3 Forzanti L, e Ma per impalcati da ponte 


Scanlan and Tomko [15] estesero lo studio di Theodorsen [17] alle sezioni da ponte, 
generalmente aventi spessori non trascurabili. Ciò che ne derivò furono le cosiddet- 
te derivate aeroelastiche o derivate di Flutter (cfr. Scanlan and Tomko [15], Simiu 
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and Scanlan [16]). La differenza sostanziale con la trattazione originaria di Theo- 
dorsen, risiede nel fatto che mentre la funzione di Theodorsen è una funzione nota 
a priori, le derivate aeroelastiche devono essere necessariamente stimate tramite spe- 
rimentazione in galleria del vento, perché rappresentano l’estensione di una teoria 
esatta fondata su determinate ipotesi che cadono nel caso di sezioni come quelle dei 
ponti. Basti pensare che la teoria di Theodorsen si basa sull’ipotesi di profili alari 
sottili assimilabili cioè ad una lastra piana, ipotesi che permette di ipotizzare che i 
disturbi indotti dal moto della sezione, e da cui trae origine la non stazionarietà del 
flusso superficiale, possano essere intesi come uniformemente distribuiti sull’intera 
sezione stessa. Tramite un espediente geniale, ossia tramite una trasformazione con- 
forme, che trasforma la lastra piana in una circonferenza, Theodorsen riuscì a rendere 
la sezione alare una sezione streamline anche del flusso non stazionario. Grazie a tale 
espediente, congiuntamente all’ipotesi di non stazionarietà uniformemente distribui- 
ta, Theodorsen poté applicare i fondamentali concetti di base della fluidodinamica e 
quindi trovare l’espressione delle forzanti aeroelstiche di lift e di momento che ab- 
biamo visto precedentemente. 

È ben evidente che gli impalcati da ponte hanno sezioni con spessori tutt’ altro che 
trascurabili, responsabili dei disturbi del flusso che non solo non possono essere inte- 
si come uniformemente distribuiti (basta pensare al distacco del flusso in prossimità 
di uno spigolo a monte della sezione, seguito da un riattacco a valle della sezione, 
oppure che si propaghi nella scia senza nemmeno riattaccarsi alla sezione stessa), ma 
variano anche tra intradosso e extradosso della sezione stessa. 

Per questo le derivate aeroelastiche andrebbero usate e, ancor prima, ricavate con 
molta cautela. Si ricorda, infatti, che le forzanti espresse in funzione delle derivate 
aeroelastiche sono lineari e poiché soltanto per valori di k — 0 tali forzanti lineari as- 
sumono il loro valore quasi-stazionario (nel caso delle sezioni alari sottili C (k) — 1 
per k — 0, come di vede dalla Fig. 14.8), allora qualora, durante la sperimentazione, 
si lavori con valori di k >> 0 bisognerebbe porre notevole attenzione alle condizioni 
iniziali impresse al sistema: ampiezze iniziali di oscillazione relativamente grandi, 
per valori di k >> 0, inducono una non stazionarietà che si traduce in forzanti tut- 
taltro che lineari. Pertanto, le derivate aeroelastiche rischierebbero di non essere 
univoche, ma dipendenti dall’ampiezza di oscillazione iniziale, visto che per defini- 
zione l’interazione aeroelastica tra fluido e struttura si instaura per effetto della non 
stazionarietà indotta dal moto della struttura. 

Le forzanti aeroelastiche proposte da Scanlan and Tomko [15], e largamente 
utilizzate oggigiorno, sono le seguenti: 


Ln 


ae: ee: een ey, 
TUB KEO FKE r TERO FER 


B 
(14.44) 
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Ma duro RESTO: IPER ETRE. 
TOR KAIK) + KAK) y + KAR Ka + KOA(K) = 
(14.45) 


dove H¥(K) e A*(K) con j = 1...4 sono le derivate aeroelastiche, mentre la fre- 
quenza ridotta in questo caso viene indicata con K e assume la seguente espressione: 


2rf B 


K = 
U 


(14.46) 


con f frequenza di oscillazione, e B larghezza dell’impalcato. 
Nel caso di lastre piane (a = 0), la corrispondenza tra le derivate aeroelastiche e 
la funzione di Theodorsen è la seguente: 


Hi (K) = a (14.47) 
At(K) = -7 ch) (14.48) 
Hs(K)= T h + F(k) 4 Pa] (14.49) 
*(K) 7 |! F(k) sui (14.50) 
Hi(K) = re FO) ai (14.51) 
3(K) = go FO) sai (14.52) 
Hi(K) = 5 h + aoa (14.53) 
i(K) = nett) (14.54) 
k= = (14.55) 
dove: 

F(k) = oe - co E (14.56) 
GET (14.57) 


con J; e Y; funzioni di Bessel del primo e del secondo tipo, di ordine j. 
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Figura 14.9: Esempi di derivate aeroelastiche per varie sezioni 
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Figura 14.10: Esempi di derivate aeroelastiche per varie sezioni 
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Figura 14.11: Derivate aeroelastiche per il Great Belt Bridge (curve 1,2 e +). —.— 
Ponte a sezione reticolare. — Lastra piana 
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14.4.4 Soluzione per il calcolo della velocità critica di Flutter classico 


Si considerano le seguenti soluzioni del problema aeroelastico: 


gt 14. 
Fg : (14.58) 
a ee (14.59) 


Andando a sostituire tali soluzione nelle equazioni di moto del sistema (Eq. 
14.75), si ottengono le seguenti equazioni: 


$ B? 5 * x h 
[1 + 2iv, Kn K + K? — Z (ik? HI (K) + K?HI(K)| ca 


B? B? 
= [Z irem) n PE em (| œo =0 (14.60) 


4 
|- ereak) + K°A;(8)| do, 


a 
2 ; i eae 
+|-K + 2i0K Ka + Ka- ik A3(K)- 
a 


B4 
z PK A;(K) Oo =0 (14.61) 


dove K, = Bur, Ka, = Boo, Vp, € Va SONO rispettivamente i rapporti di smorzamento 
associati ai due gradi di libertà, h e a. 


Definendo 
w 

X = — (14.62) 
Wh 


e imponendo che il determinante dell’equazioni 14.60 e 14.61 sia nullo, si ottiene 
un’equazione del quarto ordine complessa. Considerando l’equazione reale e quella 
immaginaria separatamente, e risolvendole per diversi valori di K, è possibile dia- 
grammare le loro radici, X, in funzione di K. Il punto in cui le due curve coincidono 
indica la frequenza ridotta dalla quale è possibile ricavare la velocità critica di Flutter 
classico. 


14.4.5 Flutter ad un grado di libertà 


La frequenza di oscillazione traslazionale e torsionale dipendono dalla dalle derivate 
aeroelastiche H} (K) e Aj(K) rispettivamente. Poiché tali derivate aeroelstiche di- 
pendono a loro volta dalla frequenza ridotta, K, allora l effettiva frequenza di oscil- 
lazione dovrebbe essere calcolata iterativamente. Talvolta è preferibile evitare tale 
calcolo iterativo, e ciò è possibile come vedremo nei successivi paragrafi. 
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Flutter traslazionale 


x 


In questo caso il moto è puramente traslazionale, pertanto le derivate H3() e 
H3(K) sono nulle. La rigidezza totale del sistema sarà dunque data dalla somma 
della rigidezza strutturale e di quella aerodinamica: 


1 
Cha = Ch- zP’ K 2 Hi(K) (14.63) 


da cui si evince che valori negativi di Hï (K) incrementano apparentemente la ri- 
gidezza della struttura. In particolare la frequenza propria traslazionale sarà data 
da: 


2 
pB? tit pB? 
satin 14.64 
fht in| Imp PP i( ) th Imp, i( ) (14.64) 
dove 

2r fn ¿B 

Kae = iB (14.65) 
U 

2r fn B 

K= mh (14.66) 


Analogamente, lo smorzamento totale è somma di quello strutturale e di quello 
aerodinamico: 


th pB? 
=A= ASHK 14.67 
Vht fi Amy, i( ae) ( 67) 


Quando lo smorzamento totale è nullo allora si può trovare il valore di Ht (Kae) 
che innesca il Flutter traslazionale corrispondente all’instabilità dinamica dovuta al 
solo moto traslazionale. Il valore della frequenza ridotta associata al valore critico di 
H¥ (Kae), permette di trovare la velocità critica. Tuttavia per ponti di grandi luce che 
si muovono di solo moto traslazionale, lo smorzamento totale risulta sempre positivo 
perché Hy (Kae) risulta sempre negativo (Figg. 14.9, 14.10 e 14.11). Pertanto tale 
instabilità è poco ricorrente. 


Flutter rotazionale 
Nel caso di moto puramente torsionale si ottiene la seguente rigidezza totale: 
Cat = Ca — =pU?B?K?AY(K 14 
at = Ca — 5pU 3(K) (14.68) 


Pertanto la frequenza naturale rotazionale della struttura verra scritta nel seguente 
modo: 


B4 f? 
fat = taj = COTE AY Koy) = fof 1 — PY As(K) (14.69) 
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Analogamente lo smorzamento totale sara dato dalla somma dello smorzamento 
strutturale e di quello aerodinamco: 


fa pB* * 
SLI <A. 14.70 
Uh,t pe: Ado, 5( ) ( ) 


Quando lo smorzamento aerodinamico risulta nullo, si ha il Flutter torsionale (insta- 
bilità torsionale) corrispondente ad un valore di A3( Kae) pari a: 


Ja , fa 
pB* fag 


’ 


A5(Kae) = 4 (14.71) 
da cui è possibile trovare il valore critico della velocità del vento che causa l’instabi- 
lità nota la frequenza ridotta per la quale Až (Kae) è critico. Questo tipo di instabilità 


non è rara per i ponti al contrario del Flutter traslazionale (instabilità traslazionale) 
(Figg. 14.9, 14.10 e 14.11). 


14.5 Fenomeni aeroelastici in presenza della turbolenza at- 
mosferica 


La risposta strutturale dovuta alla componente turbolenta della velocità del vento è 
stata ampiamente discussa nei capitoli precedenti. In questa sezione, consideriamo i 
fenomeni aeroelastici associati ai ponti in presenza delle componenti fluttuanti della 
velocità eolica, u = u(t) e w = w(t), relative rispettivamente alla turbolenza longi- 
tudinale (lungo x) e verticale (lungo 2). 

Le forzanti indotte dalla turbolenza atmosferica sono note come forzanti di buffe- 
ting, e possono essere scritte, con l’ausilio della teoria quasi-stazionaria linearizzata, 
nel seguente modo: 


u 

U 
nn poe 

M = 5 PU B° |\2CM(Qo0)XM,u 


1 
Ly = 5pU"B [2Cr (0) Xz. + (Cha(de) + Cp(eo)) XL | (14.72) 


u 


W 
ut Cuala) Xo 7 | (14.73) 


U 

ottenute linearizzando le forzanti intorno all’angolo d’attacco medio a. Come sem- 
pre vengono trascurate i termini non lineari delle componenti fluttuanti della velocità 
del vento: 


(U +u)? = U? +2Uu (14.74) 


Inoltre, passando alla densità spettrali, i termini x? indicano, come al solito, le am- 
mettenze aerodinamiche delle varie componenti delle forzanti di buffetting. 

Le componenti turbolente della velocita del vento inducono un moto della strut- 
tura. Nascono, quindi, delle forzanti aeroelastiche esprimibili con le equazioni 15.9 
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e 15.10. Pertanto, le equazioni di moto, che risolvono la risposta strutturale della 
struttura, saranno le seguenti: 


mph + Sei + cnh + Crh = Ln + Lo 


; (14.75) 
Tod + Sh + Co + Coa = Ma + My 


Un approccio immediato, per il calcolo della risposta strutturale in presenza di 
turbolenza atmosferica, consiste nel trascurare le derivate aeroelastiche che accop- 
piano aerodinamicamente i due gradi di libertà della struttura, ossia H5(K), Hj(£), 
Aj(K) e Aj(K), dal momento che la risposta strutturale è valutata in un intervallo di 
velocità del vento sufficientemente lontano dalla velocità critica che innesca il Flutter 
accoppiato. Pertanto, attraverso una decomposizione modale è possibile disaccoppia- 
re i due gradi di libertà anche meccanicamente. Si ottengono dunque due sistemi ad 
un solo grado di libertà, la cui risposta può essere calcolata con le procedure discusse 
nei capitoli precedenti per la valutazione della risposta strutturale alle azioni turbo- 
lente indotte dal vento. 

Tuttavia, per velocità del vento sufficientemente vicine alla velocità di Flutter ac- 
coppiato, non è più possibile disaccoppiare aerodinamicamente i due gradi di libertà 
(Dyrbye and Hansen [4], Simiu and Scanlan [16]). La trattazione non è qui riportata, 
ma dalla Fig. 14.12 si può infatti notare che la risposta del sistema aerodinamica- 
mente accoppiato e di quello disaccoppiato coincidono solo fintanto che la velocità 
eolica non eccede un valore circa pari al 60% della velocità critica di Flutter. 


RMS, vertical dellection 


RMS, torsional rotation 


100 


50 
Wind velocity, U|m/s] 


Figura 14.12: — Sistema aerodinamicamente accoppiato; —— Sistema aerodinami- 
camente disaccoppiato 


Risposta aeroelastica delle strutture secondo 
l’Eurocodice 1 


L’Eurocodice-1 [5] fornisce delle regole per la progettazione di strutture sensibili alle 
instabilità aeroelastiche. Come vedremo, tali regole si fondano sulla trattazione ge- 
nerale affrontata nel Cap. 14. Pertanto è possibile riconoscere facilmente le formule 
trovate precedentemente. 


15.1 Galloping 


La velocità critica di Galloping, secondo L’ECI, è data dalla formula: 
25 
Uger = ——N1,yb (15.1) 
ag 


dove b è la larghezza della sezione trasversale, nı y è la frequenza di vibrazione della 
struttura nella direzione ortogonale a quella del vento, ag è il fattore di Galloping 
riportato in Fig. 15.1, pari a 10 se non specificato, e Se è il numero di Scruton 
definito come segue: 

_ 4nume 


= TE (15.2) 


dove me è la massa equivalente per unità di lunghezza: 


do SLY (Am(2)dz 
e J, v2 (z)dz 


(15.3) 
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con il consueto significato dei simboli. 
Inoltre deve essere garantito che: 


Uger > 1.25Um (15.4) 


dove Um è la velocità media del vento calcolata in corrispondenza dell’ altezza in cui 
si innesca il Galloping, ritenuta come punto in cui si ha la massima ampiezza di oscil- 
lazione. 

Quando la velocità critica di sincronizzazione, Us, er, è prossima a quella di Gal- 
loping, Ug,er, ossia quando 


U, cr 
0.7 < Ta < 1.5 (15.5) 


è probabile che s’inneschi un interazione tra il distacco dei vortici e il Galloping. Si 
raccomandano pertanto prove in galleria del vento o consulenze specialistiche. 
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Figura 15.1: Fattore d’instabilità di Galloping a, 
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0 0,05 0,10 0,15 0,20 0,25 


Figura 15.2: Cw 0 per sezioni rettangolari 


15.2 Divergenza e Flutter 


Una generica sezione strutturale è incline ad instabilità aeroelastiche, di tipo Diver- 
genza o Flutter, se i seguenti criteri sono tutti contemporaneamente soddisfatti: 


e La struttura, o una parte sostanziale di essa, deve avere una sezione trasversale 


allungata con proporzioni tra larghezza, d, ed altezza, b, maggiore o uguale di 
4: d/b > 4. 


e L'asse torsionale deve essere parallelo al piano medio della sezione e normale 
alla direzione del vento, ed il centro di torsione deve essere almeno d/4 dal 
bordo sopravento della sezione. 


e La frequenza naturale minore deve corrispondere ad un moto torsionale, oppu- 
re la minima frequenza torsionale deve essere minore di due volte la più bassa 
frequenza traslazionale. 
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15.2.1 Divergenza 


2ko 


reg e 
si pdC mo 


(15.6) 


dove kg è la rigidezza torsionale, e Cm o è la derivata prima del coefficiente aerodi- 
namico di momento, Cm, rispetto all’angolo d’attacco 0. 
Deve essere garantito che: 


Uder 20) (15.7) 


dove Um(Zeq) è la velocità eolica media misurata alla quota equivalente, Zeq, definita 
in Fig. 11.7. 
Alcuni valori di Cm per sezioni rettangolari sono definiti in Fig. 15.2. 


15.2.2 Flutter 


Per quanto riguarda il Flutter, L’EC1 rimanda alla soluzione del problema mediante 
prove in galleria del vento. 

Sebbene non direttamente specificato, nella pratica si fa spesso uso del modello 
di Scanlan and Tomko [15] che viene nuovamente riportato: 


h+ S,à + ch + Cyh=L 
a Solo A (15.8) 


Ja + Sh + Co + Caa = Ma 


con Ly e Ma definite nel seguente modo: 


Ln oah apn BA oora forati 
E KH(K) + KHK) +K H3(K)a + K°Hi(K)z| (15.9) 
Ma * h Pi Bà d: Sus h 
(15.10) 


dove le derivate aeroelstiche, H} (K) e A}( K) con j = 1...4, sono stimate mediante 
prove in galleria del vento. La valutazione della velocità critica di Flutter avviene 
mediante la procedura descritta nel paragrafo 14.4.4. 


Sperimentazione in galleria del vento 


Sebbene l’impiego della fluidodinamica computazionale sia in rapida espansione, og- 
gigiorno rimane la necessità di effettuare esperimenti fisici per comprendere i com- 
plessi fenomeni legati ad una corrente fluida che investe un corpo tozzo. Per questo 
motivo, le gallerie del vento hanno avuto, e continuano ad avere, un notevole svilup- 
po non solo nel campo dell’aeronautica, dove furono utilizzate per la prima volta, ma 
anche nel campo dell’ingegneria strutturale. Nelle simulazioni in galleria del ven- 
to, la sperimentazione viene effettuata su modelli in scala per ovvie considerazioni 
di carattere economico e di convenienza. La scalatura del modello apre le porte ad 
una vasta gamma di teorie che si fondano su parametri adimensionali e/o su crite- 
ri di similitudine applicabili alla riproduzione del flusso e alla sperimentazione sui 
modelli. 


16.1 Criteri di similitudine 


Per analizzare il problema, conviene ragionare in termini di parametri adimensiona- 
li. Tali parametri possono essere trovati impostando le equazioni differenziali che 
descrivono il sistema fisico in questione. Successivamente, vengono divisi per quan- 
tità che hanno la medesima dimensione, così da ottenere delle quantità adimensionali 
che governano il comportamento del sistema. Mantenendo invariato il valore di al- 
cuni di questi parametri, nel passaggio da prototipo (struttura reale) a modello, viene 
automaticamente garantita la similitudine. 
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16.1.1 Analisi dimensionale 


La forza F agente su un corpo immerso in un fluido è funzione di sei parametri: la 
densità del fluido, p, la velocità del flusso, U, una dimensione caratteristica del corpo, 
D, una frequenza n, la viscosità del flusso, u, e l’accelerazione gravitazionale, g. In 
particolare si ha: 


F = p°u? D nd u‘ gS (16.1) 
Applicando il teorema 7 si ottiene la seguente espressione per la forzante: 
ô € È 
D u gD 
F = pU?D? (È 16.2 
i ( Uj) (UD) (U? Foa 


da cui si evince che il coefficiente di forza ID e funzione delle seguenti quantita 


adimensionali: 


| a nD u gD 
E E us 


In particolare si definiscono i seguenti numeri: 


vsD 
S= “= (16.4) 
ie me — “2 (16.5) 
U2 


noti rispettivamente come numero di Strouhal, numero di Reynolds, e numero di 
Froude. 


16.1.2 Considerazioni di base per la scalatura 


Per una corretta simulazione in galleria del vento è necessario che si conservino i 
valori dei numeri appena menzionati, ovvero, nel passaggio da prototipo a modello è 
necessario che si verifichino le seguenti uguaglianze: 


D D UD UD U? U? 
dadi an) Ga), 
U J; U > V J; v jy gD /, gD . 
dove i pedici r ed g si riferiscono rispettivamente alla struttura reale e al modello 
sperimentale. 


Si introducono i cosiddetti fattori di scala, Ar, Au, € An, noti rispettivamente 
come scala delle lunghezze, scala delle velocità, e scala delle frequenze: 


Au 


(16.8) 
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cosicché le relazioni 16.7 possono essere riscritte in funzione dei fattori di scala: 


AU 1 
AL = i Ar = = =e (16.9) 


ottenute tenendo presente che le quantita p, v e g rimangono invariate nel passare 
dalla struttura reale al modello di galleria. Le relazioni riportate nella 16.9 sono note 
rispettivamente come condizione di Strouhal, condizione di Reynolds, e condizione 
di Froude. 

Particolare attenzione va rivolta alla scalatura del numero di Reynolds: se il proto- 
tipo ed il modello sono entrambi investiti da una corrente d’aria, allora per conservare 
il numero di Reynolds è necessario che Ay/AL = 1, in accordo con la condizione di 
Reynolds. Tuttavia, questa condizione genera delle discordanze con gli altri requi- 
siti necessari per una corretta scaltura, come si può notare dalle altre due relazioni 
della 16.9. Segue che la simulazione del numero di Reynolds è incompatibile con i 
requisiti prioritari di scalatura delle lunghezze e delle velocità, a meno che le prove 
vengano effettuate con A; = 1, condizione ovviamente irrealizzabile nel caso delle 
strutture. 

In particolare, si può osservare che se si usasse una scala delle lunghezze Az = 
1 : 100, in galleria del vento occorrerebbe una velocità del flusso 100 volte maggiore 
della realtà. Per contro, la condizione di Froude, impone che la velocità reale sia sca- 
lata di VAL = 1/10. Pertanto le due condizioni sono discordanti. Inoltre, anche se 
fosse ipoteticamente rispettata la condizione Ay = 1, la condizione di Strouhal vor- 
rebbe che A, = 100 ossia per simulare i fenomeni in un’intervallo di frequenza tra 
0 e 4Hz (intervallo rappresentativo della turbolenza atmosferica reale), bisognerebbe 
riprodurre in galleria del vento uno spettro compreso tra 0 e 400Hz, cui corrisponde, 
per il teorema di Nyquist, una frequenza di acquisizione almeno pari a 800Hz. 

Nelle pratiche strutturali, la condizione di Reynolds risulta essere sicuramente 
molto importante, ma non risulta quasi mai soddisfatta, sebbene talvolta non è stretta- 
mente necessaria. Il numero di Reynolds delinea, come visto nei precedenti capitoli, 
diversi regimi in corrispondenza dei quali, in generale, la struttura ha un comporta- 
mento aerodinamico differente. Basti pensare al distacco dei vortici: ad esempio, nel 
regime critico non si osserva distacco dei vortici per cilindri circolari, pertanto questa 
condizione dovrebbe essere rispettata passando dalla struttura reale alla sperimenta- 
zione, altrimenti si osserverebbero fenomeni completamente diversi dalla realtà. 

Generalmente, si utilizza una rugosità tecnica da applicare al modello, in modo 
da riprodurre i fenomeni che si avverrebbero per numeri di Reynolds più elevati. 

L’aerodinamica dei corpi tozzi aventi spigoli vivi è governata dalla geometria 
e risulta essere indipendente dal numero di Reynolds, pertanto in tale situazione è 
possibile trascurare la condizione di Reynolds. Quando si hanno basse accelerazioni, 
come nel caso di simulazione di fenomeni ambientali, è importante soddisfare la con- 
dizione di Froude, che tuttavia è inconciliabile con quella di Reynolds. Infine, per la 


216 LEZIONI DI INGEGNERIA DEL VENTO 


Figura 16.1: Modello di una spire 


sperimentazione aeroelastica risulta essere molto importante la corretta riproduzione 
delle frequenze e degli smorzamenti strutturali. 


16.1.3 Simulazione della circolazione atmosferica 


Il rispetto dei criteri di similitudine appena visti, garantisce una fedele simulazione 
delle caratteristiche del flusso incidente sui prototipi, ossia del profilo di velocità me- 
dia, della turbolenza e delle scale integrali, e dunque degli spettri e dei cross-spettri 
della componente fluttuante della velocità del vento. 

Particolare attenzione viene dedicata alla simulazione dello strato limite atmosfe- 
rico del sito in esame. La generazione dello strato limite avviene disponendo, lungo 
tutto lo sviluppo longitudinale della galleria del vento, dei pannelli di rugosità (rugo- 
sità superficiale) che simulano la rugosità del sito in questione. Tuttavia, affinché il 
flusso raggiunga le caratteristiche desiderate, occorrerebbero delle gallerie molto lun- 
ghe, requisito non sempre rispettato. Sorgono, pertanto, due importanti questioni: la 
limitata lunghezza della galleria porta ad avere da un lato uno strato limite di modesta 
altezza (si ritiene che l’altezza debba essere almeno 2-3 volte più grande dell’altez- 
za del modello), dall’altro una turbolenza costituita da piccoli vortici che alterano lo 
spettro reale della componente fluttuante della velocità del vento. Per questo motivo, 
oltre alla rugosità superficiale, si pongono, nella sezione d’imbocco della galleria, dei 
dispositivi che alzano lo strato limite che si svilupperebbe con il solo attrito con la 
rugosità superficiale, e capaci di generare scale integrali di turbolenza di dimensioni 
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Figura 16.2: Curve per la determinazione della larghezza b delle spires; H è l’altezza 
della camera di prova della galleria del vento, Irwin [10] 


maggiori. Tra questi dispositivi abbiamo le griglie e le spires, disposte verticalmente, 
e le barriere orizzontali. 

Vediamo in dettaglio le spires (Fig. 16.1). Fissata l'altezza desiderata dello strato 
limite 6, nonché il profilo di velocità media del vento caratterizzato dall’esponente a, 
è possibile ottenere l'altezza h delle spires secondo la seguente relazione proposta da 
Irwin [10]: 


_ 1.396 
_ 1+a/2 


(16.10) 


Per ottenere la larghezza b delle spires si può fare riferimento alla Fig. 16.2, in cui 
con H è indicata l’altezza della camera di prova della galleria del vento. 

Il profilo di velocità desiderato viene raggiunto ad una distanza pari a 6h dalle 
spires. Come gia accennato, è necessario disporre lungo la dimensione longitudinale 
della galleria, dei pannelli di rugosità, costituiti da cubi di altezza k tale che sia 
soddisfatta la seguente relazione valida per 30 < 5D?/k? < 2000: 


k 2 D 2 1/2 
33 © (5) In (5) — 0.1161 (È) + 205) (16.11) 


218 LEZIONI DI INGEGNERIA DEL VENTO 


dove D è la distanza tra gli elementi di rugosita, mentre C's è definito come segue: 


2 
a 

Ce = 0.136 16.12 

f =] ( ) 


16.2 Simulazione aerodinamica ed aeroelastica dei corpi toz- 
zi 


In questa sezione viene discussa la dipendenza della risposta aerodinamica ed aeroe- 
lastica dei modelli testati in galleria del vento, dalle caratteristiche del flusso simulato 
e dal suo numero di Reynolds. Inoltre, verrà discusso l’effetto del bloccaggio. 


16.2.1 Effetti della turbolenza del flusso 


La struttura dei vortici e la loro dimensione può essere riprodotta fedelmente in gal- 
leria del vento. Come menzionato nei capitoli precedenti, la turbolenza atmosferica 
nasce dall’attrito con il terreno. Pertanto, se il terreno viene riprodotto nella medesi- 
ma scala del modello, è automaticamente soddisfatta la legge di Jensen che si scrive 
nel seguente modo: 


(=) -(=) (16.13) 
Zo F Zo g 


dove h è l’altezza della struttura e zp è la lunghezza di rugosita del terreno circostante. 
Analogamente al numero di Reynolds, Strouhal e Froude, si definisce il numero di 
Jensen: 


Je= — (16.14) 


Quando è soddisfatta la legge di Jensen, la struttura dei vortici e lo spettro della 
componente fluttuante della velocità del vento, vengono simulati fedelmente in gal- 
leria del vento. Se in luogo di zo si considera la scala integrale di turbolenza, L%,, 
si ottiene che anche i vortici all’interno dello strato limite devono rispettare la stes- 
sa scala utilizzata nella realizzazione del modello. Tuttavia, la conservazione della 
quantità h/L%, non è sempre possibile, anzi particolarmente difficile, vista l’impos- 
sibilità di ottenere scale di turbolenza sufficientemente grandi nelle gallerie del vento 
comunemente utilizzate. Generalmente, si ritiene che essa sia soddisfatta se si ottiene 
una buona corrispondenza con i dati reali a disposizione. 


16.2.2 Effetti del numero di Reynolds: rugosità tecnica 


Gli spigoli vivi di una struttura causano un immediato distacco del flusso, indipenden- 
te dal numero di Reynolds. Per questa ragione, si assume che le pressioni indotte su 
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Figura 16.3: Esempio di variazione del coefficiente di trascinamento (drag) con Re 


tali strutture siano correttamente riprodotte in galleria del vento se il flusso è adegua- 
tamente simulato. Tuttavia, nei corpi tozzi con un accentuato sviluppo longitudinale, 
si può riscontrare il riattacco del flusso, fenomeno che dipende dal numero di Rey- 
nolds. In queste condizioni, i valori delle forzanti indotte dal vento possono essere 
sensibilmente diversi. Inoltre, se la scala delle lunghezze è molto piccola, può acca- 
dere che il coefficiente di trascinamento sia notevolmente influenzato dal numero di 
Reynolds, dal momento che quest’ ultimo è notevolmente più piccolo di quello reale, 
e conseguentemente la forza di trascinamento è notevolmente più grande (Fig. 16.3). 
Nel caso di corpi costituiti da superfici curve, gli effetti del numero di Reynolds 
sono significativi, come si può vedere dalla Fig. 16.4 dove è riportato l’ andamento 
del coefficiente di drag per un cilindro circolare. 
Come visto nei paragrafi precedenti, una corretta riproduzione del numero di 
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Figura 16.4: Variazione del coefficiente di trascinamento (drag) con Re per un 
cilindro circolare 


Reynolds in galleria del vento, vorrebbe che Az = 1, requisito ovviamente irrea- 
lizzabile. Per riprodurre fedelmente i fenomeni relativi a numeri di Reynolds molto 
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grandi, é necessario introdurre la rugosita tecnica. In particolare, lasciando invariato 
il numero di Reynolds di galleria, che per le ragione viste è necessariamente minore 
di quello reale, è possibile cambiare regime (ad esempio da quello laminare a quello 
turbolento) applicando al modello delle costole di altezza k poste ad una distanza d. 
L'altezza k di tali costole deve rispettare le seguenti relazioni: 


k 
cha > 400 (16.15) 
V 
k z3 
D < 10 (16.16) 


dove U è la velocità media del flusso, v è la viscosità cinematica del flusso, e D è 
una dimensione caratteristica del modello. 


16.2.3 Effetti del bloccaggio 


Un corpo, posto nella sezione di prova di una galleria del vento, ostruisce parzial- 
mente il passaggio dell’aria, causando un'accelerazione del flusso, in corrispondenza 
delle regioni laterali del corpo stesso. Tale effetto è noto come bloccaggio. General- 
mente, si cerca di scegliere una scala del modello tale per cui gli effetti del bloccaggio 
siano trascurabili. Qualora ciò non fosse possibile, gli effetti del bloccaggio diverreb- 
bero considerevoli a tal punto che il flusso intorno al modello ed il comportamento 
aerodinamico del modello stesso non sarebbero più veritieri. Per correggere l’effet- 
to del bloccaggio sul coefficiente di trascinamento, è possibile utilizzare la seguente 
relazione: 


Cp 


CD. = TE KS/0 


(16.17) 
dove Cp, è il coefficiente di trascinamento corretto, C'p è il coefficiente di trascina- 
mento misurato in galleria del vento, S è l’area di riferimento per il calcolo di Cp, 
e Cp, e C'è l’area della sezione trasversale della galleria. Il rapporto S/C è definito 
come bloccaggio. Il coefficiente K può essere determinato soltanto in un limitato 
numero di situazioni. 
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